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Abstract 

A  theory  for  the  ducted  propeller  is  developed  which  is 
based  on  a  linearized  annular  airfoil  theory  and  a  lifting¬ 
line  propeller  theory*  nte  fluid  is  assumed  to  be  invlscid 
and  incompressible  and  the  free-stream  velocity  to  be  axisym- 
metrlc.  As  with  propeller  theory  it  is  not  possible  to  obtain 
a  solution  in  explicit  fora  so  a  process  of  iteration  is  used* 

The  flow  field  of  the  annular  airfoil  is  represented  by 
a  distribution  of  ring  vor  '.ces  and  ring  sources  on  a  cylinder 
and  where  necessary  a  trailing  vortex  system.  This  approach 
allows  the  airfoil  section  to  have  an  arbitrary  shape  although 
the  annular  airfoil  itaelf  is  assumed  to  he  axlsyiiinetrlc.  The 
ring  source  strength  is  shown  to  be  a  function  of  only  the  duct 
thickness  while  the  ring  vortex  strength  is  a  function  of  camber* 
thickness  and  the  radial  velocity  induced  on  the  cylinder  by 
the  pt'opeller  and  duct  trailing  vortex  system*  In  the  presence 
of  the  propeller  two  couoled  singular  integral  equations  are 
derived  for  the  vortex  strength  which  are  reduced  to  two  coup¬ 
led  Fredholm  equations  of  the  second  kind.  (If  the  propeller 
is  not  present  only  one  integral  equation  is  obtained.) 

The  flow  field  of  the  v  ropeller  is  represented  bv  a  lift¬ 
ing  line  and  a  helicoidal  trailing  vortex  system.  This  allows 
the  propeller  to  have  a  finite  number  cf  blaueci  and  an  arbitrary 
distribution  of  circulation.  By  this  approach  the  propeller 
problem  essentially  reduces  to  the  proneller  by  itself  with  the 
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V 


inclusion  of  velocity  components  from  the  duct  and  hub. 
Gonset^uently,  it  reduces  essentially  to  the  propeller  problem 
solved  by  Lerbs. 

The  hub  is  treated  by  slender  body  theory  which  allows 
it  to  have  an  arbitrary  axisyranetric  shape.  One  consequence 
of  using  this  theory  is  that  the  hub  induces  no  tangential 
velocities. 
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Notation 
duct  chord 

axial  dlatance  betwaan  leading  adga  of  duct 
and  propeller 

axial  dlatance  between  trailing  edge  of  duct 
and  propeller 

number  of  blades 

thrust  coefficient,  equations  (2.7-9)  and  (4.4-4) 

power  coefficient,  equation  (4.4-6) 

mean  line  rdlnate  of  the  duct  section  measured 
from  the  nwse-tall  line 

complete  elliptic  Integral  of  the  second  kind 

nondlmensionalized  circulation  distribution  of 
the  propeller,  equation  (4.3-11) 

(a/2R^)  chord-diameter  ratio  of  the  duct 

unit  vectors 

modified  Bessel  function  of  the  first  kind 

modified  Bessel  function  of  the  second  kind 

modulus  of  the  elliptic  integrals 

modified  Struve  function 

propeller  torque 

ring  source  strength 

hub  source  strength 

dvict  radius 

propeller  radius 

cylindrical  coordinates 

radius  at  which  the  vortex  Is  shed  from  the 
propeller  blade 
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s(s) 


X»  <P  I* 


Xo 
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ft 

A  . 
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half  thickness  ordinate  o?  the  duct  section 

thrust 

ship  speed 

axial  component  of  Induced  velocity 

local  axial  velocity 

radial  component  of  Induced  velocity 

tangential  component  of  induced  velocity 

local  wake  fraction 

rectangular  coordinates 

nondlmenslonallzed  cylindrical  coordinates 

radial  coordinate  nondlmenslonalised  by  the 
propeller  radius 

nondimenslonallzed  radius  at  which  a  vortex  is 
shed  from  the  propeller  blade 

axial  coordinate  nondlmenslonalised  by  the  duct 
chord 

axial  coordinate  nondlmenslonalised  by  the 
propeller  radius 

z  -  a^ 

angle  of  attack  of  a  duct  section 

Ideal  angle  of  attack  of  a  duct  section 

relative  angle  between  free-stream  velocity  and 
duct 

propeller  advance  angle 

duct  advance  angle 

pxvpeller  hydrodynamic  pitch  angle 

ring  vortex  strength 

circulation 


^  coordinate 

ideal  efficiency,  ecjuation  (4,4-7) 

^  advance  coefficient,  equation  (4.6-2) 

advance  coefficient,  equation  (4.4-3) 
mass  density  of  fluid 
stream  function 
angular  velocity 


Subscsipts 

d 

h 

P 

q 

r 

ir 

5? 


duct 

hub 

propeller 
ring  source 
ring  vortex 

trailing  vortex  aystepi  of  the  vortex  cylinder 


Note:  Many  functions  are  defined  in  th'.  text. 
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A.  THEORY  OF  THE  DUCTED  PROPELLER  WITH  A  FINITE  NUMBER  OF  BUDES 
1.  INTRODUCTION 

The  name  *'ducted  propeller"  means  a  propeller-annular  air¬ 
foil  combination  acting  as  a  propulsion  unit.  The  annular  air¬ 
foil  or  duct  can  be  used  either  to  accelerate  the  flov  at  the 
propeller  (Kort  nozzle)^  or  to  deaccelerate  the  flov.  In  the 
first  type  of  flow  the  ducted  propeller  is  used  where  a  pro¬ 
peller  alone  would  be  heavily  loaded.  The  duct  accelerates 
the  flow  at  the  propellei  and  thus  the  propeller  operates  at 
a  more  favorable  loadings  in  addition,  the  duct  Itself  will 
generally  produce  a  positive  thrust.  In  the  case  of  the  duct 
which  deacceleratea  the  flow,  the  annular  airfoil  is  used  to 
increase  the  static  pressure  at  the  propeller  and  thus  delay 
cavitation  on  ship  propellers  or  decrease  compressibility 
effects  on  aircraft  propellers. 

Most  of  the  work  on  the  ducted  propeller  of  the  ’ICort 
nozzle"  type  has  been  done  in  Europe.  Well  known  are  the 
experimental  results  of  Van  lisnen  for  a  systematic 

series  of  "Kort  nozzles."  A  series  for  low-pitched  three- 
bladed  propellers  in  a  duct  has  also  been  given  by  Nakonechny^. 
Most  theoretical  approaches  have  been  restricted  to  representing 
the  duct  by  a  distribution  of  ring  vortices  along  a  cylinder^ 
ahH  the  enclosed  prcpcller  by  a  uiaLrlbutlon  of  sources  over 
the  propeller  disc. 

1  -  References  found  on  page  146 
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One  of  the  first  paoers  on  the  theory  of  the  ducted  pro- 

8  #7 

peller  is  that  of  Horn  in  which  he  uses  the  work  of  Oickmann 

for  representing  the  nozzle  and  considers  the  nropeller  as 

free-runni.ng.  This  orocedure  did  not  lead  to  a  design  method, 

Q 

In  1950  Horn  and  Aritsburg'  developed  a  design  piocedure  based 
on  representing  the  ^rict  by  a  vortex  distribution  and  the  pro¬ 
peller  by  a  sink  distribution.  Later  (1955)  Dickmann  and 
Weissingcr^®  considered  the  ducted  prcpcller  as  a  propulsion 
unit  and  in  the  optimum  t  se  represented  the  duct  by  ring 
vortices  with  effect  of  the  propeller  taken  into  consideration 
by  momentum  theory.  In  1959  Gutsche^^  developed  correction 
factors  from  the  propeller  in  a  long  tube  and  then  used  simple 
jet  tlieory, 

'?he  foregoing  described  results  have  in  general  been  re¬ 
stricted  to  ducts  which  accelerate  the  flow  at  the  propeller 
although  the  theory  developed  is  in  principle  applicable  to 
the  deaccelerating  duct,  Lerbs  has  applied  the  theory,  i.e, 
representing  the  duct  by  a  vortex  distribution  and  the  propeller 
by  a  sink  distribution,  specifically  to  deaccelerating  flow  in 
the  duct,  Kuchemann  and  Weber^^  have  considered  the  ducted  pro¬ 
peller  in  general  but  only  with  simple  momentum  theory. 

From  a  review  of  the  literature  it  is  apparent  that  the 

problem  of  the  ducted  propeller  with  a  fiin'ro  nrimK-j.  blades 

and  arbitrary  distribution  of  circulation  has  not  been  developed 

previously*  but  that  such  a  theory  is  necessary  for  the  ade«^ate 
♦After  the  work  in  this  report  was  completed,  Reference  [41] 
was  received  which  considers  a  lightly  load€',d  propeller  with  a 
finite  number  of  blades  in  a  duct  of  zero  tl.ickness. 
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design  of  propellers  operating  in  a  duct.  This  paper  presents 
such  a  theory  in  which  it  is  assumed  that  the  nozzle  flow 
field  can  be  represented  by  a  linearized  theory  and  the  pro¬ 
peller  by  the  lifting  line  theory.  As  usual  for  potential  flow 
problems  a  number  of  assumptions  are  made  about  the  fluid  as 
well  as  about  the  flow  field.  For  this  problem  these  are: 

1.  The  fluid  is  inviscid  and  incompressible  and  no  separa¬ 
tion  occurs  on  the  duct. 

2.  Body  forces,  such  a  gravity,  may  be  neglected. 

3.  The  free-stream  flow  is  axisymmetric  and  steady  with 
respect  to  a  coordinate  system  attached  to  the  propeller.  This 
allows  a  radial  variation  in  velocity  and  implies  that  the 
coordinate  system  is  rotating  v;ith  the  propeller.  It  causes  no 
loss  in  generality  to  assume  the  duct  is  also  rotating  since 
the  duct  by  itself  at  zero  incidence  induces  no  tangential 
velocity. 

4.  The  annular  airfoil  is  axisymmetric  and  of  finite  length, 

5.  The  thickness  and  camber-distribution  of  the  annular  air¬ 
foil  section  can  be  expanded  in  a  Fourier  series  with  resoect 
to  the  axial  coordinate.  This  assumotion  offers  no  restriction 
to  streamline  shapes. 

6.  The  linearized  flow  around  the  annular  airfoil  can  be 
represented  mathematically  by  a  distribution  of  ring  vortices 
and  ring  sources  along  a  cylinder  of  diameter  R^. 

7.  The  trailing  vortex  system  from  the  duct  has  a  constant 
diameter  Rjj  and  extends  from  the  duct  to  infinity. 
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8.  The  influence  of  all  induced  velocities  on  the  shape  of 
the  trailing  vortex  system  from  the  duct  is  neglected.  This 
implies  that  the  pitch  angle  of  all  the  free  vortex  lines  from 
the  duct  is  the  same  and  equal  to 

f^nAi  = 

9.  The  propeller  flow  field  can  be  represented  by  a  lifting 
line  and  trailing  vortices,  i.e,  a  horse-shoe  vortex  system. 

The  trailing  vortex  syst^'m  is  directed  along  helical  stream 
lines  trailing  aft  from  the  propeller  blades.  Each  vortex  is 
of  constant  pitch  and  lies  on  a  cylinder  of  constant  diameter. 
This  implies  that  the  contraction  of  the  slip-stream  is  ignored. 

Using  these  aasumptioiis  the  linearized  boundary  conditloiis 
on  the  duct  are  derived.  The  annular  airfoil  is  first  considered 
by  itself  and  its  flow  field  represented  by  ring  sources  and 
ring  vortices.  The  strength  of  each  is  chosen  so  that  the 
boundary  conditions  are  satisfied.  With  the  strength  of  the 
vortex  and  source  distribution  known,  the  entile  flow  field 
of  the  duct  can  be  derived.  The  propeller  with  a  finite  hub 
is  then  added  to  the  flow  field  and  the  interaction  effects 
determined.  The  solution  to  the  problem  of  the  combination  of 
the  propeller  and  duct  reduces  to  a  process  of  iteration. 
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II.  T.TNEARIZED  THEORY  OF  THE  ANNULAR  AIRFOIL 
II. 1  Previous  Theories 

The  theory  of  the  annular  airfoil  has  been  discussed  by 
numerous  investigators  but  in  most  cases  was  not  developed 
sufficiently  for  ducts  of  arbitrary  section  shape.  The  first 
theoretical  discussion  of  the  annular  airfoil  appears  to  be  that 
of  Dickmann^®  who  represented  the  annular  airfoil  in  uniform 
axial  flow  by  a  distribut  m  of  ring  vortices.  This  is  equiva¬ 
lent  in  thin  wing  theory  to  representing  an  airfoil  by  a  distri¬ 
bution  of  vortices  only, and,  thus,  the  thickness  of  the  foil  is 
neglected.  In  linearized  two-dimensional  wing  theory  neglecting 
the  thickness  is  justifiable  for  obtaining  the  lift  but  not  for 
the  pressure  and  velocity  distribution.  In  annular  airfoil 
theory  it  wuld  be  expected,  because  of  the  interference  effects, 
that  the  thickness  plays  a  more  important  role  than  in  linear¬ 
ized  wing  theory. 

The  solution  for  the  induced  velocities  from  a  single  vor¬ 
tex  ring  was  given  in  the  form  of  elliptic  integrals  by  Lamb.^^ 

By  arranging  ring  vortices  of  varying  strength  along  a  cylinder 
Dickmann  represented  the  resulting  integral  of  elliptic  integrals 
by  a  Fourier  series.  Tabulated  coefficients  for  determining 
the  velocity  distribution  of  a  vortex  ring  nnd  vertex  cylinder 
are  found  in  Kuchemann  and  Weber^^  and  a  more  complete  theo¬ 
retical  development  of  singularities  useful  for  this  problem  is 
found  in  a  report  by  Meyerhoff  and  Finkelstein^^. 
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Stpwart,^^  indenendently  of  Dickmann,  derived  the  induced 
velocities  of  an  annular  airfoil,  again  represented  by  a  vortex 
cylinder,  using  the  vector  potential.  He  was  able  to  represent 
the  velocity  comnonents  bv  integrals  of  a  product  of  modified 
Bessel  functions. 

Some  work  has  also  been  done  on  flow  about  thick  annular 
airfoils  and  foils  at  an  angle  of  attack.  Specifically 
Kuchemann^^  considered  anr  'lar  airfoils  of  finite  thiclcness 
wiHiout  circulation  by  a  distribution  of  source  and  sink  rings 
and  later  Kuchemann  and  Weber  considered  annular  foils  of 
finite  thickness  x<;ith  circulation  but  infinite  length.  In 
either  case  the  theory  was  not  adapted  to  foils  of  arbitrary 
shape  but  the  shape  and  velocitv  distribution  calculated  for 
an  assumed  distribution  of  sources  and  sinks. 

Weissinger^^  has  disc’issed  the  flow  field  about  annular 
airfoils  with  zero  thickness  operating  at  an  angle  of  incidence. 
To  represent  the  flow  mathematically  he  uses  a  distribution  of 
ring  vortices  along  the.  duct  whose  strength,  T  (^,z),  at  a 
point  on  the  ring  is  dependent  on  the  angular  position  as  well 
as  the  axial.  Since  in  this  case  there  are  free  vortices  in 


addition  to  the  bound  vortices,  he  uses,  in  addition,  a  system 
of  vortices  of  strength  trailing  from  the  cylinder.  An 

integral  equation  for  the  vortex  distribution  results  from  this 
analysis  which  is  solved  approximately,  Weissinger  later 


included  the  effect  of  finite  thickness  by  using  a  distribution 


of  ring  sources.  An  approach  which  follows  Weissinger’s 
work  very  closelv  is  tliat  of  Bagley,  Kirby  and  Marccr  , 

In  their  work  use  is  made  of  standard  ring  vortex  distri¬ 
butions  which  were  tabulated  bv  Kuchemann  a\id  Weber^^.  Since 
this  method  is  restricted  to  satisfying  the  boundary  condition 
at  an  arbitrary  number  of  noints  along  the  chord,  (maximum 
of  five)  it  gives  only  an  approximate  solution.  Their  consid¬ 
eration  of  the  annular  airfoil  at  on  angle  of  attack  is  similar 
to  Weissinger’s. 

22 

Recently  Pivko  considered  annular  airfoils  t^rith  thick 
symmetrical  sections  but  his  work  is  only  applicable  to  nozzle 
length-diameter  ratios  of  much  less  than  one.  This  restriction 
comes  from  the  fact  that  he  assumed  a  svmmetrical  section  could 
be  represented  hr'  a  distribution  of  sources  and  sinks.  Because 
of  the  interference  effects,  however,  a  vortex  distribution 
must  also  be  iised  in  addition  to  the  source-sink  distribution. 
Pivko  has  also  considered  thick  cambered  annular  airfoils 
operating  at  an  angle  of  attack  and  included  the  effect  of 
the  propeller  by  a  sink  disc.  In  general  he  makes  use  of  the 
velocity  coefficients  given  by  Kuchemann  and  Weber^^  and 


superimposes  the  velocity  of  each  effect.  This  theory  is  not 
readily  adaptable  to  sections  of  arbitrary  shape. 


t  «  ^  f‘ 

In  addirion  to  the  precedinj'  work  Halavard*'"'  has  considered 


the  pressure  distribution  on  annular  airfoils  with  and  without 

25 

thickness  using  electrical  analogy  and  Hacques  has  considered 
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the  problem  of  designing  for  a  given  pressure  distribution 
also  by  electrical  analogy. 

In  the  next  section  the  linearization  of  the  boundary 
conditions  on  the  duct  is  considered.  ^Tlien  linearizing 
inviscid,  incompressible  fluid  problems,  it  is  the  boundary 
conditions  which  are  linearized  as  the  equation  of  continuity^ 
which  reduces  to  Laplace's  equation  in  this  case,  is  linear. 

11.2  Linearized  Boundary  Conditions 
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The  coordinate  system  which  3  s  adopted  is  a  cylindrical 
system  (  ^  f  ,  O  with  the  zero  axial  coordinate  ( ^  )  located 
at  the  nropeller  blade  center  line  and  the  free-stream  flow 
is  from  right  to  left  (see  Figure  L).  The  annular  aiifoil  is 
assumed  axisymmetric  and  has  a  length  (a)»  Any  radius  on  the 
foil  can  be  chosen  as  the  reference  radius  (R<j)  but  since 
this  reference  radius  is  the  radius  of  the  cvlinder  along  which 
the  vortices  and  sources  are  distributed  it  would  seem.  i.iaSon- 
able.  to  use  some  sort  of  an  average.  For  convenience,  but  some¬ 
what  arbitrarily,  the  reference  radius  will  be  taken  as  the 
inside  rad3lus  of  the  annular  airfoil  at  "the  propeller. 

The.  section  shape  is  assumed  to  be  delineated  by  the  out¬ 
side  of  the  annular  airfoil  u(0  arid  by  the  inside  b(?)  as 
shown  in  figure  2. 


Figure  2.  The  annular  airfoil  section 
Using  the  above  notation  the  mean  line  of  the  foil  as 
measured  from  the  ^  -axis  (r=0)  is 


C(fJ  =5  -f  t(ejJ 


(2.2-1) 
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and  the  half  thickness  ordinate  is 

-b^Ol  (2.2-2) 

The  outer  surface  in  terms  of  the  mean  line  ordinate,  or 
camber,  o( ^  )  and  the  half  thickness  ordinate  s(f  )  is  then 
given  by 


y  =  ccCf)  =  c(()  see)  (2.2-3) 

and  the  inner  surface  by 

=  b(e)  ^  cC()  -  H()  (2.2-4) 

If  it  is  assumed  that  the  mean  line  deviates  little  from 
the  cylinder  of  radius  and  length  a  and  that  the  thickness 
is  small,  then  the  camber  and  the  thickness  can  be  expanded  in 
terms  of  a  nondimensional  perturbation  parameter  €. 

y.u.(?;e)--[c(f;e)  ^  se<;6)J  -  Rj*e[c^'U)^  s‘'>efj]  ^  ’ 

«  ;?J  f-  £  rc‘%)  y-  (2.2-5) 

e  =  b(fi  6)  =Cc  (f;  £)  - 6[c^'^(0  - 5 %)]  f  + - 

^  Rji-e[c ‘'^CO  -  s  "^CfJj  (2.2-6) 

As  f  -  0  the  problem  reduces  to  that  of  a  thin  circular  cylinder 
of  constant  diameter. 

Since  the  problem  is  axisymmetric  and  linear,  a  total 


stream  function  can  be  written  for  the  flow  in  terms  of  the 
perturbation  parameter  e  . 
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I 


+  -h  ■  (2.2-7) 

As  f-*0  the  stream  function  reduces  to  that  for  the  free- 
stream  velocity  alone.  From  equations  (2.2-5),  (2.2-6)  and 
(2.2-7)  the  linearized  boundary  conditions  are  developed. 

First  there  is  the  kinematic  boundary  condition,  i.e.  on  the 
surface  of  the  body  =  0.  As  a  first  approximation 

the  equation  for  the  outf  surface  of  the  ring  is  obtained 
from  equation  (2.2-5). 

fj  =  /?j  ^6[c^'k)  ^  -  »'  (2.2-8) 

and  for  the  inner  surface  of  the  ring  from  equation  (2.2-6) 
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and  ^  outside  of  the  annular  airfoil 


f  /  ^  r  1  -  _  / 

'Vvv*;^  -■“/ 


(2.2-12) 


and  on  the  inside 
=  -/ 


(2.2-13) 


With  these  values  for  :he  normal  to  the  surface,  it  follows 
from  equation  (2.2-10)  and  (2.2-11),  after  multiplying  through 
by  r  V^/Cr.f  )  +  F^^Cr,  I*  ),  that  the  kinematic  boundary 
condition  on  the  outside  of  the  ring  is 


(2.2-14) 


and  on  the  inside  of  the  ring  is 

The  stream  function  $  is  given  by  equation  (2.2-7)  and 
substituting  this  value  into  equation  (2.2-14),  the  boundary 
condition  on  the  outside  of  the  ring  is  obtained  in  terms  of  the 
perturbation  parameter  6  . 

[-(/?d  f  6[c‘‘\a  i-  (R4  r  e[c‘'\^)  f  s^''(fj]) 

( e[c‘'^a)  i-  r.i)\  6 


(2.2-16) 
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The  stream  function  is  a  function  of  +  s^^^z):] 
and  is  next  expanded  in  Taylor's  series  in  terms  of  this  para¬ 
meter.  Equation  (2-16)  then  becomes 


(S.2-17) 


The  stream  function  ‘'(R^+Q.f  ;  ^  )  can  also  be  expanded  in 
terms  of  6  ,  this  equation  then  becomes 

»  6  t  '  j  e[c‘'^'c?)  -h  s^'^'C^)] 

(2-18) 


Collecting  terms  of  the  same  order  the  nrevious  equation 
can  be  written  as 


£  /-  /?J  ivo  (Pj)[c^'^U)  t  -r  fjj 

' Lc^%:)  i- }  -h  ■  .  ,  =  n 
■  ^  C  J 

(2-19) 

The  first  approximation  is  obtained  by  neglecting  terms 
involving  powers  of  g  greater  than  one.  Derivation  of  the 


boundary  condition  for  the  inside  of  the  annular  airfoil  follows 
in  similar  manner.  Considering  both  sides  of  the  ring,  the 
first  approximation  for  the  kinematic  boundary  condition  is 

m  K, C nj) [c )  t  '(Vj  (2.2-20) 

The  +  sign  refers  to  the  outside  of  the  annular  airfoil  and  the 
-  sign  to  the  inside.  This  equation  has  only  to  be  satisfied 
on  the  circular  cylinder  o*'  diameter  (R^j)  and  length  (a). 

In  addition  to  the  foregoing  boundary  condition,  the  Kutta 
condition  must  be  satisfied  at  the  trailing  edge  of  the  ring. 
This  means  essentially  that  a  stagnation  point  must  occur  at 
the.  trailing  edge  of  the  ring  and  for  this  the  radial  velocity 
at  this  noint  must  be  zero.  Consequently  at  the  ring  trailing 
edge  the  stream  function  must  satisfy  the  following  equation, 

tO.At)  ::o  =  f/Vc/  to,  At)  (2.2-21) 

Since  it  has  been  assumed  as  a  first  anproximation  that 
1/^(1)  and  similarly  for  c^^^((f  ),  and  ),'  the  Super¬ 

scripts  in  the  'last  equations  can  he  drooped  and  the  bdiihdSrir 
conditions  can  then  be  written  as 

=  f?jw,CRj)Cc'(fJ  i:  SWJ  0^  a, 

^fCRdtOiAi)  =0  (2.2-22) 

For  convenience  the  radial  velocities  will  be  nondimen- 
sionalized  by  the  free-stream  velocity  if  this  velocity  is 
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uniform  or  by  the  ship  velocity  if  the  flow  is  axis3mimetric. 
The  axial  coordinate  will  also  be  nondimensionalized  by  the 
annular  airfoil  chord  a  and  the  radial  coordinate  by  the  pro¬ 
peller  radius  Rp,  i.e. 


z  =  axial  coordinate  = 

duct  chord  ^ 

X  =  radial  coordinate  =  — 
propeller  radius 


(2.2-23) 


Xd=  duct  rad:  .s 

propeller  radius 


if  the  duct  diameter 
is  equal  the  nropelle.* 
diameter  then  x^=l 


The  following  notation  will  also  be  introduced  which  is 


consistent  with  common  usage  in  naval  architecture. 


Cl  -  v/^}  =  (2.2-24) 

If  the  velocity  is  uniform,  i.e.  independent  of  radius, 
instead  of  the  shin  speed, the  free-stream  velocity  is  used  for 
nondimensionalizing  and  the  wake  fraction  is  unity,  i.e. 

Cl  ^  =  / 

In  the  definition  of  the  half-thickness  ordinate,  equation 
(2.2-2)  and  mean  line  ordinate,  equation  (2.2-1),  the  angle  of 
attack  of  the  section  was  not  discussed.  The  mean  line  ordinate 
is  measured  from  r=0  and  it  is  convenient  to  decompose  this 
ordinate  into  a  nart  from  camber  and  a  part  from  angle  of  attack, 
see  figure  3. 
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Figure  3.  Delineation  of  the  annular  airfoil  section 

Usx  .g  the  notation  given  in  Figure  the  mean  line  ordinate 
c(z)  can  be  expressed  in  the  follotdng  form. 


»» a 


I  —  »  V/ 


(2.2»25) 

It  may  often  be  assumed  for  practical  purposes  that  c^Cz) 
and  s(z)  may  he  measured  perpendicularly  to  the  nose-tail  line. 
This  is  for  convenience  in  delineating  the  camber  and  thickness, 
since  normally  the  nose-tail  line  is  used  as  a  reference  in 
describing  section  shapes,  and  implies  that  the  angle  oi  is 
small  '^n,„oc  ).  If  fhe  angle  oc  is  too  large  for  such  an 
««cui::pcior:  cnen,  although  inconvenient,  CfCz)  and  sCz)  must 
be  taken  perpendicularly  to  the  z-axis.  For  two-dimensional 
airfoils  it  would  be  expected  that  angles  for  which 
would  be  outside  the  applicability  of  the  linearized  theory 
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but  this  may  iioL  be  necessarily  so  for  the  sections  of  an 
annular  airfoil. 

If  the  foregoing  equation  for  the  camber  is  substititued 
into  the  boundary  condition,  equation  (2.2-22)  and  the  coordi¬ 
nates  and  velocities  are  nondimensionalized  as  discussed,  the 
boundary  condition  can  be  written  as 

to, 2)  -  -XJCi  -  tV^)[c/02)  V-  i  s'Cz)J  on  At 

d 

±0,at)  =0  (2.2-26) 

If  a  propeller  is  in  the  duct,  then  the  radial  velocity 
on  th«*  duct  is  a  function  of  angular  as  well  as  axial  position. 

In  this  case  the  singularities  used  to  represent  the  duct  must 
also  be  a  function  of  both  angular  and  axial  position.  Since 
the  shape  of  the  duct  is  assumed  to  be  axisymmetric,  (this 
assumption  can  be  removed)  the  right-hand  side  of  equation 
(2.2-26)  is  independent  of  the  angle  ^  ,  however,  the  left- 
hand  side  includes  all  the  radial  velocities  and  can  be  depend¬ 
ent  on  the  angular  position,  i.e,  the  following  equations  can 
be  valid  even  though  the  right-hand  side  is  independent  of  the 
angle  <f  . 

tC  ,  <P,Z)  =  rj  (I -Wrj)[c,'(2)  i-  t  s'(z)] 

(2.2-27) 

II. 3  Derivation  of  the  Vortex  and  Source  Distribution  to 
Represent  the  Annular  Airfoil  at  Zero  Incidence 

The  boundary  conditions  just  derived  will  be  used  to 
determine  the  strength  of  the  ring  vortex  and  source  distributions 
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on  the  cylinder.  These  ring  vortices,  of  elementary  strength 
and  ring  sources  of  elementary  strength  q(^,  ^), 
are  used  as  a  mathe.matical  model  to  represent  the  flow  around 
the  annular  airfoil.  In  order  that  the  distribution  of  sources 
and  sinks  represents  a  closed  body,  it  is  required  that  there 
be  no  outflow  of  fluid  from  the  source-sink  distribution  or 
that 

I  2ff- 

j  j  f(  t  d <pdz  -  O 

to  (2.3-1) 

If  a  propeller  is  in  the  duct,  then  the  ring  vortex 
strength  is  dependent  on  the  angular  nosition  and  a  trailing 
vortex  system  exists  behind  the  duct.  This  system  has  a 
strength  of  1  ) <K  and  the  induced  velocities  from  this 

system  must  be  added  to  that  of  the  ring  vortices  and  sources 
and  those  of  the  nropeller.  As  discussed  in  the  previous 
section  the  stream  function  occurring  in  the  boundary  condition 
is  the  total  stream  function  for  the  flow,  excluding  that  for 
the  free  stream,  and  since  it  is  linear  all  the  induced  veloci¬ 
ties  from  the  various  singularities  are  added  linearly  and  the 
boundary  condition  on  the  duct  is  expressed  by  the  following 
equation. 

^  fcf 

t  SU)]  (2.3-2) 


where 
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=  radial  velocity  induced  on  the  duct  by  the 
ring  vortex  system  and  given  by  equation  (A-16) 

=  radial  velocity  induced  on  the  duct  by  the  ring 
source  system  and  given  by  eouation  (B-8) 

=  radial  velocity  induced  on  the  duct  by  the 

by  equation 

=  radial  velocity  induced  on  the  duct  by  the 
propeller  including  the  hub. 

Making  the  substitutions  into  equation  (2.3-2)  for  the 

various  velocity  component.-  ,  an  equation  is  obtained  for  the 

vortex  and  source  -strength. 

>  sr 

^  I  a -r) cos  (<p-4>')  wa?'  ,  , 

J  j  -  B  to$((P~<f»)y^2.  ^ 

O  ^ 

3fh 

2  -  ’  cos(cp.r)r^^  -  2 


+  fcYzJ  ^.fa.cKtS  'Cz)J  -C\^^  (T(J, 

where  z  z.  z -a.^  (2.3-3) 

Since  the  integrals  occurring  in  this  equation  have  only 
one  sign  and  since  the  radial  velocity  induced  by  the  propeller 
on  the  duct  does  not  change  sign  from  one  side  of  the  duct  to 
the  other,  it  must  be  concluded  that  the  +  signs  go  with  the 

?  signs  and  hence  the  source-sink  strength  q(^,2)  is  given  as 
follows : 


which  imnlies  that 


(2.3-4) 
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fZ)  =  q(2) 

or  that  within  the  linearized  theory  the  source  distribution 
is  independent  of  angle.  If  equation  (2.3-4)  is  substituted 
into  equation  (2.3-5),  a  singular  integral  equation  is 
obtained  for  the  unknotm  circulation  distribution  which  also 
includes  the  derivative  of  the  circulation.  To  solve  this 
equation  for  the  circulation  distribution  requires  a  knowledge 
of  the  form  of  the  radial  elocity  induced  on  the  duct  by  the 
propeller  and  hub.  Before  deriving  this  velocity,  it  is 
convenient  to  consider  the  duct  by  itself,  first  .iL  zero 
angle  of  attack  and  then  at  an  angle  of  attack.  In  the  first 
case  the  circulation  distribution  is  independent  of  angle  and 
the  trailing  vortex  system  does  not  exist.  Utilizing  equations 
(A-23)  (B-11),  and  (2,3-4)  the  equation  for  the  circulation 
distribution  is  obtained  as: 

I 

^  J  (2.3-5) 

=  Vfr(/- Wxj)  -  ¥AO-^^r-J)  j 


where 


=  complete  elliptic  integral 
of  the  first  kind 

E(k)  =  de  -  complete  elliptic  integral 

o  of  the  second  kind 

h  =  a 
2Rd 
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From  this  equation  and  equation  (2.3-A)  it  can  be  seen 
that  the  vortex  strength  is  a  function  of  both  camber  and 
thickness  distribution  while  the  source  strength  is  a  function 
of  only  thickness.  This  differs  from  linearized  wing  theory 
where  the  vortex  strength  is  a  function  of  only  canber. 
Equation  (2.3-5)  is  a  singular  integral  equation  of  the  first 
kind  and  the  solution  of  this  equation  will  be  discussed  in 
the  next  secticn. 


II. 4  Reduction  of  the  Integral  Equation  for  the  Vortex 
Distribution. 
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Equation  (2.3-5)  Is  a  singular  Integral  equation  for  the 
vortex  distribution.  For  convenience  this  equation  can  be 
rewritten  in  the  following  form 

t 

^HCz)  (2.4-1) 

O 

where 


g  =  i  [ ^k^(2~zf[Ka)  -E(-^)]  -  2E{A))  (2.4-2) 


/ 


O 


(2.4-3) 


The  complete  elliptic  Integral  K(k)  has  a  logarithmic 
singularity  at  k  ■  1,  (z  ■  z').  This  causes  no  difficulty  in 
equation  (2.4-2)  since  llm(z-z*)^K(k)  ■  0,  but  the  Integral 


In  equation  (2.4-3)  Is  Improper.  A  logarithmic  singularity 
Is  removable  from  the  Integrand  of  an  Integral  by  a  change  In 
variable  and  this  technique  will  be  used  later. 

Muskhelishvili^®  has  shown  how  to  reduce  an  equation  of 
the  type  of  equation  (2.4-1)  to  a  Fredholm  equation  of  the 
second  kind  whose  solution  is  known.  Using  Muskhellshvili *s 
procedure  the  term  g(z’-z')  ■  g(0)  is  added  to  and  subtracted 
from  the  kernal  g(z-z')  in  equation  (2.4-1)  and  then  this  equa¬ 
tion  becomes 


nz'J 

(Z-z') 


iz' 


I 


(2.4-4) 
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From  equation  (2-, 4-2)  it  follows  that  g(0)  ■  -2  and 
equation  (2.4-4)  is  then 

I  / 

j(7^  -  1  (2.4-5) 

O  O 

The  integral  on  the  ri^thmd  side  is  not  singular  at 
the  point  z  »  z'  but  has  the  indeterminate  form  0/0  which 
can  be  shown  to  be  zero.  By  letting  the  right  hand  side 
of  this  equation  equal  f  (z)  this  equation  becomes 


(2.4-6) 


This  is  the  well-known  Cauchy  type  singular  integral 
equation^^  and  has  a  unique  Inverse  given  by 

The  last  term  is  the  total  circulation  about  a  section 
of  the  duct  and  is  a  constant.  In  this  equation  there  is  a 
singularity  at  the  trailing  edge  of  the  cylinder  (z  -  0, 
i.e.  z  -  aj.)  and  at  the  leading  edge  (z  ■  !)•  In  order  to 
satisfy  the  Kutta  condition,  equation  (2.2-21),  the  circula- 
tion  at  the  trailing  edge,  i^(0),  must  be  made  zero.^'  This 
is  accomplished  by  picking  the  total  circulation  so  that 
this  singularity  is  removed.  For  z  =  0,  then 

(2.4-8) 

0  O 

Substituting  this  result  into  equation  (2.4-7)  removes 
the  singularity  at  the  trailing  edge  and  the  circulation  dis¬ 
tribution  is  in  such  a  form  that  the  Kutta  condition  is 


satisfied. 
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(2.4-9) 


.  Substituting  in  for  f^Cz')  and  interchanging*  the 
order  of  integration,  a  Fredholm  equation  of  the  second  kind 
is  obtained  for  the  circulation  ^(z). 


where 


/=£:'  Cs.-h^Czl'-rJ.  rCzOdz' 


2tr^^(i-2)j  V  I"  (i"-2)(z"-Z') 


(2.4-10) 


a(z“-z')  =  ^  f¥hHz"-z'f[K(^)-E(.^)]  -  2  E(^)] 

(2.4-11) 

/ 

~  hUz"-Z‘f  4- ! 

This  equation  is  still  not  in  a  form  which  can  easily  be 
solved  since  Y (z)  has  a  singularity  at  z  ■  1,  i.e.  the 
leading  edge.  To  remove  this  singularity  a  new  dependent 
variable  is  defined  as  Y*iz)  =  V 1  -  z  Y(X)i  thus  at 
2  »  1  this  new  dependent  variable  is  zero.  Equation  (2.4-10) 
can  then  be  rewritten  as  follows 

I 

y^Cf)  --sji-z  f  (z)  =  {(z)  +  [k, (I,  z')  dz' 

J  ^  (2.4-t2) 


*  Interchange  of  the  order  of  Integration  of  two  Cauchy 
principal  value  Integrals  would  normally  result  in  a 
residue.  The  residue  is  zero  here  because  of  the  form 
of  the  integrand. 
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where 


Ol 

K(Z,r)  J;r« 


Both  f(z)  and  K^(z,  z*)  are  Cauchy  principle  value 
integrals  since  the  integrand  is  singular.  Since  it  is 
desired  that  a  solution  ir  thod  be  obtained  for  an  arbitrary 
H(z)  and  since  g(z"  -  z')  is  of  such  a  form  that  K(z,  z') 
cannot  be  obtained  by  simple  quadratures,  a  change  of 
variable  is  made  and  then  certain  functions  are  expanded  in 
a  Fourier  series.  Let  z  -  j(l  +  cos9j,  z'  ■  j(l  +  cos0') 

and  z"  ■  ^(1  +  cosQ"),  then  equation  (2.4-12)  becomes 


Y^e)  =  sin  f  s  y^e)  ::  f(e)  +  \t<.(eje')  Y%')de‘ 


(2.4-13) 


where 


He)  r  -  cos^olj'  ■  H(e')cle' 

^  J  (cos  o'  -cose) 


(2.4-14) 


IT 

it  ^  2  j(coj$"^cose)[  (ci’ie"  -cose)  J  (2.4-15) 

6 

To  find  the  Cauchy  principal  value  of  the  integral  for 
f(0),  the  function  H(0!)  will  be  obtained  in  a  different  form. 
From  equation  ^2.4-14)  H(9’)  is  obtained  as 


IT 

■n(&)  =  ‘f(i- js'(6"J^[m  -E(f\)]sm&"Jej  (2.4-16) 
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where 


|2 _ V _ 

”  k^(ccse'-cosel'}‘'  +  i 


The  thickness  distribution  and  the  mean  line  shape  are 
now  expanded  in  a  Fourier  sine  series  in  9 ' .  Because  of  the 
shape  of  the  airfoil  section  it  would  be  expected  that  such 
a  series  would  converge  very  rapidly.  The  use  of  such  a 
Fourier  series  implies  tha  the  required  slopes  can  be 
obtained  by  a  term  by  term  differentiation  of  the  series. 

For  this  to  be  possible  the  mean  line  and  thickness  distribu¬ 
tions  must  satisfy  additional  restrictions  than  would  normally 
be  necessary  for  their  expansion  in  a  Fourier  series.  Specif¬ 
ically  for  the  Fourier  series  of  a  function  to  be  differentiated 

no 

tenn  by  term  the  function  must  be  every^^here  continuous  and 
possess  a  derivation  which  satisfies  the  Dirlchlet  conditions. 
For  practical  sections  this  causes  no  restriction  on  the 
shape,  even  if,  the  slopes  are  infinite  at  the  ends  since  the 
requirement  in  this  case  is  that  the  Integral  of  the  slopes 
be  absolutely  convergent.  Expanding  the  thickness  and  camber 
distribution  in  a  Fourier  sine  series  in  0,  the  following  are 


obtained. 


c,(e) 

Cm 


(2.4-17) 


(2.4-18) 
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The  slopes  are  obtained  by  differentiation  and  are 

c:m  Cn%  m  e 

mzf 

OO 

IH  Hi  C0S  ^  0 

(2.4-19) 

Introducing  these  expressions  for  the  slopes  into 
equation  (2.4-16)  and  Interchanging^S  the  order  of  integration 
and  summation,  H(e)  is  ob  dned  in  the  following  form. 


IK 


H(e‘j  - 


)  ■=(!  -  >^rj)  [tan<x.  j-  ]^c»,>nCi>s  m&J  - p ~  £'^^)Jci>si’te"sin 

ptiSf  J 

O 

(2.4-20) 

The  elliptic  integral  of  the  first  kind  K(k)  has  a  log¬ 
arithmic  singularity  at  k  -  1,  i.e.  when  cos  0"  •  cos  0'. 

This  singularity  can  be  removed  from  the  integral  by  considering 
the  change  in  variable  (cos  0*  -  cos  0")  -  t^,  then 


it' 

r 


'kltcois') 

j-^CK(l)-E(^)]s,»Md"stne"a'0"=  j - E(£)]coif» [rK^ccoi[co}o'-t^]J 

-coie') 


dt 


-  ZG(Q',m) 


(2.4-21) 


where 


This  removes  the  singularity  and  the  integral  can  He 
evaluated  numerically  without  difficulty.  To  complete  the 


28 


solution,  the  functiou  G(w*,in)  is  expanded  in  a  Fourier 

cosine  series  in  0',  i.e- 
i-ccsa') 

I  >>0 

E(J()'jcoSt*i  f  CLi'CC0s[coi6)' cos  4:^0' 

I  -  j  i^zO  ' 


-  -\j(i-cose') 


(2.4-22) 


T'rtiere 


7f 

b 

'TT 

dfo(M)  =i.  f6(0jntjs//  e'dd'  J  ( -^  =  /jHy  J  ■  '  <^) 

^6  (2.4-23) 

Using  this  representation  of  G(0',m)  and  substituting 
equation  (2.4-16)  into  (2.4-14),  the  integral  for  f(0)  can  be 


evaluated 


{(6)  -  -|(/- 


-  ^Kj) — Ae'  +cos\eZ^»<'^'\^^^^^^-cosMacJ0 
I  2  \(.Cgl6' -cose)  ^  hiS!  j(COS0'-CoS0) 


(2.4-24) 


The  integrals  occurring  in  this  equation  are  of  the 

Glauert  type  and  their  evaluation  is  given  in  Reference  [27] 

fr  ^  ^ 

[(tzCfM^cosme'de'  de'  -  \cose!mm£'  de' 

\(cose'-cose)  j(cos$'-co3e  jdose' -cose) 


-  _  McoSh^&'+  S.2.^ecosn,a‘ 


Ji  {.Cose' -Cose). 


jc/e' 


zz  irCi  -cose) 

'  s/»e 


(2.4-25) 
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Also 


IT 


I  (t  ~ cose^j^t  _ 

J(cas&'-co::o} 


Using  these  values  for  the  integrals,  f(9)  follows  as 


fCa)  =Z(l  -  jtgnoLCOS-^$  - 

L  jMi/  ms! 

~  ( I  ~  ia  n  <X-  'f' ^■^'»^jCc>S'^&  ■f~ 

MU  ^  /M=/  ”'  ' -^ 


(2.4-26) 


where 


(2.4-27) 


The  coefficients  Fjjj  and  Bjj(O)  are  independent  of  the 
section  shape,  however,  they  are  dependent  on  the  chord- 
diameter  ratio  of  the  duct,  and  can  be  tabulated.  Now 
consider  the  equation  for  the  kemal  K(0,9*),  equation 
(2.4-15).  Attempts  at  evaluating  this  improper  integral  by 
simple  quadratures  have  been  unsuccessful  so  K(9,9')  will 
be  obtained  in  a  form  which  can  be  solved  by  numerical  methods. 
In  the  following  method  the  term 


s  +  ^  (cose" -cose' 
coie"  -  cose' 
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Is  expanded  in  a  Fourier  cosine  series  in  0"  in  the  range 
0  ^  0"^  .  This  function  satisfies  the  Dirichlet 

conditions'^  so  can  be  expanded  in  such  a  series.  Further¬ 
more,  the  function  is  continuous  for  0-  but 

the  first  derivative  is  discontinuous  at  6"  ■  9*.  Expanding 
this  term  in  an  even  scries,  it  becomes 


o<> 

(case" -C0S9') 


where 


0  Vj  COSO  -COiO' 

0  ^ 

D- 

"  fTj  [  cose  -COSO' 


(2.4-28) 


(2.4-29) 


The  Integral  for  the  coefficients  must  be  evaluated 
numerically  but  some  difficulty  arises  at  the  point  9"  -  9' 
vdiere  the  integrand  has  an  Indeterminate  form.  To  determine 
the  value  at  this  point  1' Hospital's  rule  is  applied  to  the 
Integrand. 


2  ^  ^ (cos a"-' Case) 
Cose" -cose' 


COSfl  0 


„  ^  _  c),„(0koshe'  _  Q 


(2.4-30) 


Since 


-s>ne"(cos0"-coso')^l  (^£1^!^ 
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At  the  point  6”  -  S*,  therefore,  the  integrand  is  zero, 
however,  because  of  the  form  of  equation  (2.4-30),  it  must 
be  shown  that  when  Q*  =  0  or  'ft'  equation  (2.4-30)  is  still 
valid.  Following  the  same  procedure  as  above  for  9’  -  0,  then 


ti^*\ 

2  ^■t](co5e‘'~l) 

cosy^e"  — 

cjs"(<^o^9"-i)cos e"  -ns/nn0'' [2  t  <i(ce>s 

cose"-!  . 

&-*0 

—  smo" 

llhi 

dUo 


(cose"-l)  - £(^}\  *  I  ^  7gjf"  2  +^(coso"-l) 


_  IlM 

~  e"-*o 


(ns'ihhe")^Q>i  (cosd"-i)  -^-h^cosM  e"[2i-j(co$e''~i)j 
cos  e" 


O 


From  this  last  result  and  equation  (2.4-30)  it  can  be 
concluded  that  the  Integrand  of  equation  (2.4-28)  has  the 
value  of  zero  for  9'*  «  0*,  0^  0’-/)'  .  Substituting 

equation  2.4-28)  into  the  kernal,  equation  (2.4-15),  the 
integration  can  be  performed. 


it 


ft 


=  ^iCos±e'cos 


L^r.C£^  Cos^&"Jd" 
J  (cose  ^^cose) 


I  /•  ,€j^  o' 


•  ^ 
-(cos^&)bo(e')  +( (9  ')sih  n  e 

t^-.\ 


(2.4-31) 
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II. 5  Ideal  Angle  of  Attack 

In  the  Integral  equation  (2.4-10)  for  the  circulation 
distribution  a  singularity  exists  at  the  leading  edge  of  the 
ring  vdiich  made  it  necessary  to  redefine  the  circulation 
distribution.  This  singularity  cannot  exist  if  the  ring 
section  is  designed  so  that  a  stagnation  point  occurs  at  the 
leading  edge.  The  angle  of  attack  at  wii:h  the  section  is 
operating  when  a  stagnation  point  occurs  at  the  leading  edge 
is  known  as  the  ideal  an/ 'e  of  attack  Since  an  axi- 

symmetrlc  annular  airfoil  has  been  assumed,  it  is  obvious 
that  an  ideal  angle  of  attack  cannot  be  defined  if  the  radial 
velocities  on  the  duct  surface  are  functions  of  angle. 
Consequently,  in  the  presence  of  a  propeller  an  ideal  angle 
of  attack  cannot  be  defined  which  applies  to  every  section. 

To  remove  the  singularity  occurring  at  the  leading  edge, 
l.e.,  2  *  1,  in  equation  (2.4-9),  this  equation  is  rewritten 
as  follows: 


I 


where 

f.  w)  =  i-  Z(i~t^xj)rrtx^o(. 

Following  the  arguments  for  satisfying  the  K”tta  condition, 
a  stagnation  point  will  occur  at  the  leading  edge  of  the  duct 
if  the  circulation  at  that  point  is  zero.  The  circulation  will 
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be  zero  at  the  leading  edge  if  the  angle  o<  is  chosen  so 
chat  the  terms  In  the  brackets  of  the  previous  equation 


cancel  out  the  singularity  of  the  term 


at  z  «  1. 


The  ideal  angle  of  attack  is  therefore »  defined  as  follows 


tan  ^  -  -- 


^rj} J 


if/'-/  V  2' 


f,U') 
M  -z‘)Z' 


(2.5-2) 


The  term  f^^Cz’)  occurring  in  this  equation  is  a  function 
of  the  ciifculation  distribution  Y  (z)  and  this  distribution 
must  correspond  to  that  occurring  at  the  ideal  angle  of  attack. 
To  determine  the  so  called  ideal  lift  coefficient,  ^ (2) 
equation  (2.5-2)  is  substituted  into  equation  (2.5-1).  The 
solution  of  this  Integral  equation  will  give  the  circulation 
distribution  for  the  section  operating  at  its  ideal  angle 
of  attack. 


f.(z')dz' 


(z'-  Zh^jz'c  I -z.y 


(2.5-3) 


From  the  form  of  this  eqmition  it  can  be  seen  that  the 
singularity  is  now  removed  from  the  leading  edge.  Proceeding 
as  in  the  last  section  it  can  be  shown  that  this  equation 
can  be  reduced  to  a  Fredholm  equation  of  the  second  kind  for 
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the  circulation  distribution.  The  following  form  is 
obtained  for  this  equation  where  the  coefficients  are  the 
sarae  as  given  in  the  previous  section. 


r  ^  ^ 

L  /«=/  «•=/ 

ir 

+  ^\  smO'Y^ni'^'^sinhe 

Jl 

^  ( 
By  substituting  into  equation  (2.5-2)  it  follows 
that  the  ideal  angle  of  attack  can  be  written  as: 


(2.5-3) 


fo-n  = 


#•  cc 


(2.5-4) 
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II. 6  Solution  of  the  Fredholm  Equation  of  the  Second 
Kind  for  the  Circulation  Distribution. 

Both  equations  (2.4-3!^)  and  (2.5-3)  can  be  solved  by  the 
same  procedure.  The  general  form  for  either  equation  can  be 
written  as  follows 


ne)  =  0-^tj)jBo(e)  +B(e) 


op 

Bm  (e) 


L— 

»/=/ 


ftt-l 


IT 

r  oo 

■t-c(6)  UDOO  Do(e)  1,^(6')  ■hYj>n(&'Jsi>lhe 


}  y(&')de' 


(2.6-1) 


In  the  general  case: 


r(9)  -  r*(0) 

CO 

Bo(6)  “  [2  tan°c  +  ^S^  Fj^]  cos  y  0 

3(9)  -  sin  ^ 

C(0)  ■  —  sin  T  0 

ir  ^ 

D(0^  ••  cos  ^  0^ 

Do(0)  “  ■  cot  j  0 

and  for  the  Ideal  circulation  distribution 

r(0)  -  (0) 

Bo(0)  -  0 

3(0)  -  1 

C(0)  =  ^ 

D(0')  -  sin  0' 

Do(0)  -  0 
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The  Fourier  coefficients  are  the  same  as  defined  in 
Section  II. 4.  The  integral  equation  (2.6-1)  can  be  solved 
by  the  method  of  successive  approximation  or  the  Fredholm 
solution  method. In  addition,  since  the  kernel  in  this 
equation  in  a  degenerate  (or  product)  kernel,  the  special 
method  appropriate  to  this  type  of  kernel  can  be  used. 

Which  method  is  best  depends  on  the  convergence  of  the 
Fourier  series  representl 'g  the  kernal.  If  the  series 
converges  rapidly  the  method  appropriate  to  product  kernal 
is  probably  best.  Folloid.ng  is  an  outline  of  this  method, 
necessary  proof  of  the  convergence  is  given  in  Reference  [29]. 

To  apply  this  method  the  order  of  integration  and  summa¬ 
tion  are  interchanged  in  equation  (2.6-1)  and  then  the 
equation  can  be  vnritten  as 


■ir  ^ 

Y(6)  -((&)+  C(d)  ^<)(&) DO')  i-  sih& j  (t>)  D(&')i  '(b 


<(■ 


nr 

r 


+  Sin 


ze^b^Ce'} 0(e') -t- •  *  ♦  5'mn0jbrtO'jD(e')^(e')cie' 

"  (2.6-2) 


\diere 


r  r  ^  oo  ■ 

f(e)  =0 - +  B(e)\-z'^t»ryisih me  B^^Ce) 

^  ^  ms./  Ws/ 


For  convenience  let 
nr 


A  =  J b„(e') D(0')  e'  ^  (c  =  Oj  Ij  2,  •  •  -  n) 

o 

and  then  the  above  equation  can  be  rewritten  as: 


)r(e)  =f(&)  +  c(d)  DJe)  Ao  ■hC(0)A,sin0+  •  •  •  -f 


(2.6-3) 
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By  substituting  this  equation  into  equation  (2.6-2)  the 

A  's  are  obtained. 

J 

it' 

/lo-  Jors‘jb„(a'j[fM  +c(a') pje')Ao-f-c(6»)A,s)ne'-h 


a(6')A,sih&‘-i- 


=.J  D(e')bn(e‘)[f(&')-hC(e')Oj60Ao  +  c(&')A.S'bd' - 

^  (2.6-4) 

In  deriving  the  equation  for  Aq,  it  was  necessary  to  divide 

both  sides  of  the  equation  by  Do(0)  and  if  Dq(0)  is  zero  this 

is  not  valid.  For  this  case  then,  i.e.  ideal  circulation 

distribution,  A^  must  be  taken  as  zero  and  the  above  set  of 

(n  +  1)  equations  reduces  to  a  set  of  (n)  equations.  For 

convenience  the  following  notation  is  Introduced  into  the 

preceding  equation. 
ir 

c<y  =  ^O(e')bi(e')c(0ODj(e')de'  ; 

o  ir 

=  ^jbi  (a')s/b e'de'  (2 . 6-5) 

O 

where 

OM^-cot^e' 

Df(e')  =  ,  (^  -  l,2.j  '  ■  ‘  h) 

ir 

=J  D(e‘)  biCeOHeOde* 


and 
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This  system  of  equations  represents  an  algebraic  set 
of  simultaneous  equations  for  the  unknown  A^’s.  The 
existence  of  a  unique  solution  depends  on  the  determlnent 
of  the  coefficients  of  on  the  left  hand  side  being 
different  from  zero  (Cramer’s  rule).^®  The  number  of 
simultaneous  equations  depends  on  the  number  of  terms 
needed  in  the  Fourier  series  so  it  satisfactorily 
approximates  the  kernel,  fortunately  this  can  be  determined 
once  and  for  all  since  the  kernel,  K(0,  O'),  Is  Independent 
of  section  shape. 
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I I . 7  Pressure  Distribution  and  Forces  on  the  Annular  Airfoil. 


The  velocity  field  of- the  ducted  propeller  Is  found  by 
summing  the  free-stream  velocity,  the  velocity  Induced  by 
the  duct  Including  the  trailing  vortex  system,  and  the  velocity 
Induced  by  the  propeller  and  hub.  Since  the  flow  field  has 
been  asstmed  to  be  Irrotatlonal,  steady  and  incompressible 
and  the  body  forces  have  been  neglected*  Bernoulli's 
equation  can  be  written  as  follows: 


=  -  ([-t  - 

[%0' 


■t 

J, 

**  i 

(2.7-1) 


The  axial  velocity  w  (x,  (f ,  z),  radial  velocity  Vr(x,<p,z) 

A 

and  tangential  velocity  w^(x,  (p  ,  z)  are  the  total  velocities 
Induced  by  the  various  singularities  In  the  flow  and  are 
commonly  called  perturbation  velocities.  The  pressure  p^ 

Is  the  pressure  infinitely  far  ahead  of  the  propeller 
while  p(x,<f^,  z)  Is  the  local  pressure.  Tills  pressure 
distribution  has  been  nondimenslonallzed  by  the  ship 
velocity  times  the  wake  (1  -  Wx)  which  is  the  local 
free-stream  velocity. 

If  this  equation  is  linearized  by  the  same  method  as 
used  for  the  linearized  boundary  conditions,  the  squared 
terms  will  be  neglected  and  the  linearized  prressure  distri¬ 
bution  Is  then 


-fo  _  2  W^(Xj<P.z) 


(2.7-2) 
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The  pressure  on  the  annular  airfoil  Itself  is  found 
by  substituting  the  duct  radius  for  x  In  the  preceding 
equation . 

-  g  <P,z) 

(2.7-3) 

The  velocity  w-(x»if,  z)  is  the  axial  velocity  induced 
by  all  the  singularities  'n  the  flow.  Also,  It  should  be 
noted  that  since  the  velocity  induced  by  the  vortex  cylinder 
Is  discontinuous  across  the  cylinder,  equation  (A-15), 
the  pressure  changes  from  the  inside  to  the  outside  of  the 
duct.. 

In  the  problem  of  a  ducted  propeller  in  an  luviscld 
fluid,  the  only  net  force  on  the  duct  Itself  Is  the  so-called 
Induced  drag  or  force  in  the  axial  dlrextlon.  There  Is  a 
radial  force  on  each  section  which  contributes  to  hoop  stress 
but  because  of  symmetry  of  the  flow  this  net  force  Is 
zero.  Since  the  net  lateral  force  Is  zero,  there  Is  no  moment 
on  the  duct,  however  a  moment  on  epoh  Individual  section 
could  be  defined. 

The  force  F  on  any  section  of  the  duct  is  given  by 
the  Kutta-Joukowskl  law^^  which  can  be  expressed  as 

F-/oVr  (2.7-4) 

The  velocity  Y  Is  the  velocity  by  the  annular  airfoil 
section  perpendicular  to  the  direction  of  the  force  and  P 
is  the  total  circulation  about  each  section.  The  velocity Y 
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does  not  Include  the  self  Induced  velocities  so  does  not 
contain  the  velocities  induced  by  the  vortex  and  source 
rings.  Assuming  that  both  the  veludcy  and  circulation 

z)  are  nondimensionalized  by  the  ship  speed  and  the 
axial  coordinate  z  by  the  chord,  the  lift  at  each  section  is 


dL 


(2.7-5) 


The  Induced  velocity  v  (xj,  <p  , 


is  the  axial 


velocity  induced  on  the  cylinder  of  radius  x^  by  the  propeller 
hub  and  duct  trailing  vortex  system.  If  the  total  lift  on 
the  vdiole  ring  is  taken  to  be  some  arbitrary  direction  with 
respect  to  a  propeller  blade,  normally  in  the  vertical  direc¬ 
tion  and  positive  upward,  the  contribution  to  the  lift  of 
any  arbitrary  section  of  the  ring  would  be 

I 


<iLcesi>  =  (cos 


(2.7-6) 


> 


The  angle  is  measured  from  the  arbitrarily  taken 

direction  which  is  attached  to  one  of  the  propeller  blades, 

i.e.  rotating  coordinate  system.  If  this  equation  is  integrated 

completely  around  the  circumference  of  the  ring,  the  net  lift 

will  be  zero  since  the  flow  is  axl83mmetric.  If  the  annular 

airfoil  is  at  an  angle  of  attack  to  the  flow,  the  flow  is 

% 

no  longer  axlsymmetrlc  and  there  is  a  net  lift  force. 

The  thrust  force  on  the  duct,  which  is  equal  and 
opposite  to  the  induced  drag,  is  also  given  by  the  Kutta- 
Joukowski  law,  equation  (2.7-4),  but  in  this  case  the 
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velocity  at  the  duct  is  the  radial  velocity,  again  not 
Including  the  self  5.nduced  velocities.  Denoting  the  duct 
thrust  by  j  the  following  equation  results  fcr  the  thrust 
force  on  each  section. 

o  ^  ^  (2.7r7) 

and  the  total  thrust  is  given  by  Integrating  this  equation 
around  the  ring. 

I  tir 

(2.7-8) 

This  equation  car.  be  put  In  the  form  of  the  thrust 
coefficient  used  In  propeller  design  by  dividing  by  ^fff,irVs\ 


- 


(2.7-9) 


The  subscript  "1"  in  this  equation  means  Inviscld  fluid. 
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II . 8  Linearized  Theory  of  the  Annular  Airfoil  at  an 
Angle  of  Incidence. 

In  addition  to  the  assumptions  made  previously  concerning 
the  annular  airfoil,  It  will  now  be  assumed  that  the  free- 
stream  velocity  Is  a  constant  but  at  an  angle  ^  to  the  ring. 
The  angle  will  be  assumed  to  be  small  enough  so  that  8ln<^j. 

“  tan*^  j.  j.  “  1.  As  previously,  a  cylindrical 

coordinate  system  (x^ep ,  £)  will  be  used  with  the  zero  axial 
coordinate  (z)  zero  at  the  ‘.railing  edge  and  (R^j)  the  reference 
diameter  will  be  located  at  the  propeller  centerline. 


Figure  5.  Free  stream  velocity  at  an  angle  to 
the  duct 

Applying  the  perturbation  theory  as  in  Section  11.2,  and 
assuming  the  angle  of  incidence  o< ^  to  be  small,  It  can  be 
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shown  that  the  linearized  kinematic  boundary  condition  Is 

^(Rjto^ip^z)  +  Wo  <^^C0S(p  =  -  Wo[c'(^)  t  s'(^)] /os-z^  I)  (2.8-1) 

In  addition  to  this  boundary  condition,  the  Kutta 
condition  must  be  satisfied  at  the  trailing  edge  of  the  ring. 

fz(R^tO,cf,o)  -^0  (2.8-2) 

One  difference  betwe  t  equation  (2.8-1)  and  the  boundary 
condition  for  the  annular  airfoil  at  zero  Incidence  Is  the 
addition  of  a  radial  velocity  term  from  the  free-stream 
velocity  which  Is  dependent  on  the  angular  position.  This 
then  implies  that  the  radial  velocity  Induced  by  the  annular 
airfoil  must  also  be  a  fxmctlon  of  the  angular  position  of 
with  reference  to  the  mathematical  model  the  ring  vortex 
strength  Is  a  function  of  angle.  As  for  a  three-dimensional 
wing^^  with  a  spanwlse  change  in  vortex  strength,  the  fact 
that  the  vortex  strength  Is  a  function  of  angular  position 
leads  to  a  trailing  vortex  sheet.  This  vortex  sheet  is 
assumed  to  be  cylindrical  In  shape  and  to  extend  to  Infinity 
behind  the  annular  ring. 

The  radial  velocity  w^  In  equation  (2.8-1)  Is  mathematically 
conceived  as  being  a  sum  of  the  radial  velocity  induced  on  the 
ring  by  the  r^ng  vortices,  sing  sources  and  trailing  vortex 
system.  The  Integral  equation  for  the  ring  vortex  and  source 
strength  vaa  derived  In  Section  II. 3,  equation  (2.3-3),  and 
further  the  soiurce  strength  was  shown  to  be  Independent  of 
the  angular  position,  eq\!atlon  (2.3-4).  The  velocity  Induced 
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by  the  trailing  vortex  system  from  the  duct  at  an  angle  of 
attack  follows  from  the  law  of  Biot-Savart,  and  is  derived 
in  Appendix  C,  specifically  equation  (C-12).  If  this  equation 
is  substituted  into  equation  (2.3-3)  along  with  equation 
(2.3-A)  an  integral  equation  for  the  vortex  circulation  is 
obtained.  In  this  section  the  flow  is  aisimed  uniform  so 
(  1 “  ^xd^  “  1 » 


h(z-z'J 


sir 


/  cos(f-  ')/Cip,z')  d(p' 

[  J  -zf-h  Hsin*i(tp-(p‘)y^^ 


dz' 


I  sir 


L  2  T/J  r.  ,2  ^ 


o  0 


-  +  / 


rr[am  H  <=>c^cos(fj  +  h  fs'(7f)J^[KrJl) 


=  -  H(E)  -  i^TToCr  C0S(p  =  Hz  (Z,(p) 


(2.8-3) 


where 

The  term  is  the  same  as  given  by  equation  (2.4-3) 
and  the  evaluation  of  the  integral  in  H(z),  which  has  a 
singularity  in  the  integrand,  is  discussed  in  Section  II. 4. 

Equation  (2.8-4)  is  a  singular  integral  equation  of  two 
variables  for  the  circulation  distribution.  The  inversion 
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of  an  Integral  equation  of  this  type  is  not  in  general 

10 

possible.  Weissinger  encountered  this  Integral  in  the 
problem  of  the  thin  annular  airfoil  at  incidence  and  x<!duced 
the  two  variable  integral  for  the  circulation  to  a  one- 
dimensional  Integral  in  terms  of  the  Fourier  coefficients  for 
the  circulation.  Following  his  procedure  the  circulation 

and  the  function  H2(s, <^)  are  expanded  in  a  Fourier 
cosine  series  in  cf  .  This  Involves  a  restriction  on  the 
circulation  that  it  must  ’  continuous  in  the  angular 
direction,  but  from  a  practical  point  of  view  this  presents 
no  difficulty. 


where 


oO 

mo 

OOi 

Ur>(E)cosfi<P 

«s<? 

ir 

Jn 

'tr 

=  ^J  rlf,r)cosh(py<P' 


(2.8-4) 

(2.8-5) 

(2.8-6) 


ir 

-r 

^  (2.8- 

u^Cz)  =;l^jHtC(P,E)ccrsncpJ(f 

-TT 

A  cosine  series  is  used  since  the  flow  is  symmetric 
about  a  vertical  plans  in  the  direction  of  the  z-axls  (See 
Figure  5)  and  hence  the  circulation  and  radial  velocities 
are  even  fimctions  of  , 
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Substituting  equations  (2.8-4)  and  (2.8-5)  into  equation 
(2.8-3)  and  Interchanging  the  order  of  suomation  and  Integra** 
tion,  the  following  is  obtained 


t 


i 


U)2k(z~z‘) 


ir 


Jz' 


-r 


S  in  n  'O' cot i  (<P-<P^)  d  cP‘  1^1 

CvaYz 


/  fr 


ir 

'  j  sinh<P'coti((p-r)cl<p''^dz] 

-ir  ^ 


(2.8-8) 


but,  see  Reference  ^271 , 

0  (2.8-9) 

and  making  a  change  of  variable,  0'  »  (^^^)  and  using  the 
trigonometric  Identity  cos  n<f*"  cos  2n0^costa^-  sl.n  2n0*  sin  n^ 


TT 


I 


tr 

”E 


/  cos(<p  ~<p')cosn(py<p'  _ ,  /  coiB9*co5n<p'de' 

J +  "  J [‘fhHz-zT i- 

-2crinr[  coiSf^6'cosB$'de'  slmhe'coszs'de 


y  Qg  Pc/Mnifi  ^fncn^  ucf 

^j[^h^C2  -zf  f  V sm^0'p't- 


% 

=  cossne'cossfi'dtp' 


(2.8-10) 
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The  integral  Involving  sin  2n9'  is  zero  because  it  is  an  odd 
function  of  9*.  Again  introducing  the  new  variable  9'  - 
(  ^  )  and  the  trigonometric  identity  sin  n  *  • 

sin  2n9'  cos  n  <p  +  cos  2n0'  sin  n  (f 

r  ^ 

I  s/tin>p  cot^(cp~(p')j(p'  ^  ^  t^r.r>inepf sihane'cofe'tJe 

-IT  n 

(2.8-11) 

Substituting  equation  (2.8-9),  (2.8-10)  and  (2.8-11) 
into  equation  (2.8-8),  the  follov?ing  equation  is  obtained. 


t 

Y^un  (2)cosnep 

hsO  naa  \J 

o 


^‘i(eosah  $'coszb'+  nzimne'sim0*)eld* 
[‘fh*(z~2')  -h  <fsin^&‘J^ 


^<r^LS/4  I  iifiene'cote'do' 

a-jrl 

»  i  I 

(2.8-12) 

Since  a  Fourier  series  is  unique,  the  coefficients  can  be 
equated  and  an  Integral  aquation  of  one  variable  is  obtained 
for  the  Fourier  coefficients. 

'  J 

Un(2)  ~2hl <j^(r}2h(E  -zo\  I 


‘f(cosahe'cosee'  -t-  ns'/nzh^^smss')  j^, 


^  h^Cz-Z'f  _ 

0 


i-Ztfhh  Un(Z')dz' 


(2.8-13) 

Before  proceeding  further  with  the  solution  of  eq«aation 
(2.8-  3)  the  Fourier  coefficients  will  be  examined. 
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Substituting  for  z)  into  (2.8-7)  the  Uj^(z)'8  are  given 

as 


u  (z)  -  -  H(z) 

(2.8-14) 

0 

Uj^(z)  «•  -47ro^y 

(2.8-15) 

Uj,(z)  -  0  (n  -  2,3...co) 

From  these  equations  it  can  be  concluded  that  all 
for  n  2,  3  ...  oe  are  zero  and  therefore  only  g^  and  gj^ 
exist,  also  that  gg  is  a  f  notion  of  shape  only  and  gj^  of 
angle  of  attack.  First  examine  the  circulation  coefficient 
g(x  ,  letting  n  «  0  in  equation  (2.8-13). 


This  equation  is  of  course  exactly  the  same  as  (2.4-1), 
thus  it  can  be  concluded  that  g^(z)  »  (z)  where  K(£)  is 

the  circulation  distribution  for  the  annular  airfoil  at  zero 
incidence.  This  eqtaation  (2.8-16)  is  solved  in  Sections  II. 4 
and  II. 6.  It  should  be  mentioned  that  on  applying  the  Kutta 
condition,  equation  (2.8-2),  the  equation  for  the  Fourier 
coefficients  g^,  equation  (2.8-6),  implies  go(0)  «  8^(0)  ■  0. 

As  has  been  stated,  gj^(z)  is  independent  of  the  shape  of 
the  airfoil  and  dependent  only  on  the  angle  of  attack  and 
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chord-dianeter  ratio  h  .  For  n  »  1  (2.8-13)  becomes 

/  _  Ji 

,  .yif /„!/./-  ,/)  /  sin^£0')de' 

=Y/>  g,(zJh{z-zJ 


- g£P>sV  de'  +srrh  g,U0dz' 


-iv<r  »  Z 


_  t^o-4*)^K(^)dz'  +  Zfrhj^,(2')dz^ 

o  O 

/  I 

=  2  V,  C2  -ridr  +zvrh^g,  (Z')dz' 


where 


(2.8-17) 


-  fd-  (2.8-18) 

J?*-  ' 

^  '  h^(2-Zf  ^  / 

It  can  be  easily  shown  that 

.'ij  =  I 


By  adding  and  subtracting  - - 7“  equation  (2.8-17)  can 

(z  -  z') 

be  obtained  In  a  form,  similar  to  equation  (2.4-5). 

/  I 

J='  I  u-z)  J 


where 


//>H  f  W  CZ  -Z')  -  I  1  =:<9 
z'-^/L  rf-7'J  J 
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and  this  eqtiatlon  can  be  solved  by  the  same  procedure  as  given 
in  Sections  II. 4  and  II. 6.  Since  the  solution  follows  exactly 
that  of  the  previous  sections  with  only  changes  h\  the  Kernel 
K(z,  z*)  and  in  f(z)»  a  complete  derivation  will  not  be  given 
here.  Following  the  procedure  of  Section  II. 4  yields  for 
equation  (2.8-20)  the  following  Fredholm  equation  of  the 
second  kind^  where  z  -(  -j)(l  +  cos  9). 

0  ' 


where 


■h(Sln  j3i^bm(9')sintne 

msi  •’/ 


(2.8-20) 


ir 

o 

•r 

L  (O')  -  2  /[  V>  ■'<^0^')  - 

Uaso"  ^cose')  J 

O 

(2.8-22) 

From  Section  II. 6  the  solution  of  (2.8-20)  is  obtained 
in  the  following  form. 


^*(ei  =  B°^r  CoS^Q  -  ^s'i>^-^9^AgCot-^9  t  A,s'in9 


+  •  *  '  *  1 


(2.8-23) 


The  ^'s  are  obtained  from  a  solution  of  the  following  set 
of  algebraic  equations 
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Ao(l-Coo)-  ...... 

“Vio  +  A^(l  -  -  A2C^2  •  •  •  -  -  d^ 


"^o^mo  "  "  ''^2Cm2  - 


+Am(l-C^  -  «*„ 

(2.8-24) 


The  coefficients  C  are 
nn 


C^m  =  Cii  =  - 4- j '3i(&‘) D^(6‘)sm 0'd» ',  (c,f  -o,ht  •  •• 


or  more  specifically 


ir 

-  +  -h ^jl^\e')(i  -t-coseOde' 

ir 

Co,  =  -fA  ^(0')sm^9de' 

ir  ^ 

o 

ir 

c<o=  :lj,jk(e')(/i-coi0')de'  -j 

If 

Ci^-  jbi(e')s'/hfe's’m0'e/&'j 


if  *£.  '  •"«/ 


(2.8-25) 


(f  =  ' '  '/»»y 


(i,^  = 


In  the  preceding  equations  for  C^j  the  coefficients 
b(l>(0)  and  bjL(0)  are  given  by  equations  (2.8-21)  and 
(2.8-22)  respectively.  The  coefficients  dj,j  occurring  in 
the  set  of  equations  (2.8-24)  are  dependent  on  the  angle  of 
Incidence  of  the  annular  alrfoii.  They  can  be  shown  to  be 
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dm  =  dc  =:J '.OS  ^9'biM  {Zlfoif  COS  ^  e')cJ&' 
o 

fr 

=  b^(e‘)coi^0'd0 


U-  O,  !  ,2,  •  •  < 


(2.8-26) 


where 


fr 
.  r 


i  =  i  jhi:hB>fc_od^e'j0' 

C) 

fr 

fc  =  J W  c.oi^ e'Je'  j  a  -  i,  £j 


(2.8-27) 


As  can  be  seen  from  the  form  of  the  coefficients >  the  set 
of  equations  for  the  A^^'s  are  completely  Independent  of  the  axial 
coordinate  8.  Once  the  A^'s  are  determined  from  equation  (2.8r24) 
they  are  substituted  back  Into  equation  (2.8-23)  for  determining 
g^*(0)  or  as  easily  g^(0).  The  circulation  distribution  Is 
readily  calculated  from  equation  (2.8-A),  l.e. 


s  S!)  -  go(z)  + 


-  ^(2)  +  g  (e)  cos  (2.8-28) 

In  this  section  the  linearized  flow  field  about  an 
annular  airfoil  has  been  derived  and  has  been  shown  to  be  a 
linear  combination  of  the  airfoil  at  zero  Incidence  and  a  term 
Involving  the  angle  of  attack  but  not  the  section  shape.  It 
should  be  noted  however  that  both  ^erms  are  dependent  on  the 
chord-diameter  ratio  of  the  duct.  Since  the  circulr*:lon  term 
gj^(z)  is  independent  of  the  section  shape,  It  can  be  tabulated 
for  different  angles  of  incidence  and  chord-diameter  ratios. 
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The  presstire  and  velocity  distributions  of  the  annular 
airfoil  at  an  angle  of  incidence  follow  from  Section  II. 7 
and  as  can  be  seen  from  equation  (2.8-28>y  the  effect  of  the 
angle  of  attack  is  to  add  a  term  to  the  coefficient  for  the 
airfoil  at  zero  incidence.  For  instance  the  linearized 
pressure  distribution  follows  from  equation  (2.7-1),  (A-15), 
(B-10)  and  (2.a-2&)  as. 


Rj)] 


f 


Jtr  _ 

I  ^ 

jO  -  4*cof^((p-(p')]^ 

vdiere 


a  £ 


+  2 

[K  J 

f 


(2.8-29) 


induced  velocity  from  vortex  distribution 
of  the  annular  airfoil  at  zero  incidence, 
equation  (A-22). 

Induced  velocity  from  source  distribution, 
equation  (B-10). 


The  integrand  of  the  Integral  has  a  logarithmic  singularity 


at  k  •  1  and  a  square  mot  singularity  at  z'  ■  1. 


To  remove 
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these  first  the  variables  z’  +  cos  9')  and 

z  «(  ^)(1  +  cos  9)  are  introduced,  then  the  variable 
(cos  15  -  cos  O')  -  cos^O, 

AYCUSC^iT) 


(2.8-30) 


where 


1^--^ _ 


The  lift  of  the  annular  airfoil  at  an  angle  of  incidence 

is  given  by  equation  (2.7-6)  and  eqtiatlon  (2.8-28). 
tir  I 

^^,(z')cos<fyzycp 

c 

=  jco/fdcp  j<//' 


/  " 
zr-&rjg,(£')<iz'  =  -BIT 


(2.8-31) 


The  Induced  drag  follows  from  equations (2.7-8),  (2.8-28) 
and  (C-12)  as 


c  _ H 

■  fK 


I  zir 


i  i  ^ 
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,T  tr 

=  -lifgfc&)c0s^0 


-£'«j  (2.8-5J) 


where 


g«  V _ 

^  -cos&‘)‘  +  ‘t 


The  mocent  on  the  annv ' ar  airfoil  about  the  leading  edge 
Is 

®  (2.8-33) 

Introducing  the  Fotirler  aeries  for  ^  and  the  variable 
z  «(?)(!  +  cos  0))  the  moment  can  be  shown  to  depend  only 
on  angle  of  attack  and  not  the  section  shape. 


C/^  =-j  Jg'!{e)s/n^0cojg0c/0 


(2.8-34) 


19 

Some  of  these  coefficient**  have  been  tabulated  by  Welsslnger  . 
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III.  i^ffect  of  the  Hub 

'  Assumiiig  that  the  ordinate  of  the  surface,  of  the  hub  is 
denoted' by and  a  polar  coordinate  system  attached  to  a 
propeller  blade  is  used,  the  linearized  kinematic  boundary 
condition  which  must  be  satisfied  on  the  hub  follows  directly 
from  section  II. 2  as 


z)  =  -  (I  -  ‘XhCzJ  (3_1) 

Vj 


As  before  the  velocities  have  been  nondimensionalized  by 
the  ship  speed  and  the  radial  coordinate  by  the  propeller 
diameter.  The  function  (z)  is  the  slope  of  the  hub  surface 


Figure  4.  Notation  for  the  hub 
Since  the  hub  is  a  symmetrical  body  and  is  assumed  to  have 
no  angle  of  attack  with  respect  to  the  free-stream  velocity, 
the  shape  of  the  hub  can  be  reprasented  bv  a  distribution  of 
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sources,  doublets  or  vortices  over  the  surface  of  the  hub. 

If  the  radial  velocity  can  be  shown  to  be  independent  of  angle 
and  the  hub  is  not  too  blunt,  then  the  flow  can  be  represented 
by  a  distribution  of  sox'rces  and  sinks  along  the  z-axis  only, 
as  indicated  in  Figure  4.  The  velocity  z)  in  equation 

(3-1)  represents  the  whole  flow  field,  i.e.  it  is  composed  of 
the  radial  velocity  induced  by  the  annular  airfoil  with  its 
trailing  vortex  system,  the  propeller  and  the  hub  itself.  In 
terms  of  the  nondimensional  ^ed  velocities  of  the  individual 
singularities  equation  (3-1)  can  be  rewritten  as  follows 

-  -  fJXh,Cp,z)^ 

=  -0 -  ^  (Xk, (pi z)^  (3.2) 

Since  nonrially  it  would  be  expected  that  the  radial 
velocity  induced  by  the  propeller  and  the  annular  airfoil  at 
the  hub  is  small,  it  will  be  assumed  that  the  hub  source 
distribution  satisfies  only  the  average  velocity  at  the  hub. 

If  the  hub  radius  is  zero,  this  must  be  satisfied  exactly. 
Consequently  the  radial  velocities  induced  at  the  hub  will  be 
expressed  as  follows: 


and 
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tTC 

The  use  of  these  induced  \reIocltles  means  essentially 
that  only  the  average  radial  velocity  at  the  hub  iw  considered 
and  that  now  the  hub  can  mathematically  be  represented  by  a 
source  distribution  along  the  z-axis.  Using  equation  (3-3) 
the  boundary  condition  equation  (3-2),  now  becomes 

For  the  derivation  of  the  induced  velocities  from  a 
distribution  of  sources  and  sinks  along  a  line,  consider  first 
a  single  three-dimensional  source  of  strength  qj^(  ^  ')  at  the 
point  f  =  f *»  r=0.  In  polar  coordinates  the  stream  function  of 
such  a  point  source  is^^  (sec  Figure  4) 

I  I  ^  f  ~  ^ - 1 

I  ITf ^5-5) 

The  stream  function  of  an  axisymmctric  body  is  obtained 
by  integrating  a  distribution  of  sources  of  strength  qj^(z') 
per  unit  length  along  the  f  -axis  from  the  after  end  (  f  =bj^) 
to  the  nose  (f  »  b2). 

bi 

The  induced  velocities  follow  from  equation  (2.2-10)  by 
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differentiating  this  equation. 


(3-8) 


If  the  induced  velocity  is  nondiraensionalized  by  the  ship 


speed  (Vg)  the  source  stren  th  by  the  shin  speed  times  the 
propeller  tip  circumference  (27rRpVg),  the  axial  coordinate 
^  ,  by  the  duct  chord  (a)  and  the  radial  coordinates  by  the 

propeller  radius  (Bp),  the  induced  velocities  can  be  written 
in  nondimensionalized  form  as  in  section  II.  TJiis  source 
strength  has  dimensions  of  length^  ner  unit  time  comnared  to 
the  elementary  strength  of  the  ring  vortex  which  has  the 
dimensions  of  length  ner  unit  time. 


=  Jl  (JL)  / _ df' _ 

Vs  hjlJ  {hhHz-Z')*- 


where 


I  _  iu.tt  ehofJ  _ 

duet  JitLmefer  Z^jRp 


at  the  propeller  z=o  and  the  equation  (3-9)  becomes 


(3-9) 


(3-10) 


(3-11) 
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Sometimes  it  may  be  more  convenient  to  nondiraensionalize 
the  axial  coordinate  by  the  propeller  diameter,  i.e.  let 
Zp  =  then  the  axial  induced  velocity  can  be  witten  as 
fol lows 


(3-12) 


The  problem  now  is  to  find  the  source  distribution  qj^(z') 
which  represents  the  shape  of  the  hub.  For  this  the  boundary 
condiffon  given  by  (3-4)  muot  be  used.  From  equation  (3-10) 
the  radial  velocity  at  the  hub  is 

A 


AIza)  I _ ih(i')  dl' _ _ 


and  then  the  boundary  condition  becomes 


A  f. _ dz'  -  (>  -  W^,)'X.{ (z) 


(3-13) 


This  is  a  Fredholm  integral  equation  of  the  first  kind  for 
the  unknown  source  distribution  qh(z).  Direct  inversion  of  this 
type  of  integral  equation  usually  is  not  possible  end  a  solu¬ 
tion  normally  involves  an  infinite  series  of  eigenvalues  of  the 
kernel  and  eigenfunctions.  An  iteration  procedure  for  the  solu¬ 
tion  of  this  type  of  integral  equation  has  been  developed  by 
Lanaweber^^.  In  general,  since  the  diameter  of  the  hub  is 
normally  small  compared  to  its  length  (the  hub  is  often  assumed 


to  be  infinite  in  length),  and  the  velocity  distribution  some 
distance  from  the  hub  and  not  on  the  hub  itself  is  desired, 
further  simplification  can  reasonably  be  made.  The  nrevious 
statement  about  the  effect  of  the  hub  essentially  describes 
the  assumptions  involved  in  slender  body  theory^®  and  it  is 

this  theory  which  will  be  applied  now. 

32 

Laitone  derived  the  slender  body  theory  by  expanding 
the  source  distribution  in  a  Taylor  series.  The  first  term 
of  this  series  gives  the  rt  iilt  of  slender  oody  theory.  Apply¬ 
ing  the  result  from  Reference  [32]  directly  to  equation  (3-13), 
the  followii^g  equation  for  the  source  strength  qj^Cz)  results. 

a  -  -  »vi  J  (Kh  ; ZU  (  3- 14 ) 

where  A(z)  =  Tt  x.h*^(*)»  i.e.  the  hub  cross-sectional  area. 

Equation  (3-14)  shows  that  as  a  first  approximation  the 
source  strength  representing  the  hub  at  a  point  is  a  function 
of  the  change  in  the  cross-sectional  area  of  the  hub  at  that 
point  and  the  radial  velocity  induced  at  the  hub  surface  by 
the  annular  airfoil  and  propeller  blades.  The  velocity  induced 
by  the  hub  at  any  point  in  the  surrounding  flow  field  is 
obtained  by  substituting  equation  (3-12)  into  (3-9)  and  (3-10), 
If  the  hub  is  of  constant  diameter,  then  the  slope  is  zero 


^^(z)  *  ’~(l  "  ^0^  ’~ZhU) 
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^^'(z)=0)  and  the  source  strength  is  a  function  of  onlv  the 

velocities  induced  at  the  hub  by  the  duct  and  propeller.  If 

the  hub  is  assumed  to  be  infinitely  long,  then  bi=- oo  and  b2 

=  in  equations  (3-9)  and  (3-10),  From  a  practical  point  of 

view  this  presents  no  orohlem  since  the  point  for  which  the 

induced  velocity  is  desired  6c, z)  will  be  at  or  close  to  the 

*# 

propeller  and  therefore  at  large  values  of  z  the  integrands  of 
equations  (3-9)  and  (3-10)  become  small  very  rapidly,  as  -p 

for  the  axial  induced  veloci  /  and  1  for  the  radial  induced 

z* 

velocity ,  and  the  integral  can  be  showri  to  converge  uniformly. 
Furthermore  the  slope  of  the  hub  must  be  either  zero  or  undulate 
some  distance  from  the  propeller.  Consequently  the  hub  shape 
some  distance  from  the  propeller  has  no  affect  on  the  flow 
through  the  propeller.  It  would  be  expected  that  in  the  normal 
case  the  velocity  induced  by  the  hub  at  the  duct  would  be  negli¬ 
gible. 
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IV.  Circulation  Theory  of  the  Propeller  in  the  Duct 
IV, 1  Introduction  and  Asaumptions 

The  circulation  theory  as  apolied  to  nronellers  is 
analogous  to  lifting-line  theory  of  vn.ngs  of  finite  length. 

The  main  difference,  which  considerably  comolicates  the  flow 
field,  is  that  the  trailing  vortices  which  were  assumed  to 
lie  lay  in  the  plane  of  the  wing  now  lie  along  helices.  In 
the  rotating  coordinate  system  these  helicoidal  vortex 
sheets  are  stream  surfaces  s  are  the  vortex  sheets  behind 
wings . 

Lerbs  has  developed  the  theory  of  the  moderately  loaded 
propeller^^  and  his  general  approactv  will  be  used  here.* 

In  this  theory  Lerbs  considers  the  influence  of  the  induced 
velocities  on  the  shape  of  the  helical  vortex  sheet  at  the 
lifting  line  but  neglects  the  effects  of  centrifugal  force 
and  of  the  contraction  of  the  slip-stream.  In  addition,  the 
change  in  shape  of  the  vortex  lines  are  neglected  in  the 
axial  direction,  i.e.  they  are  of  constant  pitch.  These  same 
assumptions  will  be  made  here  and  further  it  will  be  assumed 
that  the  influence  of  the  duct  on  che  change  in  the  shape  of 
the  helical  vortex  sheet  in  the  axial  direction  can  be  neglected. 
It  should  he  mentioned  here  that  the  vortex  sheets  are  not 
necessarily  true  helicoidal  surfaces  since  the  pitch  may  vary 
along  the  radius  but  each  vortex  line  is  assumed  to  be  of 
constant  pitch. 


any  of  the  derivations  g 
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In  addition  to  the  effect  of  the  duct  on  the  propeller, 
the  radial  velocities  induced  on  the  cylinder  representing 
the  duct  by  the  propeller  must  also  be  considered.  For  this 
reason  the  induced  radial  velocity  of  the  propeller  must  be 
derived  in  more  general  terms  than  was  done  by  Lerbs, 

In  the  following  development  the  free- stream  velocity 
will  be  allowed  to  have  a  radial  variation  but  must  be 
axisymmetric  and  the  propeller  may  take  any  axial  position 
in  relation  to  the  duct.  Tr„  other  major  assumptions  arc 
stated  more  explicitly  in  Section  I. 

The  boundary  conditions  imposed  on  the  bound  circulation 
is  that  it  be  zero  at  the  hub,^^  P (Rh)=0  and  if  R^>  Rp  then 
the  circulation  at  the  blade  tip  is  zero,  p(Rp)=0,  If  the 
diameter  of  the  duct  and  propeller  are  equal  (R(j-Rp)»  then 
the  circulation  at  the  tip  need  not  be  zero.  This  last 
statement  comes  from  the  fact  that  no  equalization  of  pressure 
takes  place  around  the  blade  tip  if  there  is  no  clearance 
and  the  duct  is  sufficiently  long.  From  a  practical  point 
of  view,  because  of  the  boundary  layer,  no  equalization  of 
pressure  will  take  place  if  the  tip  clearance  is  small.  The 
determination  of  how  small  is  sufficiently  small  requires 
an  analysis  from  boundary  layer  theory.  This  is  not  treated 
here. 
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IV. 2  Induced  Velocities  from  the  Vortex  Lines  of  one  Propeller 
Blade 

The  flow  field  of  each  propeller  blade  is  considered  to 
be  made  up  of  a  system  of  horse-shoe  vortices  lying  along  a 
helix.  The  elementary  system  used  will  consist  of  three  parts; 
a  single  helical  vortex  line,  the  bound  vortex  lying  along  a 
radius,  and  another  free  vortex  line  along  the  negative  ^  -axis. 
This  system  is  shown  in  the  following  figure. 


Figure  6,  Propeller  vortex  system 
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The  point  P,  is  the  point  at  which  the  induced  velocity 
is  desired.  It  is  assumed  that  one  blade  is  located  on  the 
x-axis  and  another  at  the  angle  ^o.  Letting  (b)  be  the  number 
of  blades  the  angle  is  given  by 

dL  ~  -/)Ut  ,  f>  = 

^  b  ^  '  (4.2-1) 

The  location  of  point  P  is  taken  as  arbitrary.  Only  the 
radial  velocity  is  needed  tt  satisfy  the  boundary  condition  on 
the  hub  and  duct  so  only  the  radial  velocity  will  be  calculated 
at  the  arbitrary  point  P  (r  cos  »  r  sinc^,  ^  ).  The  axial 
and  tangential  velocities  induced  by  this  system  are  needed 
only  at  each  blade.  Since  the  blades  are  assumed  to  be  identi¬ 
cal,  it  suffices  to  consider  only  the  blade  along  the  x-axis, 
i.e.  at  the  point  P(x,0,0)  =  P(r,0,0), 

In  Figure  6  the  other  singularities  representing  the  hub 
and  duct  could  be  considered  but  since  the  strength  of  these 
singularities  are  dependent  on  the  induced  velocities  from  the 
propeller  vortex  system  this  procedure  is  not  practical  and  a 
method  of  iteration  must  be  used.  The  other  singularities 
affect  the  hydrodynajnic  pitch  angle  shown  in  this  figure. 

The  velocity  induced  by  a  single  helical  vortex  line 
follows  from  the  Biot-Savart  law,  and  is  derived  in  Appendix  D, 
Designating  this  contribution  to  the  induced  velocity  bv  a  super¬ 
script  (1),  the  radial  velocity  induced  at  an  arbitrary  point  P 
by  a  I.elical  vortex  line  Leaving  from  the  point  F’CrgCOS  djo  .r^sin 
p,0)  follows  from  equation  (D-8)  as 
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H[s>'^(<P-<Pf>  V-  ^cos((p-<Pt>-’<)f^>s^i I ^0, 

'  VTT  J  1  +(z-  t(,oi  f^H/ii}^  -ZKK,  COS(<P  -<f>f  -•<)]*'*  I 

(4.2-2) 

The  axial  and  tangential  velocities  induced  on  the  pro¬ 
peller  blade  at  P(r,0,0)  by  the  same  vortex  line  is  obtained 
from  equations  (D-7)  and  (D-^). 


=  Ejk  / _ ro  -l'cos((p^  +e<) _ ^ 

<*•  \[t^  -il'K,  Cos((pf,  -ht)  -h  ri 


/(■')_  Pn^n/lii  I  I'-to  COSCtf^-f-Uj  -  tj  oi  S/tlCc^ -f-otj 


(4.2-3) 


J  -erncosCfft^)  +  ti‘ +a»/ic  j  (4.2-/') 

The  velocity  induced  by  the  vortex  line  along  the  z-axis 
again  follows  from  the  law  of  Biot-Ravart.  Using  the  desig¬ 
nation  superscript  (2)  the  induced  velocity  is 


r  /  RiXJs  ,  r  /  Unos(p)T  +  (  rsih((>)t  -f-4]X Lot +dt  v-i 

•mOO  ^OO 


O 

_  r  ff(i'sifi4’jr-(t'coi(p)pl  Je  ^ _ r_  (n)i^<P)T -I 

~tfir j  [f-a  +  1  '  tffr  [ 


(4.2-5) 


From  this  equation  it  can  be  seen  that  the  induced  axial 

(2) 

velocity  from  the  vortex  along  the  z-axis  is  zero  i.e,  (w^  =0) 

and  also,  as  would  be  expected,  the  induced  radial  velocity  is 
zero.  This  follows  from  equations  (4.2-5)  and  (A-11). 

(~  -  f'ciPS  <fsin  (p)  -  O 


The  tangential  velocity  follows  from  the  same  set  of 


equations  and  is 


-  rcos^cp)  a  -C. 

^  ^ '  ‘trrt' 

(4.2-6) 

The  induced  velocity  from  the.  radial  vortex  line  (lifting 
line)  is  also  obtained  from  the  Biot-Savart  law  and  is 

*'’1 

JL  f  -  K‘cosip^)T  ^ (ncscp j.  d]Xfcos(&T  j-skAT+o 

li 

\  -  a  j  (4.2-7) 


= - Zl _  _  fe  -  t'cos(<p  -<Pjt]  _  rcosC^'A) 

y*s/n^(f-<p^iil  [vf'  -aT/jcosff 

.  jf-f  •f’ (fees^e)^ j 

^  J(4.2-8) 

The  axial  induced  velocity  of  the  lifting  line  is  given 
by  the  comr  'nent  in  the  ^  direction. 


,(3)_  Ptsin  ( 


H,  -  i'cos((p  -<pfg) 


j(f^  j((f  i»n,  Mao^l.z...)  (4.2-9) 

.or  *  (n— 0,l2.,.)  it  can  be  shown  that  w  ^^^=0. 

At  the  reference  blade  f  =0  and  z=0  and  the  above  equation 
reduces  to 


ra  -fCOS 


(4.2-10) 
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and  if  =  1  nYT,  (n=0,l,2...) 

The  radial  induced  velocity  for  an  arbitrary  point 
P(i-,  (p,z)  is  obtained  from  equation  (4,2-8)  cs 

_ _ 

(4.2-11) 

For  ^  =C  or  (p  =  <Pf  i  «rr.  <’/?  =  <?.  i,^  •“) , 

The  tangential  velocity  at  the  propeller  blade  is  zero 
(w^^^^=0).  This  is  obvious  since  at  the  blade  =0  and  by 
equation  (4.2-8)  both  the  i  and  5^  components  are  then  zero. 

No  singular?.tiee  occur  in  equations  (4,2-10)  and  (4,2-11), 

This  can  be  seen  by  examination  of  equation  (4.2-7), 

The  induced  velocities  at  a  point  P  from  the  three  vortex 
lines  is  obtained  by  summing  and  at  the  blade,  P(r,0,0),  are 

(w^)  =  w  +  w  w  ^^^=0  (4.2-12) 

a  a  ®  a 

where  w_^^^  is  given  by  equation  (4,2-3)  and  is  given  by 

A  A 

equation  (4.2-9) 

(w^)  =  ^  (4.2-13) 

where  is  given  by  equation  (4,2-4)  and  w^^^^  by  equation 

(4.2-6),  '■  . 

The  radial  velocity  induced  at  an  arbitrary  point  P(r  sinq?, 
r  cos  cy’ ,  O  is  given  by 


fe  -  rcOS(<P  -’pft)  ^  KC6i(tP-'ii») 
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(4.2-U) 


where  w^'  is  given  by  equation  (4.2-2)  and  w  is  gi 


given  by 


equation  (4,2-11).. 
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IV. 3  Induced  Velocltiea  from  the  Vortex  Lines  from 
All  the  Blades. 

As  discussed  previously  the  blades  are  assumed  to  be 
evenly  spaced  Td.th  one  vertically  upward  along  the  positive 
x-axis  and  the  helical  vortex  lines  shed  at  a  given  radius 
are  asstased  to  leave  all  the  blades  at  the  same  pitch.  The 
velocity  Induced  at  one  of  the  blades  by  the  horseshoe 
vortices  from  all  the  blades  Is  desired.  It  Is  sufficient 
for  this  purpose  to  consider  the  blade  along  the  vertical 
x-axls  as  the  reference  blade  and  calculate  the  velocities 
induced  at  a  point  F(r,  0,  0)  on  this  blade.  Also  desired 
besides  the  velocity  Induced  at  a  blade  Is  the  radial 
velocity  at  an  arbitrary  point  P(r  cos<^,  r  sln<^,  ^). 

First  consider  the  velocity  components  Induced  by  the 

helical  vortex  lines  shed  from  the  blades  at  the  radius  r  . 

o 

The  total  contribution  to  the  axial  and  tangential  components 
at  any  one  blade  Is  obtained  by  summing  equations  (4.2-3)  and 
(A. 2-4)  over  the  number  of  blades. 

b  r-r  T 


i-oc)  -t 


(4.3-1) 


4.  ¥fr  A  / 


f 


r-'l'g[ct>S((Pfi-oc)  ^Oi)] _ 

-Zi't'o  r.os  t-ot)  f  UhA]  ^ 


(4.3-2) 


The  radial  velocity  Induced  at  an  arbitrary  point 
P(r  co8<f,  r  sin<^,  f)  by  all  the  vortices  shed  from  r^,  is 
obtained  from  equation  (4.2-2)  as 
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(•?  -ff)-’*)  /-  «  COS  (<f  -<fp~  ^<)]  fitrt/ft 


fC^  ~ro<x-t^Lh/3c)^  -Bi^K>cosCf-<f!jo-p<)'J  ' 

(4.3-3) 


The  velocity  Induced  by  the  hub  vortex  ccnslsts  only 
of  the  tangential  component.  Summing  this  component  results 
In  equation  (4.2-6)  being  multiplied  by  the  number  of 
blades k  l.e. 


-i£ 

^  ^  4,  v/rr  (4.3-4) 

Finally  the  velocity  components  Induced  by  the  vortex 

along  the  lifting  line  will  be  discussed.  At  the  reference 

blade  Itself  only  the  axial  component  will  be  considered  as 

on  the  blade  both  the  tangential  and  radial  velocities  are 

zero.  Referring  to  equation  (4.2-10)  It  can  be  shown  that 

there  is  no  contribution  from  this  component  either »  l.e. 

(3) 

(w  '  )t  “0.  In  the  first  place  there  is  no  contribution 
from  the  reference  blade  Itself  nor  from  a  blade  opposite 
It,  l.e.  for  n  »  0  and  1.  The  effect  of  the  other  blades 
must  cancel  in  pairs  since  -  -  ^2’^b  "  1  "  "  ^3» 
Consequently,  since  sin  *^p  “  "  sin  (  -  ^^p)  and  cos<^p 
»  cos  (  -^p),  It  follows  vdien  equation  (4.2-9)  is  summed 

over  the  number  of  blades  that  (Wg^^^)b^=0. 

The  radial  velocity  Induced  at  an  arbitrary  point 
by  the  radial  vortex  lines  equation  (4.2-11)  does  not  cancel. 
This  velocity  follows  by  summing  equation  (4.2-10)  over 
the  number  of  blades. 
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re  -  rcosC<P-<fi>) 

+  K*  +  fs*  -2  rn  cos(<f-^) 

For  f  ■  0  or  <f’  p  +  n  7?',  (n 


0.  1,2, 


..),  (wi^>3 


(4.3-5) 

=0.. 


The  Integrands c£  the  Integrals  for  the  axial  and  tangen¬ 
tial  velocities  on  the  blade,  equations  (4.3-1)  and  (4.3-2), 
are  singular.  To  facilitate  numerical  calculation  of  these 
Integrals  the  "Induction  factors"^^  are  Introduced.  These 
are  nondlmenslonal  qualities  of  the  Induced  velocity 
components  and  are  a  function  of  geometry  only.  They  are 
defined  by 


r 


CLnJ 


(4.3-6) 


From  equations  (4.2-3)  and  (4.2-4)  it  follows  that  the 
axial  and  tangential  Induction  factors  at  the  blades  are 


f'-i. 


4  -  cos  ((ff,  i-et) 


_ Joi. 

[fn-mf  -  2  (a)  ’■‘J 

(4.3-7) 


=(n> 


I  -(^)[cos((f>,-ho()  -ot  -hot)] 


(4.3-8) 
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In  the  limic  as  r 

o 

/a  ‘*(kJ  " 

and 

The  Induction  factors  have  been  tabulated  by  Morgan"^^ 

OO 

and  are  given  in  graphical  form  in  Lerbs’  paper. 

Details  of  the  method  of  calculating  the  induction  factors 

33  35 

are  given  by  Lerbs  and  Wrench.  . 

The  total  velocity  Indi  ?d  by  (b)  vortex  lines  is 

obtained  by  summing  the  different  velocity  components. 

Ihls  summing  does  not  represent  the  total  velocity  induced 

on  a  blade  by  the  system  of  vortices  but  only  from  a  single 

horse-shoe  vortex  on  each  blade.  To  obtain  the  total  velocity 

induced  at  a  point  from  all  the  horse-shoe  vortices  the 

summed  equations  must  be  Integrated  from  the  hub  to  the  tip 

of  the  blade.  The  circulation  strength  T  of  a  single 

trailing  vortex  is  related  to  the  total  circulation  along 

the  blade.  From  wing  theory^^  it  follows  that  the  strength 

of  this  elementary  vortex  is  ^  where  /^(r  )  is  the 

total  strength  of  the  bound  vortex  on  the  blade.  With  this 

notation  it  then  follows  from  equation  (4.3-6)  that  the  axial 

velocity  Induced  by  the  propeller  at  one  of  its  blades  is 


and  from  equations  (4.3-4)  and  (4.3-6) 


(4.3-9) 
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'-'r  ‘trldro 


hlD. 


O'-  n) 


f  i  Arc 


=  - f  bPO?.)  .  /  im  ii  (^)  ,J 
,ir[  dir 


(4.3-10) 


If  there  Is  clearance  between  the  blade  tip  and  duct, 
then  as  discussed  in  sectlo  (IV. 1),  the  circulation  at  the 
tip  is  zero,  r(R^)  »  0,  and  the  first  term  does  not  exist 
For  convenience  both  of  these  equations  will  be  non- 
dimensionallzed;  let 


X  • 

reference  radius 

propeller  radius 

X  » 

la  = 

radius  at  which  vortex  is  shed 

o 

propeller  radius 

\  - 

hub  radius 

propeller  radius 

and  for  the  circulation 


,  r 


(I  -  ^r)G 


(I  -  Wr)r 
SrrR/fiVfCt) 


(4.3-11) 

If  the  free-stream  velocity  is  constant,  then  the  free-stream 
velocity  w^  is  used  for  nondime ns ionalizfng  instead  of  the 

ship  speed  V  . 

8 


Introducing  this  notation  into  equations  (4^-9)  and 
(4.^^-10)  the  induced  velocities  are  given  in  nondimenslonal  fo?m. 
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/ 

(X-Xc) 


CtL.  JXo 


and 


(4.3-12) 


or 


(4.3-iaf) 


The  Integrand  of  the  Integrals  In  these  equations  are 
singular  at  x  Xq.  Lerbs33  has  discussed  both  these  equations 
except  that  the  equations  were  derived  for  a  propeller  without 
the  duct  so  that  6^(1)  -  0.  This  term  cAuses  no  difficulty, 
however,  If  It  Is  rewritten  on  the  left-hand  side  as  shown 
since  Lerbs  deals  with  the  same  Integral. 
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IV. 4  Thrust  and  Torque  Developed  by  the  Propeller 

The  thrust  on  the  propeller  follows  from  the  law  of 
Kutta- Joukowski t  equation  (2.7"4).  For  each  element  of 
each  blade  this  law  states 


d  T  V njr  =  /O  rM  ‘  (4.4-1) 

where 

2  Tf  (rps)  Is  e  angular  velocity  of  the 
propeller,  w^.^  Is  the  tangential  velocity 
Induced  at  the  propeller  blade  by  the  duct  and 
Includes  the  tangential  induced  velocities  from 
the  ring  vortices,  ring  sources  and  duct  trailing 
vortex  system,  w^^  is  the  tangential  Induced 
velocity  by  the  trailing  vortex  system  of  the 
propeller,  equation  (4.3-13). 

For  the  thrust  the  velocity  V  is  the  total  tangential 

velocity  by  the  propeller  llftlng-line  excluding  self-induced 

velocities.  It  should  be  noted  that  by  assumption  In  Section 

III  the  tangential  induced  velocities  from  the  hub  are  zero 

(3) 

and  for  a  lifting  line  the  self  Induced  velocities  (w^  ) 

arc  zero.  The  total  thrust  Is  obtained  by  Integi.'atlng  this 
equation  from  the  propeller  hub  to  the  tip  and  sunomlng  over 
the  number  of  blades.  Since  the  thrust  of  each  blade  Is  the 
same,  the  thrust  of  one  blade  Is  multiplied  by  the  nimtber 
of  blades. 


81 


T^-  =/shj  r'  -  wij(f')  -  (rjjd't' 


(4.4-2) 


This  equation  Is  ncndlmenslonallzed  using  the  same 
notation  as  for  equation  (4.3-11)  with  the  added  definition 
of  the  advance  coefficient 


iOoRp 


(4.4-3) 


and  the  thrust  coefficient  of  the  propeller  ic 

I 

-(f*- 1')]“'* 


(4.4-A) 


If  the  free  stream  velocity  is  uniform  over  the 
radius  then  it  is  used  for  nondimensionalizing  instead  of 
the  ship  speed.  The  subscript  "i"  is  used  to  denote  that 
the  thrust  is  the  thrust  in  an  inviscid  fluid.  The  total 
thrust  of  the  ducted  propeller  is  given  by  adding  this 
equation  to  equation  (2,7-9), 

The  torque  is  also  obtained  using  the  law  of  Kutta- 
Joulcowski,  This  law  gives  a  tangential  force  at  each  radius 
which  when  multiplied  by  the  radius  and  integrated  over  the 
blade  length  gives  the  torque  per  blade. 


where 


w  (r)  is  the  free-stream  velocity 
o 

is  the  axial  velocity  induced  by  the  hub 
Wfld  axial  velocity  induced  by  the  duct 


Wgp  is  the  axial  velocity  induced  by  the  trailing 

vortex  system  of  the  propeller,  equation  (4.3~9) 

Integrating  this  equation  from  the  hub  to  the  blade  tip 

and  multiplying  by  the  number  of  blades  gives  the  total 

torque  of  the  propeller. 


r 

ci=  byC  I  rr(i')[v/^M  1 1^4  +  i^apje/r 


"  (4.4-5) 

Nondlmensionalizlng  as  before  but  defining  a  power  coefficient 


r  —  ^  ^ 

"f 

the  nondlmenslonalized  fozm  of  this  equation  is 
> 

where  (1"W„)  wake. 


(4.4-6) 


Tne  ideal  efficiency  (inviscid)  is  given  by  the  ratio 
of  the  total  thrust  to  the  power,  i.e. 


=  (EmJc/  + 

^psi 


(4.4-7) 


where  (C(.gj^)d  is  the  duct  thrust  given  by  equation  (2.7-9) 

*tip  "  «  ideal  efficiency  of  the  propeller 
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IV. 5  'Che  Integral  Equation  for  the  Circulation 


The  free  helical  vortices  lie  on  stream  surfaces  and 
are  assisned  to  have  the  pitch  of  the  resultant  flow  angle 
at  the  lifting  line.  The  velocities  at  a  radius  x  are  shot<m 
In  the  following  velocity  diagram. 


-^je)  S  V/,(X) 


Figure  1.  Velocity  diagram  of  the  ducted  propeller 
The  helical  vortex  is  shed  at  the  angle  which  Is 
commonly  called  the  hydrodynamic  pitch  angle.  The  angle 
is  the  propeller  advance  angle. 

From  this  velocity  diagram  the  following  equation  Is 
obtained  for 

wkj -f- 


(A. 5-1) 
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Substituting  in  for  and  ,  equations 

(4.3-12)  and  4.3-13)  this  equation  is  obtained  in  terms  of 
the  unknown  circulation 


fa.yy3i  = 


or 


Cfikn^i 


+•  ifi. 


dXo 


% 


^j(l) 

zx 


-0  -  ^r) 


(4.5-2) 


with  boundary  conditions 

C  o  i'f  Rfi  ^  Rj 

Gt(Xk)  =  0  And  Sj(l) 

This  last  equation  represents  an  integro-differential 
equation  for  the  unknown  circulation  distribution  which  will 
give  the  desired  thrust.  Lerbs^^  gives  a  nvnnerical  method 
for  solving  this  equation  so  only  a  few  remarks  will  be 
made  concerning  it.  The  left  hand  side  of  this  equation  is 
the  same  as  Lerbs  discussed  but  the  right  hand  side  contains 
velocities  induced  by  the  hub  and  duct  and,  in  addition,  the 
circulation  at  the  tip,  if  it  is  not  zero.  Since  the  right 
hand  side  is  assvnned  known  in  either  the  free-running  or 
ducted  propeller  case,  this  does  not  affect  the  solution  method. 
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In  s  strict  sense  the  induced  velocities  of  the  duct  and 
hub  at  the  propeller  depend  on  the  propeller  circulation  but 
because  of  the  complexity  of  the  problem  these  velocities 
must  be  assumed  known.  To  consider  that  they  are  dependent 
on  the  propeller  circulation  in  equation  (4.5-’2)  implies 
that  the  circulation  distribution  representing  the  nozzle, 
the  source  distribution  representing  the  hub,  and  the  propeller 
blade  circulation  must  be  determined  simultaneously.  This 
of  course  is  not  possible  '>nd  it  is  necessary  to  resort  to  a 
method  of  Iteration.  For  Instance,  the  duct  and  hub  problem 
could  be  solved  separately  without  the  propeller  and  then 
using  the  resulting  Induced  velocities  the  propeller 
problem  solved.  This  process  is  then  repeated,  using  each 
time  the  last  derived  induced  velocities,  until  satisfactory 
convergence  is  obtained. 

Equation  (4.5-2)  is  in  a  general  form  as  the  free-stream 

velocity  may  vary  radially  (wake-adapted'  propellers)  and  the 

circulation  distribution  need  not  be  optimum.  In  addition 

this  equation  applies  to  moderately  loaded  propellers*  as 

well  as,  of  course,  to  lightly  loaded  propellers.  These 

33 

various  cases  are  discussed  by  Lerbs.*^ 


*  Tne  difference  between  a  moderately  loaded  and  a  lightly 
loaded  propeller  is  that  for  a  moderately  loaded  propeller 
the  velocities  induced  by  the  helical  vortices  are 
Included  in  the  calculation  of  the  angle  at  vdiich  these 
vortices  are  shed  while  for  a  lightly  loaded  propeller 
their  effect  is  ignored. 


IV . 6  The  Optlminn  Circulation  Distribution  of  the  Propeller 


iU 


For  Che  free'-runnlng  propeller  a  question  arises  as  to 
what  is  the  circulation  distribution  so  that  a  propeller 
produces  a  given  thrust  with  a  minimum  aniount  of  power.  This 
optimum  circulation  distribution  Is,  of  course,  based  on  Che 
llftlng-llne  theory  and  an  Invlsc id  fluid.  For  the  ducted 
propeller  a  similar  question  arises  but  the  force  on  Che 
duct  Itself  enters  the  problem.  The  problem  could  also 
Include  determination  of  th.  shape  of  the  duct  as  well  as 
the  propeller  circulation  distribution.  Dlckmann  and 
Welsslnger^®  have  considered  this  problem,  the  optimum 
shape  of  Che  duct,  but  for  ducts  of  zero  thickness  anv.  u 
simplified  representation  of  the  propeller. 

The  combined  problem  of  optimum  duct  shape  and  optimum 
circulation  distribution  along  Che  propeller  blade  Is  a 
formidable  one  since  It  Is  not  possible,  within  the  concepts 
of  the  theory  developed  here,  to  obtain  the  interference 
velocities  In  explicit  form.  This  can  be  seen  by  referring 
to  equation  (^.5-2)  in  which  it  Is  necessary,  in  order  to 
obtain  a  solution  to  assume  that  the  Induced  velocities  from 
the  nozzle  and  hub  are  known  and  not  functions  of  the  circula¬ 
tion.  For  the  same  reason  It  Is  not  feasible  to  take  Into 
consideration  the  total  thrust  of  the  ducted  propeller  system 
but  to  consider  only  the  propeller  thrust. 

The  problem,  which  can  reasonably  be  solved, reduces  to 
the  determination  of  the  circulation  distribution  on  the 
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propeller  blade  so  that  the  propeller  produces  a  given 
thrust  with  minimum  torque.  This  approach  doss  not 
consider  whether  or  not  for  this  optimum  circulation 
distribution  the  complete  ducted  system  produces  a  given 
amount  of  thrust  with  minimum  torque.  For  this  reason  the 
problem  posed  is  somewhat  academic  and  is  discussed  further 
only  for  the  sake  of  completeness. 

In  the  following  analysis  it  will  be  assumed  that  the 
circulation  at  the  blade  tl  is  zero,  Gg(l)  *  0,  and  that 
the  frce-stream  velocity  is  a  constant,  Wq(x)  >  Wq  or 
(1  *■  w^^)  •  1.  This  essentially  reduces  the  problem  to  a 
free -running ,  moderately  loaded  propeller  in  uniform  flow 
with  the  addition  of  velocities  Induced  by  the  duct  and  hub. 
Nondlmensionallzing  with  the  free-stream  velocity,  the  thrust 
coefficient,  equation  (4.4-4),  and  power  coefficient,  equation 
(4.4-5)  become 


and 


where 


(4.6-1) 


(4.6-2) 


r - - -JJlf-xG  M 


A  = 


si  Rf> 


G  = 


n 

ZtrRfiK 
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The  propeller  induced  velocities,  and  w. ,  in  these 
equations  can  be  replaced  by  their  values  in  terms  of  the 
circulation  distribution,  equations  (4.3-12)  and  (4.3-13). 
After  Integration  by  parts  the  thrust  and  power  coefficients 
can  be  written  in  terms  of  the  derivative  of  the  circulation 
distribution. 


(4.6-4) 

In  this  form  the  thrust  and  power  coefficients  are 
fxinctions  of  G*(x)  and  not  of  both  G'(x)  and  G(x).  The 
problem  is  to  find  the  circulation  distribution  G(x)  so  that 
the  power  coefficient  (Cp^)  is  a  minimum  while  the  thrust 
coefficient  remains  unchanged.  This  is  a  problem  in 

the  calculus  of  variations. 

A  small  variation  is  now  taken  of  the  slope  of  the 
circulation  distribution,  i.e.  G’(x)  +  cT G'(x)  represents 
the  value  of  the  slope  of  the  circulation  distribution  in  a 
small  region  surrounding  the  point  x  and  at  the  end  points  the 
variation  is  zero,  /g'(xjj)  •«</g*(1)  -  0.  If  the  circulation 
distribution  is  an  optimum,  then  in  this  small  region 
CTi(G'(x)  +<^G’(x)jp  and  C-j-£[G'(x)]  have  the  same  value. 

If  second  order  terms  are  ignored,  the  difference 
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CrjIG-W  +  -  CrjfG-(x))^  l,  obtained  as  follows; 


/  I 


(<fc^  =  C..[G'(x)^fs'rx)]^- Cr.[G'(x)]^=8bj  J gWx 

ri  r  r  '* 

Xh  ih  n 

i  H,  i 


I  ^ 


/G'Cryx 


=  0 


(4.6-5) 

ihe  variation  (cfCj^)^  is  taken  as  zero  since  within 
the  variation  In  G'(x).  OtJG'(x)  H-*-(x)lp-  C^i[G-(x)]  Is 
taken  as  aero.  A  first  order  variation  /c  ,  may  also  be 


defined  as 


1 

II 

1 — 

"4? 

_ 

-^XA  J 

-f 


XA 


-fj 

XA 


xa 


+  -^djx'c/x'  -  lUx,) 


r.  X 


^  ^  *^'+  2  ^  '(Xol 


XA  -A), 

'  2^ 

X' 


XA 


x> 


Xa 


jt/x»  JG'Mchf 

; 


=  0 


(4.6-6) 

The  variation  (fc^^  must  be  approximately  sero  for  a 
sufficiently  small  variation  in  order  for  Cpj  to  be  a  mlntamn. 
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Actually  by  this  analysis  can  be  a  maximum  or  a  mlnlmun 

and  to  show  that  C  ^  for  the  optimum  value  of  G'(x)  Is  a 

pi 

minimum  it  should  be  shown  that  Cpj^[ /G' (x)]  >  0  for 

cff  0.  This  is  not  possible  without  knowing  the  form  of  G'(x), 
however,  it  would  be  expected  that  Is  a  tiilnlmum. 

By  equations  (4.6-5)  and  (4.6-6)  both  ^^d  <fC^f 

are  zero  for  all  (/g'(x)  (  cf  sufficiently  small)  and  this 
can  be  true  if  and  only  If  the  Integrands  of  both  these 
equations  are  proportional..  The  equation  that  the  circula¬ 


tion  distribution  must  satisfy  is  then  the  following: 


(4.6-8) 


This  is  an  Integral  equation  of  the  first  kind  for  the 
circulation  distribution  in  terms  of  the  constant  A.  As 
discussed  in  section  III  direct  Inversion  of  an  integral 
of  this  type  is  usually  not  possible,  however  Landweber 
has  discussed  an  iteration  procedure  for  this  type  of  Integral. 
The  solution  of  this  eviuatlon  will  not  be  discussed  further. 
Once  the  form  of  the  circulation  is  obtained  then  the  value 


1 


“1 

I 
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of  the  parameter  "A"  can  be  calculated  from  the  eqtiatlon  for 
the  thrust  coefficient  equation  (A. 6-3),  and  the  value  will 
depend  on  the  value  of  the  thrust  coefficient. 

If  the  propeller  Is  a  free-running  propeller,  then 

- P  and  -  are  zero  and  this  equation  (4.6-8) 

reduces  to  that  for  the  optimum  circulation  distribution  for 
a  moderately  loaded  propel  r*.  (The  velocity  induced  by  the 
hub  Is  usxially  neglected.) 

By  interchanging  the  order  of  Integration  and  then  integrating 
by  parts, equation  (4.6-7)  can  be  written  in  the  following  form: 


After  differentiating  and  Inverting,  this  equation  becomes 


If  the  Integrals  are  zero  in  this  equation,  then  from 

Figure  7  it  can  be  seen  that  this  Implies  x  >  tan/?  1  -  constant. 

This  is  mentioned  in  light  of  Betz's  theorem  which  states  that 

*  This  particular  form  was  obtained  by  Prof.  J.V.  Wehausen  in 
his  unpublished  class  notes  for  "Hydrod}rnamics  of  Ships". 
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for  a  fr«e-runnlng  llijhtly  loadvd  propeller  che  circulation 
dletrlbutlon  i$  optimum  if  the  «ngle  eetieflee 

X  tan/^i  "  constant.  Zt  ie  not  obvious  that  equation  (4.6-9) 
will  result  in  thlH  theorent  even  if  the  ]>ropcll«r  is 
assumed  to  be  free-nmning  and  lightly  loaded. 


93 


V.  Interaction  Effects 

When  discussing  the  duct  and  hub  in  Sections  II  and  III, 
the  form  of  the  propellei-  induced  velocities  was  not  included 
because  they  had  not  yet  been  derived,  'fhia  section  then  will 
deal  mainly  with  the  effect  of  the  propeller  induced  velocities 
on  the  duct,  however,  the  equations  for  the  duct  and  propeller 
induced  velocities  at  the  hub  will  be  given  in  a  more  explicit 
form  than  in  Section  III, 

V,1  The  Duct  Trailing  Vortex  System 

As  stated  in  Section  II,  when  an  annular  airfoil  is  sub¬ 
jected  to  a  radial  velocity  which  is  dependent  on  the  angular 
position,  a  trailing  vortex  of  strength  ^  is  shed  from 

each  point  (<P,z)  on  the  duct.  This  trailing  vortex  system  is 
shed  at  an  angle  equal  to  the  flow  angle  and  follows  a  stream 
line  in  the  rotating  coordinate  system.  This  implies,  as  in 
case  of  the  propeller,  that  the  induced  velocities  from  all 
the  components  in  the  flow  field  have  an  effect  on  the  trailing 
vortex  system.  This  represents  a  problem  considerably  more 
difficult  than  the  propeller  problem  since  the  helical  vortices 
arc  shed  from  all  over  the  duct  rather  than  along  a  line.  To 
obtain  the  equation  in  a  form  which  is  amenable  to  solution 
and  yet  which,  it  is  felt,  represents  the  flow  field  rather 
well,  it  will  be  assumed  that  the  helical  vortices  arc  all 
shed  at  the  advance  angle  of  the  duct.  This  is  the  angle 


94 


given  by  the  following  equation 

thu/SJ  -  '  u^^Hj 


(5.1-1) 


It  will  further  be  assumed  that  these  vortices  maintain  a 
constant  pitch  angle  and  form  a  cylindrical  vortex  sheet  of 
diameter  extending  from  the  duct  to  minus  ini'inity  in  the 
axial  direction. 

The  velocity  induced  l'>j'  a  single  helical  vortex  line  is 
derived  in  Appendix  D.  The  v  ^ocity  comnrnents  induced  by 
the  cylindrical  vortex  sheet  shed  from  tne  duct  are  obtained 
by  integrating  equation  (D-7),  (D-8)  ana  (D-9)  over  the  sur¬ 
face  of  the  duct.  Each  vortex  shed  from  the  duct  has  a  radius 
Rj,  pitch  angle  and  strength  ^  ^  so  these  comnonents  can 
be  written  in  the  following  nondimensionalized  form 

t  err  0  . 

"  2rJ  J  I  (5.1-2) 


-  ±.f  /I ^  UJ  - 


(5'.  1-3)^ 


jr  ill  R3 

(5.1-4) 


where 


/?*  =  /  +[zh(2-Z')  ~o<tan4if 
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At  tlie  oropeller  reference  blade  z=0  (i.e.  and 

^  =0,  so  the  above  equations  are  considerably  simplified 
at  that  point.  In  order  to  obtain  the  circulation  distri¬ 
bution  on  the  duct  it  is  necessary  to  have  the  radial  velocity 
induced  on  the  duct  by  the  trailing  vortex  system.  This  is 
given  by  letting  x=Xjj  in  equation  (5.1-3). 


o  o 


l2h(s~z'^ 

[  fajj  ■  +[2h(f-z') f 

(5.1-5) 


The  uniform  convergence  of  the  infinite  integral  in  this 
equation  is  discussed  in  Appendix  D, 
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V . 2  The  Radial  Velocities  Induced  at  the  Duct  by  the  Propeller 
and  Hub 

The  radial  velocity  induced  on  the  duct  surface  by  the 
nrooeller  is  given  by  integrating  equations  (A. 3-5)  and  (A. 3-5) 
along  the  radius  and  sianming  the  two  equations.  The  cii'cu- 
lation  denoted  by  P  for  a  single  vortex  when  referred  to  the 
circulation  at  the  lifting  line  is  •  Substituting 

this  for  P  in  equation  (A. 3-3)  and  integrating  from  the  hub 
to  the  tip,  the  trailing  voi  ax  system  induces  the  following 
radial  velodity  on  the  duct. 


^  stn((p- J 

/  j\ f/?j  +  r/-h(^- rpotUh/Sc)^  -Z Rj roCos(<f~^ 

f“Loa 


Joicjli 


If  this  equation  is  nondimensionalized  as  previously, 
equation  (Ac3-10)j  then  it  can  be  written  in  the  following 
form. 


(5.2-1) 


r-<0j 

(■r 


fan/it  sin 


(5.2-2) 


As  discussed  in  Apnendix  D  the  integrand  has  a  finite 
jump  discontinuity  v^hen  the  point  P(.x,  z)  lies  on  the  helix 
and  z^  or  <f  4  ^  . 
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This  can  only  occur  in  the  foregoing  equation  when  the  duct 
has  the  same  diameter  as  the  oiopeller  (x^j-i-Xg)  and  2hz= 

taa/3i  .  At  this  point  the  integrand  has  a  finite  jump 

since 


-  + 


't'a.r./ii 

3(l  ‘ 


Z40  or 


If  z=0  and<?=<^p,  then  the  integraind  is  singular  as  dis¬ 
cussed  in  Appendix  D.  For  calculating  purooses  it  is  nrobably 
best  to  nondimensionalize  the  ixial  comnonent  with  the  propeller 
radius,  since  this  equation  is  only  dependent  on  the  duct  chord 
through  the  non-dimensionalization  of  the  axial  component.  In 

this  form  then  equation  (5.2-2)  is 
b  ° 


where  Zp  is  the  axial  coordinate  nondimensionalized  by  the 
propeller  radius  Rp. 

This  equation  is  independent  of  the  duct  chord  and  once 

tabulated  on  the  basis  of  Zp  the  values  can  easily  be  changed 

to  those  for  the  axial  coordinate  z  by  z=  .  It  should  be 

noted  that  this  factor  ij.  is  not  the  normal  induction  factor^^ 

i^  which  is  i  =(l-*o)  i  .  If  the  velocity  induced  at  the  hub 
”  ^  XT'.  ^ 

by  trailing  vortex  svstera  xj^  is  introduced  in  the  foregoing 
equations  in  place  of  x^j,  then  this  equation  gives  the  radial 
velocity  at  the  hub. 

The  radial  velocity  induced  by  the  line  vortex  itself  is 
obtained  from  equation  (4.2-10),  In  nondimensionalized  form 
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this  velocity  is^ 
Vs 


^ _ 


.  +  £0^ 


Uxo 


(•5.2-4) 


For  Zp=o  or  ^=(j(?p  , (n=0,l,2. . . ) 

This  equation  is  also  written  independently  of  the  duct 
chord  by  nondimensionalizing  wit-h  the  propeller  radius,  and 
as  before  Zp=2hx^^z.  This  equation  is  in  a  general  form  and 
can  be  applied  at  the  oroDeller  hub  or  the  duct  by  using 
either  the  hub  radius  or  duct  radius  instead  of  x  in 
the  above  equation.  If  a  factor  defined  as  for  the 

free,  vortex  sheets,  then  the.  above  equation  can  be  written  in 


a  more  simplified  form 


.f-V 


/ 

Xh 


(3. .2-5) 


where 


T%  a,  1  /Vf  \  (¥1  - 


f=l 

(5.2-6) 

And  combining  this  equation  with  equation  (5.2-5),  tne,  total 

radial  induced  velocity  by  the  nroneller  is  given  as 

I 

“'"fr/'tf ; inz)  - t!.%','P,z)]jx, 


^  n, 


(5.2-7) 


or  at  the  duct  x=Xj[ 


V$ 


fG/%) 


^(3h 


(i%6 


Xh 


99 


=  ^ J” Qs(^o)(ctf)^ciXo 


(5.2-8) 


For  any  ducted  svstem  the  circulation  distribution  Gg  is 
a  function  of  the  duct  radius  the  nitch  angle  /Sc,  the  number 
of  blades  b  and  the  duct  and  hub  shane.  It  is  written  here  as 
only  a  function  of  the  radiu  since  for  any  one  configuration 
it  can  only  be  a  function  of  the  radius.  It  should  be  noted 
that  the  factor  (ir)p  contains  both  odd  and  even  terms. 

The  radial  velocity  induced  at  the  duct  by  the  hub  is 


obtained  from  equation  (3-10)  by  letting  x=Xjj 


(5.2-9) 


The  function  q(z’)  is  the  source  strength  of  the  line 
source  representing  the  hub  and  is  given  by  equation  (3-14) 
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V. 3  Integral  Equation  for  the  Circulation  Distribution  of 
the  Duct  Ring  Vortices. 

It  was  shown  in  Section  11,3  that  the  ring  siource  distri¬ 
bution  was  a  function  of  the  aiinular  airfoil  thickness  only. 
Using  this  fact  and  substituting  into  eqimtion  (2.3-3)  the 
radial  velocities  induced  at  the  duct  by  the  duct  trailing 
vortex  sheet,  propeller,  and  hub,  equations  (5.1-5,  5.2-7 
and  5.2-9),  the  integral  eqi  tion  for  the  vortex  distribution 

is  obtained  as 
/  rr 

O  -IT 

t  'ir  o  » 

^Lf{u.(<fiA  rrr2A>(/-/'j -^un/3j]coi«>-<p'-oi) - -^wi 
JJ<^  +[ih(X‘Z')  J  I 


where 


'  1 
U(<Pj^)  BH'[c'(l)  +tati<x]  +■  Sh  [ sW)-4[K/JI)  -EWJcJz* 

i  J  i 


(5.3-2) 


This  equation  io  a  singular  integro-diffcrential  equation 
in  two-dimensions  for  the  circulation  distribution.  The  oart 
of  the  equation  to  the  left  ■  of  the  equal  sign  is  assumed 

:  J-  . 

known,  i.e.  the  duct  shape  and  the  radial  velocities  induced 
by  the  piopeller  and  hub.  It  should  be  mentioned  that  it 
wuld  be  exoected  with  the  normal  configuration  the  radial' 
velocity  induced  by  the  hub  at  the  duct  would  be  negligible. 
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The  integral  equation  can  be  reduced  to  a  one-dimensional 
one  if  it  is  assumed  that  the  ring  vortex  strength  can  be 
expanded  in  a  Fourier  Series  in  (f  i.e. 


ao  oo 

Yi-O  «  =  / 

ir 

yz)  =  i-jrl 

ft' 

3nC-^)  *•  ~j  r((p,  Z)  CuSh  cpje^ 

-ft 

ir 

hJZ)=^j)!’(‘PjZ)sinr>(f>clf 


(5.3-3) 


(5.3-4) 


For  convenience  the  oart  of  the  integral  equation  denoted 
by  U(<p,2)  will  also  be  expanded  in  a  Fourier  series  in  . 
This  involves  no  assumotions  on  the  form  of  U(<^,z)  since  all 
the  functions  in  U(<?^,z)  are  continuous  with  resnect  to  . 


Ij(cf^z)  -  (zjcos/i  (f  ■ry^vA(z)s}f^n(f 


(5.3-5) 


If 

Uc(Z)  = 

TT 

U„(Z)  =  jpju((f,z)cosnCpcj(P 
■r 

Vn  fz)  =  :^J  U((P,  z)  5  in  n  {p<i(P 


(5.3-6) 
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Fi’om  equation  (5.3-2)  it  is  apparent  that  u^(z)  and  Vj^<z) 
for  n  S  1  will  conLain  terms  involving  only  the  propeller  in¬ 
duced  velocities.  The  Fourier  coefficient  u^Cz)  is  a  function 
of  the  duct  shane,  radial  induced  velocities  frota  the  hub  and 
the  average  radial  induced  velocity  from  the  oropeller.  It  can 
immediately  be  xirritten  as 


Uo(z)  = 


! 

-  2  (/  -  v^t/) |‘fr[c'(zi  t  j  +  hj s'W)-4[K('^)  -  ECJil  Jc/ij 


(5.3-7) 


The  function  H(z)  is  the  same  as  that  given  bv  equation 
(2.A-3).  Since  the  integrals  in  this  equation  are  functions 
of  hub  shaoe  and  circulation  distribution  on  the  propeller 
lifting  line  and  are  assumed  known  for  the  duct  problem  they 
can  be  evaluated.  The  integral  dependent  on  the  hub  shape 
is  simple  enough  that  it  can  easily  be  solved  numerically,  if 
necessary,  and  the  integral  involved  the  blade  lifting  line 
can  be  reduced  to  a  more  simple  form  for  numerical  solution. 
Consider  the  following  integral 
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I  IT  I  ir 

-  I  iUx.)  Jn, 

”•  ‘f  (5.3-S) 

The  integral  involving  the  lifting  line  can  easily  be  shown 
to  be  zero.  From  equation  (5.2-5) 


The  integral  involving  the  helical  vortex  sheet  can  be 
reduced  to  a  function  of  elliptic  integrals.  From  equation 
(5.2-2)  rr  f 


The  infinite  integral  in  this  equation  is  uniformly  con¬ 
vergent  with  respect  to  and  Kq  if  zio  and  consequently 

the  order  of  integration  can  be  interchanged  (see  Annendices  D 
and  E).  By  changing  the  order  of  integration,  the  integration 
with  respect  to  <><•  and  ,  can  be  carried  out. 
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(5.3-11) 

The  change  in  variable  <^=  -c(  has  been  made  in  this 
equation  and  also  the  angle  of  the  blade  does  not  appear 
since  it  maltes  no  difference  'hat  comnlete  cvcle  of  a  periodic 
function  is  integrated.  The  integral  in  the  brackets  is  dis¬ 
cussed  in  Appendix  E.  The  reduced  form  is  obtained  by  taking 
the  value  of  equation  (E-18)  at  n=0.  This  gives 


where 


^  f cose^ 

0 


(5.5-12) 


After  integrating  part  of  the  integrand,  this  integral  has 


the  following  form. 


(5.3-13) 
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where 

The  double  integral  given  by  equation  (5,3-13)  does  not 
exist  for  k=l,  i,e.  x^j=Xjj=l  and  z=0.  The  asymptotic  behavior 
is  siich  that  the  integral  is  infinite  as  l,ij  oo  at  this  point 
which  is  also  the  value  obtained  after  interchanging  the 
order  of  integration,  ^^en  k=0  the  value  of  the  double  inte¬ 
gral  is  zero  since^^ 


lim  (K-  E  1- 


=0 


The  singular! tv  at  k=l  of  the  right-hand  side  of  equation 
(5,3-13)  is  due  to  the  elliptic  integral  of  the  second  kind 
K(k),  however  this  causes  little  difficultv  since  this  singu¬ 
larity  occurs  under  the  integral  with  respect  to  and  a 
logarithmic  singularity  integrates  out.  If  equation  (5,3-13) 
is  introduced  into  equation  (5,3-8)  tmd  a  change  of  variable 
is  made  the  logarithmic  singularity  can  be  removed. 

Assuming  x^=l  and  z=Q  then  let  (1- 


jm 


o 
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U07 


The  factor  is  given  by  equation  (5.3-12)  however  when 
z=0  and  x^=l,  equation  (5.3-14)  should  be  used  for  this 
equation,  i.c. 


.3-17 

where 

^2  = 

(2-t3)2 

It  is  obvious  by  the  form  of  U(f  ,z),  equation  (5.3-2), 
and  the  above  form  of  Uo(*)  that  the  Fourier  coefficients  of 
higher  order,  g^(z)  and  h^(z)  for  n  i  1,  are  functions  only 
of  the  propeller  induced  velocity.  Prom  equation  (5.2-8) 
and  equation  (5.3-6)  the  Fourier  coefficients  and  Vn(z;) 

then  follofa  as 

ir  I 

(t]  s  -  ij(cosn(P)U6'CY.)[c^  fg ; cf, Z)]  Jxo]  Jf 

(5.3-18) 

and 


Uo(o)  =  -2H(0j 


../T  fht2'}(iz' 


3  b  1^0 [(^  -Kwv  -2c(^»)p 

,  -ft 


VnU) 


V  / 

5<A 


and  substituting  for  (i^.)^  from  equation  (5.2-8), 


(5.3-19) 


Tloa 


/  ir 


Uf(z}  =  -  (to)  J (cosnrp)[ 


(■5.3-20) 


7h 


I  'I 

,/„(z)  =-if OUx»)J  (&inn<(’)TZl‘^-Zy^^]AceJ?^i 


f c_^-OI) 


X/i  -ff 


The  integral  for  equation  (5.3-20)  will  be  dis¬ 

cussed  first.  ■  Substituting  in  from  equation  (5.2-6)  this  can 

be  written 

7  ^,jj 
J(£{)Sn^)Tr  Jcf 

ir 

fcosncfl  [  (k)  cosK-_^ 


(%)  _  cos(4>-<&)]\ 

f  -zmc»s(^<fp)  1  ]l  I 


The  equation  has  the  value  zero  since 


for  any  number  of  blades. 


(5-3-23) 


If  the  Fourier  coefficients  v^(z)  arc  being  considered 

then 


■’1X)9 


where 


b  if  n  =  *>*6 


fzi  |o  if  n=.>»k 

% 

r  ,,i  .n-hv.\^n\  •iinShfiii 


(5.3-25) 


.s^-hro\k  ol  <i» 

■  n  [^)  ^  sin^effjy'i 

,  ft  I  si^^>^^sins»cosee  de 

I  //^‘A*2*  f  sh*&^  -yji 


(5.3-26) 


L  ■  *  I  s}n£nd  s}nt9C0SSBd& 
VyA*'*^/ j  /VA*  V.  i)/** 


^ '/5ztfw 

An  interesting  observation  from  equations  (5.5-22)  and 
(5.3-24)  is  that  the  induced  velocity  from  the  bound  vortex  is 
an  odd  function  and  that  many  of  the  odd  terms  are  zero.  All 
"h”  terms  which  are  not  multiples  of  the  number  of  blades  are 
zero. 

The  contribution  of  the  free  vortex  system  of  the  nropcller 
to  the  Fourier  coefficient  u^(Z)  is  obtained  from  equations 
(5.3-20)  and  (5.2-3)  as 


Tr 

J(cosn<p)c^  J<P 


=  (cos/ 


iTjKTjl  ^Jlp  +  (2hz-^^  -  Z(^Jccs(f-/fp -/)p2]  1 


^  ir  ^ 
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(5.3-27) 


where 


rf 

^jcQsne^ 

-fT  o 

(5.3-28) 

By  a  similar  analysis  it  can  be  shown  that  the  contribution 
of  the  free  vortex  to  the  coefficient  v^(z)  is 


. 

J  lb  ^  (fAz  f  «  ^ UylSi)^-  2(^}Cos(b I 


re 


(5.3-29) 


where 


(5.3-30) 


r  sifin$l 


^[Zhz-f<xfaHM^)]cos(5-hx)  -('^}fa.t,/iiSin{0Hx)  \ 


\Jo(. 


Je 


From  equations  (5.3-20),  (5.3i-,22)  and  (5,3-27)  and  equations 
(5.3-21),  (5.3-24)  and  (5.3-29),  the  Fourier  coefficients  for 
n^O  are  obtained  as  follows 


'^n(Z)-  = 


0  if  n  /  mb 


if  n  =  mb 

m  =  1,2,3... 

(5.3-31) 


Ill 


and 


if  h=:  m6 


^  O)  if  n 

(5.2-32) 


It  will  be  noted  that  the  free  vortex  system  of  the  propeller 
contributes  to  both  the  odd  and  even  terms  while  the  lifting 
line  vortex  contributes  to  only  the  odd  term. 

The  coefficients  irn^^^  can  easily  be  obtained  but  the 
coefficients  ^ns^^^  must  be  discussed  farther.  These 

last  two  coefficients  are  a  s necial  form  of  the  equations  dis¬ 
cussed  in  Appendix  E,  The  coefficient  jnc^»  equation  (5.3-28), 
is  given  by  equation  (E-18)  where  ^  «  — rrr-y  “  — 2 - and  jns 


is  given  by  equation  (E-9)  X'd-th  the  same  value  of  C  ,  It  should 
be  noted  that  has  a  logarithmic  singularity  at  x^=Xjj,z=0. 

This  singularity  can  be  removed  by  the  method  used  for  equation 
(5.3-14). 

With  these  Fourier  coefficients  u^^(z)  and  ^^^(z)  the  inte¬ 
gral  equation  for  the  strength  of  the  ring  vortex  distribution 
representing  the  duct,  i^((/^,z),  can  be  solved.  If  the  Fourier 
expansion  of  the  ring  vortex  strength  (5.3-5)  is  substituted 
into  equation  (5.3-1),  an  integral  equation  is  obtained  for  the 
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Fourier  coefficients,  gn^z)  and  h^^Cz),  of  the  vortex  ring 
strength. 

j,  -t,  f ]  j 

JJ  1^“ 

o  ~rr 


(5.3-33) 


i  r 

o  ^ 


I  r  o«9 


where 

z=*-a^ 

The  integral  of  each  series  in  the  above  equation  will 
now  be  considered  separately.  First,  using  the  trigonometric 
identitv;  1-cos  (  cP- 'J?»)=2sin^  X 


<iA 


Now  make  the  change  of  variable  =  P  and  use  the  trigono- 

2 

metric  identity 


cos  n«/^'=cos  n^cos2n0  +  s  in 

The  last  term  integrates  out  over  the 

_  in  ^  jr 

2=^=2 


n  s  1  n  2n0 

integration  range  of 


\sr^(z')[Zh(i-m 


f 

^  ^Ccosn<Pcoi8n§  -i-  smhcfs/MBh§)coi2&  j 

\_J  Ofh^(l-Z')^  + 


-f-l 
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CO  ^ 

-  ‘f/^5h(pLf,(zO[8f,(Z-z‘}] 
»|  =  ^  6 


cosenBcoszG  da 


The  integrand  of  the  integral  is  periodic  with  respect 
to  0  so  it  makes  no  difference  whether  the  integration  is 
carried  out  over  the  range  -j-^=  5  =  -j-^or  ^  j 

and  since  the  integrand  is  an  even  function,  the  integral  over 
the  range  from  0  to  ^  is  twi'-e  that  from  ~^to  j  .  The  integral 
in  the  brackets  can  be  reduced  to  a  function  tabulated  by 
Riegel by  making  a  change  in  variable  0=0+  ~  and  using  the 
trigonometric  identity 

2cos2n9cos26  =  cos2(n+l)0  +  cos2(n-l)6 


QO 


CosHcf  I  ^n(i')[sf>(z-t)] 

i  I  » 


% 


[cos2(nH)a  ^cosa(t',  -O^lds 


J  [¥h^(Z-2')^  fV  5/^(5 
0 


Jz' 


hTi'  y  \J~  I 

"  '  O  o 

00  / 


n=o  ^ 
where 


k2  r: 


1 


and 


h^Cz-^z*  )2+i 


3 


(5.3-35) 


o 


♦  37 

Rxegels  has  given  the  expansion  of  this  function  Gt^dc) 
as  k  -*  1,  which  is  when  z'-*  2,  as 


GnCk)  = 


(5.3-36) 

thus  G^Ck)  has  a  singularity  as  k  1.  This  singularity  is 
removed  by  multiplying  by  (i  -  z*)^. 


From  this  limit  it  can  be  seen  thrt  the  function  within 
the  braces  of  equation  (5,3-34)  is  not  singular  as  z*  -*■  i, 
but  the  complete  integrand  has  a  singularity  at  z*  *  2.  Values 
of  GQ(k)  and  Gi(k)  can  easily  be  obtalne.d  in  the  form  of  complete 
elliptic  Integrals.  These  two  values  are 

ft  /i)  - 

(5.3-37) 

The  asymptotic  expulsion  of  Gj^(k)  for  k  •«■«  1  ia  given  by 
Riegels^^  as 


(5.3-38) 


The  integral  involving  the  scri-es  fhn(z')  sin  n  in 


115 


equation  (5.3-33)  reduces  in  a  similar  manner  to  the  follow¬ 
ing  equation 

t  "tf  « 

■  r  C0i{cp-<fi')]L,t^^(Z'}5in>^Cp'  j  ^  ^  f 

j  j[j^h-{z-z'}‘  +  £[/ 

O  -IT  ' 


/ 


K=o  ^  (5.3-39) 

The  other  terms  occurrin  '  in  equation  (5.3-33)  are  some¬ 
what  more  complicated  since  the  helical  vortex  sheet  shed 
from  the  duct  gives  rise  to  both  odd  and  even  terms.  This  is 
shown  in  the  following  reduction  of  the  integral  involving 
gj^(z’).  In  this  reduction  the  order  of  integration  is  inter¬ 
changed,  the  change  of  variable  made,  and  the  order 

of  integration  changed  back  again, 

f.k LL.fi  j 

o  -'if  o  ' 

I  ft  oo 

L  f  j\  ^c<  (0 1-;;’ J 

h-Iq  ~it  o 

I  ^  I 


where 
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•r 


-ff"  o 


and 


ir 


,(S; 


(z  -ij  =n[zh(z  -Z‘)]\s  ihAB  I 
-ir 


'[eh(z~Z)-l-0(ttt»/y]£es(S-fct)  -  titn/tjsin(0-Ks<) 

.  jakd-z) +<xt^»/^J^  +  f  s}zi^\(e  +-;]’'2  , 


(5,3-41) 
^oclc/fi 


(5,3-42) 


The  integrals  in  in^^^  flnd  in^®\  equations  (5,3-41)  and 
(5,3-42)  respectively,  have  Lesn  reduced  to  a  simnler  form  in 
Appendix  E  by  evaluating  the  infinite  integral.  Using  Appendix 
E  equations  (5,3-41)  and  (5,3-42)  become 


il^\z-Z')  =  •^[zkCz-lO]!*^  (5.3-43) 

i!,^\i-Z‘)  =  ^[2  h(z  (5.3-44) 

where  and  are  given  by  equations  (E-9)  and  (E-18)  re¬ 
spectively  with  x-=l. 

The  integral  from  the  vortex  sheet  involving  the  Fourier 
coefficient  hj^(z’)  is  reduced  in  similar  manner  and  is 


I  ir  t)£> 
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If  the  Fourier  series  for  U(c^,z)  is  substituted  into 
equations  (5.3~53)  along  with  equations  (5^3-34),  (5.3-39) 
(5.3-40)  and  (5.3-45),  the  equation  is  written  completely  in 
terms  of  a  Fourier  series.  Since  a  Fourier  series  is  unique 
and  linear,  the  coefficients  can  be  equated.  The  result  is 
two  linear  singular  integral  equations  for  the  Fourier  coef¬ 
ficients  of  the  ring  vortex  strength,  g^^z)  and  hj^(z),  in 

terms  of  the  knoxm  coefficients  u„(2)  and  v  (z),  equations 

n  n 


(5.3-31)  and  (5.3-32). 


vAz)  -2|z-Z' 


elz' 


\  (5.3-46) 


where 


n^(i-z')=4\cka-2‘)l%„il!l  f  <5.3-47) 

From  the  equations  for  Uxi(z)  and  Vjj(z)  it  is  known  that 
u^(z)  =  vn(z)  =0  if  n^b,  (m=l,2....).  Since  equation  (5.3-44) 
must  hold  for  all  values  of  Z,.  it  must  be  concluded  that 

gj^(S)  =  hj^(z)=0  n-^.b,  (!n=l,2, . . . . ) 

because  the  coefficients  K^(Z-z')  and  i^®^(Z-z')  obviously 
are  not  zero  for  nj^b.  From  a  practical  point  of  view  this 
fact  greatly  reduces  the  number  coefficients  which  must  be 
calculated.  For  instance  for  a  4-bladed  propeller  in  a  duct 
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only  g^(z),  h^Cz),  gg(z),  etc.  exists,  (It  should  be 
noted  that  exists  for  all  number  of  blades  and  that 

ho(z)  is  identically  zero).  From  the  decrease  in  the  order 
of  magnitude  of  Fourier  coefficients'®  with  increasing  n, 

It  IS  not  difficult  to  conclude  that  the  series  for  the  circu 
lation  distribution  converges  very  rapidly  since  so  many  of 
the  terms  are  identically  zero. 
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V.4  Reduction  of  the  Integral  Equation  for  the  Duct 
Ring  Vortex  Strength. 

The  method  of  solution  of  the  Integral  eq’.:ation8  (5.3-46) 

will  now  be  considered.  First,  however,  the  solution  when 

n  «  0  will  be  discussed.  For  this  case  h^(2)  »  Vq(z)  « 

(e) 

1'  '(z  -  z')  -  0  and  the  system  of  equations  reduces  to  a.’ 
o 

singular  Integral  equation  for  g^(z/). 


(5.4-1) 


where 


Tliie  function  Uq(2)  is  given  by  equation  (5.3-16).  This 
equation  (5.4-1)  is  identical  to  the  one  solved  in  section 
II. 4  with  the  additional  terms  from  the  Induced  radial  veloci¬ 
ties  from  the  hub  and  propeller.  The  solution  method  is  the 
same  except  that  additional  terms  now  occur  in  equations 
(2.6-2)  and  (2.6-6)  for  the  function  f(0).  Following  the 
same  procedure  as  in  section  II. 4,  the  functions  Hj^(^  ^nd  Hp(2)^ 

which  occur  in  the  equation  for  Uq(z) ,  are  expanded  in  a 
Fourier  cosine  series  in  0',  for  0  =  0  =  TT ,  where 
i  =■-  z  -a^.  -|l/2)(l  +  cos  0) ,  then 


(5.4-2) 
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where  ^ 

o 

ir 

fm~f Jl'Hh(e)  +  Hp(e‘)]co5Mi  e'Je' 

0 


(5.4-3) 


and  Hp(d)  and  %(6)  are  obtained  from  equation  (5.3-16)  with 

the  change  of  variable  mentioned  previously.  Substituting 

equation  (5.4-2)  Into  equation  (2.4-14),  the  new  function 

f  (9)  is  obtained 
P 

^ (e)  -  Ue)  +  0^^ sm (5 ^ 4_4) 

vdiere  f(Q),  Is  given  by  equation  (2.4-26). 

This  function  fp(9)  Is  used  In  equations  (2.6-2)  and 
(2.6-6)  In  place  of  f(0).  This  equation  (5.4-4)  Is  for  the 
general  case,  l.e.  when  a  singularity  occurs  In  the  circulation 
distribution  at  the  leading  edge.  It  also  seems  reasonable  to 
describe  an  Ideal  angle  of  attack  of  the  duct  section  when 
the  propeller  Is  In  the  duct.  Obviously,  each  section  of  a 
symmetrical  duct  cannot  operate  at  an  Ideal  angle  of  attack 
in  the  presence  of  the  propeller  but  an  Ideal  angle  of  attack 
can  be  defined  In  presence  of  the  average  velocity.  Making 
the  definition  of  Ideal  angle  of  attack  with  the  propeller  in 
the  duct  as  the  angle  of  attack  In  which  the  singularity  In 
8q(.z)  does  not  occvir  at  the  leading  edge  of  the  duct,  the 
following  is  written  for  the  function  fp^^j  (9) . 


kjtol  - 


oc 


(5.4-5) 
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where  is  obtained  from  equation  (2.5-3). 

To  obtain  the  ideal  circulation  distribution  [gg(z)]j^^ 
in  presence  of  the  average  hub  and  propeller  Induced 
velocities,  the  foregoing  equation  is  used  for  f(G)  in 
equations  (2.6-2)  and  (2.6-6)  and  the  calculations  carried  out 
as  described  in  section  II. 6  for  the  ideal  case. 

By  the  procedure  just  described  the  Fourier  coefficient 
gjj  or  g*,  where  g*(6)  “  8^(6)  can  obtained  by 

the  methods  already  deacribeu  in  sections  II. 4,  II. 5  and 
II. 6.  The  Fourier  coefficients  gj^(»)and  hj,(z)for  n  ft  0 
cannot  be  obtained  quite  so  easily  because  they  are  defined 
by  two  linear  singular  integral  equations  and  this  involves 
the  solution  of  a  system  of  singular  integral  equations. 

The  system  of  singular  integral  equations  given  by 
equation  (5.3-46)  can  be  reduced  to  a  system  of  non-singular 
Fredholm  equations  of  the  second  kind  which  can  be  evaltiated 
by  known  methods  .  The  method  of  reduction  is  the  same  as 
given  previously  in  that  the  singular  integral  equation  is 
reduced  to  the  airfoil  equation  and  then  inverted.  Starting 
with  the  first  equation  of  (5.3-46)  the  term  Mjj(z-2)  - 
is  added  and  subtracted  from  the  integrand.  It  can  easily 
be  shown  from  equations  (5.3-43)  and  (5.3-47)  that 
Mn(0)  *  2,  and  in^'^O)  -  0. 


/  / 


This  equation  is  now  inverted  as  In  section  II. 4.  The 
Kutta  condition,  equation(2.2-21),  is  satisfied,  by  making 
gjjCz)  zero  at  the  trailing  edge  of  the  duct.  Inverting  and 
interchanging  the  order  of  integration  the  previous  equation 
can  be  written  as  a  Fredholm  equation  of  the  second  kind. 


I  I 

_  m f  A, /HZ  [  dz"  iziJz' 

O  o 

(5.4-7) 

A  singularity  occurs  at  z  ■  i  which  is  the  leading  edge 
of  the  duct.  As  discussed  previously,  the  singularity  is 
removed  by  Introducing  a  new  variable  g*(z)  “Vl  "  ^  ^(2) 
also  h*(z)  »  Vl  ”  2  If  the  change  of  variable 

z  -  (l/2)(l  +  cos  0),  z'  -  {l/2)(l  +  cos  0'),  etc.,  is  made, 
the  preceding  equation  becomes 


g;;caj  =  sm^e  =  -f-  /  f  Je 

J  (5.4-8) 


where 


(5.4-9) 
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cos 


V 

1  e'cosi  ~  [g  -  (cosd"- 

“  *'  j  C<:osa"-Ci>Jff)l(coso"  -Ct^Jff'J 


(5.4-10) 


CCS  j6‘ cole 


altl  -CoiO")  L^(C0S6‘‘-  CcSO'jl  i  // 
jCccsfi’^-Ccsfi)  (rcsO"'~  CeseO 


(5.A-11) 

The  last  three  integrals  are  evaluated  by  expanding 
the  Integrands  in  a  Fourier  cosine  series  in  9'  or  0"  and 
using  the  Cauchy  principal  value  integrals  given  by  equation 
(2.4-25),  then 


f„f(9;  s  55:  -dno  Cos^B-t  sih  (5.4-12) 

L  ms( 


where 


0«0  =  \^»(&*)de' 


(5.4-13) 


(&')  s'm  m  S>' dd' 


(5.4-14) 


f«o 

-  br^o  (.^'kcs^  9  +  (5ii^ie^b,,„(e')s'//i  ^0 

(MS! 


where 


1,  ('a')  s  —  2  -  Afw  _ Ids" 

fr  [coSB'' 


(5.4-15) 


(5.4-16) 


12A 


[oo 

•"S/  tct. 


(5.4-18) 


vrtiere 


ir 

Ar.c(S')  ^  K  (<^ose"-CfiSff'Jd&'> 


dnmM  »  ^  jl'^(coi&“ -cosfi')coin^e‘de' 


(5.4-19) 


(5.4-20) 


The  function  is  given  by  equation  (E-18)  with  5c  =  1. 
It  has  a  logarithmic  singularity  at  9"  >■  9'  which  can  easily 
be  removed.  The  integrand  of  the  integrals  for  the  Fourier 
coefficients  b^j^  and  eqiwtion  (5.4-16)  and  (5.4-17), 
are  of  indeterminate  form.  Referring  to  equation  (5.3-47), 
it  can  easily  be  shown  that 


///M  ^  ~ 

g'Uff'  eosa" -coi9' 


The  Integral  is  given  by  equation  (E-9) ,  Appendix 
E,  with  ic  >  1  and  thv?  proper  change  of  variables. 

The  second  Integral  equation  of  the  system,  equation 
(5.3-20),  can  be  reduced  to  a  Fredholm  equation  of  the 
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second  kind  for  the  coefficient  h*(0)  in  the  same  manner. 

n 

Repeating  equation  (5.4-7)  the  coupled  Fredholm  eqtiations 
of  the  second  kind  are 

ir 


ir  ir 

h*l,(e)  =i(^}  +  jHn(a,a1^n(e')de' 


(5.4-21) 


whe"e 


r  Qg 

S//»A»*S  (5,4-22) 

L  ws/ 


tr 

&no  (5.A-23) 
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IT 

Anm  -  (5.4-24) 

o 

also 

^n(e)  f'n(B) 

The  system  of  equations  (5.4-21)  in  the  inter'/al  (0,-rr) 
can  be  reduced  to  a  single  equation  in  the  Interval  (0,2 tr) 
as  discussed  in  Reference  [29].  This  single  equation  is  not 
necessarily  more  convenient  to  use  than  applying  the  various 
solution  methods  directly  to  the  system  of  equations  ^.j.,4-21). 


To  obtain  a  numerical  solution  to  the  system  of  aquations, 
they  are  rewritten  in  the  following  form 
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=;  j^-  At,ct>sj  fi  f  A,  si*i  e -f-  •  -  '  ‘-i-  Am 

(5.4-25) 

hUoi  -  deceive -h  e,s/»e+-  •  •  •  • 

(5.4-26) 

where 


W 


(5.4-30) 


The  Fourier  coefficient  is  given  by  equation  (5.4-16) 
or  (5.4-17)  and  since  they  are  functions  of  geometry  only,  they 
can  be  tabulv'^ted.  Each  of  the  equations  (5.4-25)  and  (5.4-26) 
can  be  solved  by  the  method  given  in  section  II. 6.  The  problem 


now  exists  that  Fj  and  F2  are  xmknown,  consequently,  it  is 
necessary  to  use  the  method  of  successive  stibstitution. 
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Equation  5.4-25  is  first  solved  assuming  Fj^(6)  ■ 

then  the  resulting  g$(d)  used  to  solve  equation  (5.4-26). 

This  process  Is  repeated  using  each  time  new  values  of  Fj (9) 
and  F2(9)  until  satisfactory  convergence  Is  obtained.  All 
the  Fredholm  theorems  apply  to  the  system  of  equations  so, 

In  general,  convergence  Is  assured. 

Once  the  Fourier  coefficients  g*  and  h*  have  been 
determined,  the  circulation  distribution  of  the  ring  vortex 
strength  y  can  be  calculated  from  equation  (5.3-3)  and  the 
Induced  velocities  at  the  propeller  and  hub  can  be  determined. 
In  general,  It  will  be  found  advantageous  to  make  use  of  the 
Fourier  series  expansion  of  the  ring  vortex  strength. 


120 


(5,5-L) 

where 


1  etf 

CT,=jrJ 

(5.5-2) 

1  sir 

0  0 

(5.5-3) 

1  ftr 

(5.5-A) 

The.  radial  velocity  induced  by  the  hub  et  the  duct  is 
given  bv  equation  (5.2-9).  If  this  equation  is  substituted 
into  equation  (5.5-2),  the  contribution  to  the  duct  thxTJst 
by  the  hub  is  given.  This  equation  con  be  considerably  simpli¬ 
fied  i  C  the  Fourier  expansion  of  }^  is  used,  then 

I 

0 

Since  go(z)  may  be  singular  at  the  leading  edge,  the 
change  of  variable  z=  j'd+cos©)  is  introduced  and  then  the 


L29 


equation  can  be  written  in  terr..s  of  the  coefficient  g*  (6) 
which  was  discussed  in  the  last  section. 


U 

C =  a  j  B)f,  cos^bJb 


(5.5-5) 


where 


g*(e)=sin  j  e  go(e) 

The  radial  velocity  induced  on  the  duct  by  the  propeller 
is  given  by  equation  (5.2-8),  If  the  Fourier  series  for  both 
^  and  the  radial  velocity  induced  by  the  propeller  is  used, 
equation  (5.5-3)  is  simplified  to 

Cr2=  -irr  -  2? /[&"  cos^edB 


5  Us/ 


(5.5-6) 


The  function  Uq(9)  is  given  by  equation  (5.3-15)  with 
the  change  in  variable,  z=  j (1  +  co80)  +  a^.  The  functions 
g*(0),  h*(9),u^(9)  and  v^(8)  a^re  the  same  as  discussed  in 
the  previous  section  with  the  necessary  change  in  variable. 

The  radial  velocity  induced  on  the  duct  by  the  duct 
trailing  vortex  system  is  given  by  equation  (5.1-5).  After 
introd'icihg  the  Fourier  expansions  for  and  making  use  of 
the  evaluation  of  the  infinite  integral  in  Appendix  E,  the 
contribution  to  the  nozzle  thrust  by  the  trailing  vortex  system 
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£oi-lrS'cl9*^cos^6'jd0 

(5.5-7) 


The  functions  I*  and  I‘  »re  given  by  equations  (E-9)  and 
X  ^ 

(E-11)  with  the  variable  » j^^(cos0-cos9*)»  /3, -/Hj  and  5c:=l. 

Once  the  contribution  to  the  duct  thrust  by  the  various 
components  have  been  obtained  the  total  duct  thrust  follows 
easily.  It  should  be  remembered  that  only  terms  which  are 
harmonic  with  the  number  of  blades  appear  in  the  series  in 
these  equations,  consequently  a  majority  of  the  terms  are  zero. 


V.6.  Velocities  Induced  by  the  Duct  and  Its  Trailing 
Vortex  System  at  the  Propeller  and  Hub, 
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Only  the  axial  anJ  tangential  velocities  induced  at  a 
propeller  lifting  lin«*  are  desired.  Since  the  flow  is 
symmetrical i,  it  is  sufficient  to  consider  only  the  lifting 
line  at  <jp  »  0.  The  velocity  induced  at  the  propeller  by 
the  duct  can  be  considered  as  that  due  to  the  ring  vortices > 
ring  sources  and  duct  traili  2  vortex  system. 


-  ^fx,oL  4.  i^(x,OjCL+ 

Vs  ^  Vs  f  Vs  ’Vs 


i'f 


(5.6-1) 


and 


^(yAo) 

Vs 


=  )^(ysO,o)y  +  Y±lx,o,o)i 


Vs 


r. 


(5.6-2) 

The  ring  source  system  does  not  induce  a  tangential 
velocity  since  its  strength  is  independent  of  angle.  The 
axial  velocity  induced  by  the  source  ring  is  given  by 
equation  (B-12) .  Introducing  the  source  strength  in  teirms 
of  the  thicknes.s  slope  of  the  duct,  equation  (2.3-4),  the 
axial  velocity  induced  by  the  source  ring  is  obtained  as 


(U 


dz' 


(5.6-3) 


where 


€ 


JiM. 


(5.6-4) 


The  integrand  of  this  integral  is  singular  at  the  point 
X  »  awd  z'  -  -a^.  The  axial  position  z*  -  -a^.  is  at  the 
lifting  line  and  x  -  x^  »  1  only  occurs  at  the  blade  tip.  When 
the  propeller  is  not  within  the  duct  then  the  singularity 
does  not  exist,  i.e.  z*  -a^  does  noc  lie  between  0  and  I. 
When  the  singularity  exist  this  equation  must  be  treated 
as  a  Cauchy  principal  value  lntei;ral.  Considering  the  case 
f-  '  X  ■  x^  then  the  precedl  j  equation  can  be  expressed  as 


I 

^(h oL  =  ^(1  -  Jz' 

^  5  ^  /  *(AfhZ') 


(5.6-5) 


where 


+  f 


If  the  change  of  variable  is  made,  z'  *  (l/2)(l  fi  cos  0)  then 


where 


^0,o)a  =-^(i  -  — :  Sinede 

Vi  'f  (cose  -  COSQi) 


(5.6-6) 


COS0t  -  "(SA* 

®  ”  k^Ccose'  -cos0^)*  +  'i- 


The  term,  s‘(9)k2  E(lc2)  sin9,  is  expanded  in  a  Fourier 
cosine  series  in  9  o7er  the  range  0  to  IT .  This  involves  no 
additional  assxsnptlons  on  s’ (9)  than  in  Section  11.4.  If  this 
is  done,  and  the  series  substituted  into  the  above  equation 
and  the  integration  carried  out,  the  induced  velc-city  is 
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finally  given  as 


where 


Vs  ^  Sine<  Z." 


It 

5^  -  3  '(B)4e  £(-(^J  sfM  ecosndde 


(5.6-7) 


The  axial  velocity  induced  by  the  ring  yorcex  at  the 
propeller  lifting  line  is  given  by  equation  (A-18)  as 

^(xoo)  --  H  i dcp'dz' 


(A-18) 


If  the  Fourier  series  for  y  is  introduced,  then  this 
equation  becomes 

<  OO  ^ 

)^(X0  0)  -  ^  f  fVq  (A 


(5.6-8) 


The  function  k  is  given  by  equation  (5.6-4)  and  the 
2 

Fourier  coefficient  hjj(2')  of  the  ring  vortex  strength  inte¬ 
grates  out.  By  algebraic  manipulation  this  equation  is 
reduced  to  a  function  of  elliptic  and  Rlegel  functions. 


f(W0l  = 


I  ^ 

-'VfrJ  (i)»|4|"  1 J  0  -  J/ 

i 

_ q,(z')\m(  )  -  EIS]- 
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-  ^  *0,., «))]</*' 


(5.6-9) 


The  G^(k2)  «re  fimctions  tabulated  by  Riegels^'  and 
are  given  by  equation  (5.3-35).  There  are  a  number  of 
singularities  appearing  In  the  Integrands  when  the  axial 


coordinate  z'  and  x 


1,  but  all  are  removable. 


As  In  the  case  of  eqtiatlon  i  .6-3)  this  difficulty  only 
occurs  when  the  propeller  has  the  same  diameter  as  the  duct 
(x^-1),  the  axial  coordinate  Is  at  the  lifting  line  (z'^-af.), 
and  the  velocity  is  desired  at  the  point  (x->l) .  In  addition 
to  the  singularity  of  the  Integrands  at  this  point,  It  must 
be  remembered  that  8^(2*)  oay  also  have  a  square  root  singu¬ 
larity  at  the  leading  edge  of  the  duct.  For  this  reason 
the  change  of  varxablc  z'*  ^  (l+co^0*)  Is  introduced  into  the 

preceding  equation  and  the  function  g*(9)  obtained, 

•fr 

<1  ,  (5.6-10) 


where 


-(2av+l) 
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(5. 6-11) 


If  ■  X  ■  Ijthen  equation  (5.6-10)  reduces  to  the  follow¬ 
ing. 

II 

o 

-im  (5.6-12) 

^ere 


)^(0)  is  the  ring  vortex  strength  at  z  -  0 


h^(cos&‘  -cos f  ¥ 


The  ring  vortex  strength  /  (0)  arises  from  the  proper¬ 
ties  of  vortex  sheets^^  (see  Appendix  A).  In  this  equation 
the  integrand  of  the  first  integral  has  a  logarithni'  oingu- 
larity  at  k2*l  arising  from  the  elliptic  integral  K.  The 
integrand  of  the  second  integral  also  has  a  logarithmic 
singularity  arising  from  the  functions  Gn  at  k2“l.  To  show 
this,  the  expansion  of  the  functions  Gf)  near  k2*l,  equation 


(5.3-36),  is  substituted  into  equation  (5.6-12).  The  singu¬ 
larity  of  the  type  (  ^im^  cancels  out  and  a  logarithmic 

singularity  is  left.  The  logarithmic  singularity  causes  no 


difficulty  since  it  is  a  removable  singularity  and  Integrates 
out.  For  convenience  in  numerical  calculations  the  change  in 
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variable 

COS0'  -  COS0J.  ■  C08% 

may  be  Introduced  into  equation  (5,6-12)  to  remove  the 
singularities. 

The  tangential  velocity  induced  by  the  ring  vortex 
distribution  at  the  propeller  lifting  line  is  given  by 
equation  (A- 19).  If  the  Fourier  series  for  the  circulation 
distribution  is  Introduced  1  to  this  equation,  it  is  obvious 
that  the  even  series  Integrates  out  and  only  the  odd  series 
is  left,  consequently 


13V 


The  Integrand  of  this  integral  does  not  have  any  singu¬ 
larities  and  is  zero  when  k2  ■  1. 

The  axial  velocity  induced  by  the  duct  trailing  vortex 
sheet  at  the  propeller  lifting  line  is  given  by  equation 
(5.1-2)  with  «0  and  z  «  0.  Again  introducing  the  Fourier 
series  expansion  for  the  circulation  distribution  this 
equation  becomes 


where 


(5.6-15) 

The  integrand  of  the  infinite  integral  is  singular  at  the 

point  X  ■!  xj,  o<»  cjp*  -  2h(a>  *g*)  and  is  only  singular  when 

tan/?/ 

the  propeller  is  the  same  diameter  as  the  duct.  The  infinite 
Integrals  occurring  in  this  Integral  are  almost  the  same  as 
given  by  the  propeller  induction  factor,  equation  (4.3-7). 

If  X  yt  x^,  the  order  of  integration  can  be  interchanged  in 
this  equation  and  the  infinite  Integral  can  be  reduced  as  in 
Appendix  E. 
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cosQ^  -  -(2at+l) 

If  the  behavlrr  of  the  double  Integral  in  equation 
0.6-14)  is  considered  at  then  it  can  be  shown  that 

ifli  and  i^2  are  valid  even  for  since 


Cl  -  =  z  s/«y 
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The  tangential  velocity  Induced  by  the  duct  trailing 
vortex  sheet  at  the  propell  t  lifting  line  is  given  by  equation 
(5.1-4).  Introducing  the  Fourier  series  expansion  for  the 
circulation  distribution  this  equation  becomes 
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(5,6-21) 


cos0{.  ■  -(2at  +  1) 

These  last  two  eqxiatlons  have  no  singularities  In  the 
Integrands  of  the  Integrals.  There  are  some  Indeterminate 
forms  of  0/0,  which  are  evaltiated  as  follows 


If  f  =  0 
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with  the  preceding  equations  the  axial  and  tangential 
velocities  Induced  by  the  duct  can  be  determined  at  the 
propeller  lifting  line  by  equations  (5.6-1)  and  (5.6-2). 

These  equations  allow  the  propeller  to  have  any  axial  posi¬ 
tion  in  relation  to  the  duct  and,  in  addition,  the  diameter 
of  either  the  duct  or  propeller  is  arbitrary. 

In  addition  to  the  velocity  induced  by  the  duct  at  the 
propeller,  the  radial  velocities  induced  at  the  hub  are  needed. 
Because  of  the  way  the  hub  problem  was  treated  only  the  average 
radial  velocity  is  desired.  Ihe  average  radial  velocity 
induced  by  the  propeller  at  the  hub  was  essentially  derived 
earlier  in  Section  V.3  and  is 


I 


^  '  (5.6-22) 

'X 

The  function  iy^  (  ^  ,z  )  is  given  by  equation  (5.3-12) 
with  x^  replaced  by  The  radial  velocity  induced  at  the 

hub  by  the  ring  source  distribution  is  derived  in  Appendix  B 
as  equation  (B-13)  and  will  not  be  repeated  here.  The  radial 
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velocity  induced  at  the  hub  by  the  ring  vortex  distribiition 
can  be  derived  from  equation  (A- 14) .  If  the  Fourier  serlee 
for  the  ring  vortex  strength  Is  substituted  into  this  equation, 
the  average  ra41al  velocity  induced  at  the  hub  Is  found  to  be 


it 

1^,  fsaVjcoi (cos 0  -COS 

J  \ 


(5.6-23) 


The  modulus  k  is  always  different  than  one  so  the 
Integral  is  always  regular.  The  average  radial  velocity 
Induced  by  the  trailing  vortex  system  Is  obtained  from 
equation  (5. 1^3).  If  the  Fourier  expansion  for  the  circu¬ 
lation  distribution  is  cubstltutcd  Into  this  equation  and 
then  the  equation  integrated  with  respect  to  <f  from  0  to 
2rr,  it  is  easy  to  show  that  the  average  radial  velocity  at 
the  hub  is  zero,  l.e. 


Many  of  the  coefficients  derived  in  this  section  are 


dependent  on  geometry  only  so  can  be  tabulated.  Specifically 
the  coefficients  I^.^,  and  should  be  mentioned. 
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VI,  Conclusions 

A  theory  for  the  ducted  propeller  is  developed  which 
can  be  used  for  design  purposes.  The  method  is  based  on 
the  assumption  of  an  inviscid  fluid  and  that  the  propeller 
can  be  represented  by  a  lifting  line.  These,  among  others, 
arc  the  assumpcions  normally  made  in  propeller  theory^^  and, 
consequently,  corrections  must  be  introduced  to  allow  for 
the  difficulties  in  these  essumptions.  These  corrections  are 
not  considered  here. 

Assuming  that  lifting  sux-face  and  viscous  corrections 
can  be  made,  the  adequacy  of  the  theory  mvist  still  be  based 
on  experimental  results.  A  linearized  theory  similar  to 
that  used  for  treating  the  annular  airfoil  has  been  very  use¬ 
ful  in  two-dimensional  airfoil  theory  and,  likewise,  the 
treatment  of  the  propeller  by  lifting  line  theory  has  been 
effective.  It  is  not  self-evident,  however,  that  in  combination 
t!ie  resulting  theory  of  the  ducted  oTOoeller  will  be  satis¬ 
factory.  It  is  presumed  that  it  will  be. 

An  attempt  has  been  made  not  to  restz'ict  the  problem 
more  than  the  basic  assumptions  which  are  given  in  the  intro¬ 
duction,  consequently,  somewhat  cumbersome  equations  are 
obtained.  In  general,  these  equations  have  been  reduced  to 
coefficient  fonn  which  are  dependent  only  on  geometry  so  can 
be  tabulated. 

Accepting  the  basic  assumotions  in  Section  1,  many 
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important  observations  can  be  made  about  the  anriulax-  airloil 
and  the  ducted  propeller. 

1.  The  linearization  of  the  duct  boundary  conditions 
results  in  the  radial  velocities  induced  by  the  singulari¬ 
ties  in  the  flow  being  equal  the  duct  surface  slope.  In 
addition  the  boundary  conditions  are  only  satisfied  on  a 
cylinder  of  constant  diameter  and  chord  equal  to  the  duct 
chord. 

2.  The  strength  of  the  d  ,t  source  ring  distribution  is 
independent  of  angular  position  and  dependent  only  on  the 
duct  fiiickness  distribution  even  in  the  presence  of  the 
propeller. 

3.  The  strength  of  the  duct  ring  vortex  distribution 
depends  on  both  the  camber  and  thickness  distribution  of 
the  duct  •!».s  well  as  the  radial  velocities  induced  by  the 
oropeller.  If  the  duct  is  at  zero  incidence  and  the  pro- 
oeller  is  not  present,  the  ring  vortex  distribution  is 
indeoendent  of  angular  position  and  no  free  vortices  are 
shed  from  the  duct. 

4.  When  the  duct  is  at  an  angle  of  attack  (no  propeller 
present),  the  duct  vortex  strength  is  a  sum  of  two  terms. 

One  term  is  the  airfoil  at  zero  incidence  and  the  other  is 
dependent  on  only  the  angle  of  attack,  duct  chord -diameter 
ratio,  and  the  angular  positi.on. 

5.  The  problem  of  the  design  of  the  propeller  in  the 
duct  reduces  essentiatly  to  the  nroblem  of  the  propeller  by 
itself  with  the  inclusion  of  the  velocities  induced  by  the 
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duct.  This  means  ihe  indnrrion  fat'stors  which  have  been  calcu¬ 
lated^^  are  applicable  to  this  problem. 

6.  On  expanding  the  radial  velocity  induced  by  the  propeller 
on  the  duct  in  a  Fourier  series,  it  is  found  that  only  terms 
exist  which  are  harmonic  with  the  number  of  blades  except  for 
one  zero-order  term.  This  means  that  the  Fourier  coefficients 
of  duct  vortex  strength  are  harmonic  with  the  number  of  blades 

cpt  £'07*  one  term.  This  zero-order  terra  embodies  the  pro¬ 
peller  average  radial  velocity  and  the  duct  thickness  and  camber. 
All  the  higher  order  terms  of  the  duct  vortex  strength  are 
functions  of  only  the  propeller  induced  velocity. 

7.  The  duct  ring  source  distribution  induces  no  tangential 
velocity  at  the  propeller  nor  does  the  source  distribution 
representing  the  hub. 

8.  The  duct  trailing  vortex  system  induces  no  average  radial 
velocities  at  the  hub. 

9.  The  induced  drag  of  the  duct  is  zero  if  the  duct  is  by 
itself  and  at  zero  incidence.  In  the  presence  of  the  propeller, 
the  induced  drag  is  dependent  on  the  radial  velocities  induced 
by  the  hub,  propeller  and  duct  trailing  vortex  system. 
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Appendix  A 

^og.^ty  Induced  by  a  Vortex  Distribution  on  a  CvIinHpri3 

First  consider  the  stream  function  and  velocity  distri¬ 
bution  induced  by  a  single  vortex  ring  of  diameter,  R^, 
located  at  ^  .  The  figure  shows  such  a  rin,;;  where  both  a 

cylindrical  coordinate  system  (r,f  )  and  a  cartesian 
coordinate  system  are  used. 


The  coordinate  f  is  in  direction  of  tht  axis  of  the  vor 
tex,  the  element  of  the  vortex  filament  is  at  the  point,  P', 
and  the  stream  function  and  velocity  distribution  is  desi-ed 
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at  the  point,  P.  The  Biot-Savort  law  gives  the  velocity 
induced  by  a  vortex  filament.  In  vector  notation  and  for 
an  clement  of  filament  ds,  this  law  is  given  by  the  following 
equation: 


where 


JV  =  - 


HIT  n" 


(A-1) 


is  the  induced  velocity  at  P  from  the  element 
R  is  the  vector  from  point  P*  to  P 

rfs  is  the  incremented  vector  ’^ds  tangent  to  the  v^ortex 
ring  at  p* 

^(<P)  is  the  circulation  at  P  and  dependent  on  * 

If  the  unit  vectors  T,  arc  in  direct  ion  of  the  x,y,^ 
axes,  the  unit  vector  ds  can  be  written  as 

3? 


Js  _  fflichp' 
ds  ds 


t  cos 


Since  R(jd<^^'  =  ds  the  above  equation  becomes 


(A-2) 


(A- 5) 


From  the  figure  it  can  be  seen  that  the  radius  vector  R 
from  P'  to  P  is 


R  =(t'cos<f  -  Cos<f')c 


Bv  vector  multiplication  it  follows  that  Uxd^  is 

3? 
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RX^  ^  i+  f+-  \rJ  -rcos(<f>-<^4 

(A- 5) 

The  tnagnitude  of  R  is.  given  by 
R  =  l^j  sVff/  -  Zr  Rj  co5C<(^-cf')  (A-6) 


Substituting  equations  (A-5)  and  (A-6)  into  equation  (A~l) 
and  integrating  around  the  ring  gives  the  velocity  induced  at 

P  by  a  single  vortex  ring. 
etr 

iv.i  -Ji  [![((- 

\  [rJ  +  (A-7)  / 

For  the  induced  velocity  from  a  vortex  cylinder,  the 

induced  velocity  from  the  single  vortex  ring  is  integrated 

along  the  cylinder  from  the  trailing  edge  (f  =  a^)  to  the 

leading  edge  (f  =  aj^).  Since  each  ring  v/ill  not  necessarily 

have  the  same  strength,  the  vortex  distribution  will  also  be 

a  function  of  f.  It  follows  then  that  for  the  vortex  cylinder 

the  induced  velocity  is 
aj  itr 

°  (A-8) 

As  expressed  by  equation  (2,2-23)  the  axial  coordinate 

(^  )  will  be  assumed  to  be  nondimer sionaliaed  by  the  cylinder 

chord  (a),  the  radial  coordinate  (r)  by  the  propeller  radius 

(Rp),  the  velocity  and  also  the  vortex  distribution  by  the  ship 
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speed  or  in  the  case  of  uniform  flow,  the  frec-stream  velocity. 
The  duct  chord  diameter  ratio  Ch  =  a/2R^)  will  also  be  intro¬ 
duced  and  then  making  the  change  cf  variable  (2”  =  2’ 
and  noting  ttat  in  the  nondimensionalized  form  u^-u^=l,  the 

preceding  equation  can  be  written  as 
I  err 


^  (A-9) 

If  there  is  no  clearance  between  the  propeller  and  the 


duct,  then  =  1  arid  the  above  equation  becomes 
I  Sir 

W-ra,]  -  A  r  -rjcoitf’jZ  +  [Sh(z-Af£')c,fS(p07  -  'U^  \  y^y,. 

Vs  J  [/ 

°  (A- 10^ 

From  the  figure  at  the  beginning  of  this  section,  it  can 

be  seen  that  the  velocity  in  the  axial  direction  (w^)  is  given 

by  the  component  in  the  direction.  The  radial  velocity 

component  (wj.)  and  tangential  velocity  component  (w^)  are 

obtained  using  the  relationship  of  velocities  in  cylindrical 

and  cartesian  coordinates.^®  The  velocity  (w„)  is  the 

component  of  velocity  in  the  1  direction  and  (w^)  is  the 

component  of  velocity  in  the  j  direction,  the  radial  and 

tangential  velocities  are  given  by  the  following  equations 


Wr  =  Wjjcos  (P  +  Wysin  <p 

Wt  =  -WjjSin(^  +  WyCos 

Applying  these  relationships  to  equation  (A-9)  the 
different  velocity  components  are  found  to  be 


(A-11) 


1'56 


I  etr 

_h_f  f _ “  ?  ^(ip'jZjci^yz' 

o  o 


!  efr 

h  f  f  2  -aj-  -Z')coS (Cf- Cp')  /(f,z')J<f>‘Jl‘ 

J  if  (A-13) 

o  o 


/  2-jy^ 

'dtH<JPz\  =--t'-f  -Z')sm(cf>-»f‘)if'(<f>'tr)  Jtf'Jz' 

V,  "  /  -  2rr/  J[,  -  £lfjc^s(cf^‘)]^‘ 

b  o 

The  velocity  induced  on  the  vortex  cylinder  itself  is 
found  by  considering  the  properties  of  vortex  sheets, The 
velocity  across  a  singular  vortex  sheet  has  a  discontinuity 
in  the  tangential  velocity  component  of  magnitude  As  a 
result  the  value  of  the  induced  velocity  depends  on  the  side 
of  the  sheet,'  Since  only  the  axial  and  radial  velocities 
will  be  needed  on  the  ring  itself,  it  follows  then  from 
vortex  theory  that 


'STT 


I  sir 

L_r  r  CcoH<f-cP')  - 1]  m:z'}  Jcp'Jz'  +  i.  yM.) 


o  o 
sr 


(A-.15) 


M 

✓ 


2h(z-af  -r}coi(ce-cpjy((f]z>)Jcp'J/,i 

t  2  -  ZCo!(ce-cf')f^i  ■ 


<=>  o  (A-16) 

The  plus  sign  refers  to  the  outside  of  the  duct  and  the 
minus  sign  to  the  inside-  The  radial  velocity  will  also  be 
needed  on  the  nroneller  hub.  Denoting  the  nroneller  hub  radius 
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radial  velocity  is  given  by 

/  STT 

-  ~^~  f  /}l(^  -z‘)cos{^. r)nY,^o  d  v>' 

(A-17) 

Th«  «lal  and  tanganti<l.  vaUoitv  induced  at  each  nronelXer 
blade  Uiu  alao  be  needed.  For  thia  nurnoae  it  ia  sufficient  tc 
consider  the  propeller  blade  in  its  vertical  position  thus  It  is 
sufficient  to  consider  the  an  '.e  f>  as  aero  and  since  z=o  at  the 
propeller,  tlie  axial  and  tangential  velocities  induced  by  the 
vortex  cylinder  at  the  propeller  are 

I  iir 

^  ^±.f  I Acp'Jz‘ 

(A~18) 

I  srr 

l-AhlAtt^Osincf'm'.zl  dcp'dz' 


If  the  circulation  distribution  s)  is  a  function  of 
the  angular  position,  then  a  vortex  system  is  shed  from  the 
nossle.  The  velocity  induced  by  this  system  must  also  be 
considered  and  is  discussed  in  Appendix  0  and  D. 

If  the  circulation  distribution  is  independent 

Of  the  »,gle  f  ,  tpan  the  velocity  Induced  from  oolnt  to  point 
is  independent  of  the  angular  poaitlon  of  the  noints.  It  is 
then  sufficient  to  consider  that  the  angle f  is  sero.  From 
the  figure  at  the  beginning  of  this  section  it  can  be  seen  that 


YZ 

(A-19) 
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for  <f  =0  the  radial  velocity  at  P  is  given  by  the  ortmnor.or.t 
of  velocity  in  the  i  direction  and  the  tangential  velocity 
component  is  given  by  the.  comoonent  in  the  j  direction. 

From  considerations  of  symmetry  it  is  obvious  that  the  tangpn~ 
tial  component  of  velocity  is  zero.  This  also  follows  from 
the  fact  that  the  coefficient  of  j  is  an  odd  function  and  the 
integral  of  this  term  vanishes.  Utilizing  these  results, 
poijation  (A-9)  reduces  to  the  following  equation. 


I  sir 

JL  fmlJll  fl -£'jCoscy]T-[^C‘»^''>]  1  ]  J[ ^ 

J  0  I  / 


(A-20) 


where 


"  H^(z-a^-z')^  f  (I 

The  integral  in  the  brackets  can  be  reduced  to  the  form 

of  complete  elliptic  integrals  by  introducing  the  change  in 

variable  co8^(P'  -  sinO,  then  d.'f'  ~  ~2d0  and  when  =0,  0=  % 

and  when  0‘  =  2fr  ,  0  =  3  IT  , 

2 

I  %  ^ 


I 


+  Kri)  -  Em 


_ -I)  ^(i) 

HHz-a,-p)^  ^(f.  -/)2 


(A-21) 


where 


I^Ji  "-4>w='<9 


comnlete  elliotio  inteFrnl 
first  kind 


EfiJ  =  Nl-^fs^f^dJe  = 


comnlete  elliptic  inte^^ral 
second  kind 


of  the 


of  the 


ihc  axial  component  of  velocity  is  given  hy  the  coefficient 
of  the  unit  vector  k  and  the  radial  component  by  the  coefficient 
of  the  unit  vector  i. 

As  discussed  previously  the  axial  comnonent  of  the  induced 
velocity  is  discontinuous  across  the  vortex  cylinder  therefore 
on  the  vortex  cylinder  itself,  i.e.  (x  =  a^  -  z  ^  ),  the 

induced  velocities  are 


(A- 22) 


where 


a 


Appendix  B 

Velocity  Induced  by  a  Source  Distribution  on  a  Cylinder 
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In  the  derivation  of  the  velocity  field  from  a  source 
cylinder,  first  the  velocity  induced  by  a  single  source  ring 
of  diameter  located  at  f ^will  he  considered.  As  for  the 
vortex  ring  a  cylindrical  coordinate  system  )  and  a 

cartesian  coordinate  system  «^,y, ^ )  are  considered  as  shown 
in  Figure  8,  The  velocity  inc  ’ced  by  an  element  ^  of  the 
source  ring  of  diameter  R^j  at  the  point  P  is^ 


“  vtr 


(B-l) 


The  strengt'i  of  the  source  at  the  point  (R^,  is 
denoted  by  q(^ ,  ).  The  vector  R  and  its  magnitude  R  rave 
the  same  value  as  given  by  equations  (A-4)  and  (A-6),  there¬ 
fore  the  induced  velocity  for  a  single  ring  is  given  bv 


fir 


^ ^  -ar  RjC os (cp-cp')]^'e 


fV/iJ 


(B-2) 

The  induced  velocity  is  given  by  integrating  single 


vortex  rings  along  the  cylinder  from  the  trailing  edge  (f  =a^) 
to  the  leading  edge  (f  =  a^_  ).  Since  each  ring  will  not 
necessarily  have  the  same  strength,  the  source  distribution 


will  also  be  a  function  of  ^  .  For  a  source  cylinder  it 
follows  then  that  the  induced  velocity  is 
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As  for  the  vortex  ring  the  axial  component  will  be  non- 
dimensionalized  by  the  cylinder  chord  (a)  and  the  radius  by 
the  propeller  diameter  (Ro)»  The  induced  velocitv  nnd  sourrp 
strength  will  be  assumed  to  be  nondimensionalized  by  the  ship 
speed  or  in  the  case  of  uniform  rlovv  by  the  free-stream 
velocity.  The  chord  diameter  ratio  of  the  cylinder  (h  =  a/2Rd) 
will  be  introduced  and  the  change  of  variable  (z"  -  z*  +  ^). 
Since  in  the  nondimensionalized  form  (a^  -  =  1),  the  pre¬ 

ceding  equation  can  be  written  as- 


°  (B-4) 

The  axial  component  of  velocity  is  given  by  the  component 
in  the  5  direction  and  the  radial  and  tangential  velocities 
follow  from  eejuation  (A-11) 


(B-5) 


.jlT 

V,  ^  eir  J  j  t ^ 

^  <’  '  CB-6) 


I  sir 

V,  t  2fr  J  J I  -tTTc^/  -  2(%}cosCf- 

O  O  '■  ' 


(B-7) 
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The  velocities  induced  on  the  source  cylinder  itself  are 

found  by  considering  the  oroperties  of  source  surf aces. In 

this  case  the  tangential  velocities  are  continuous  across  the 

surtace  >>u!‘  the  normal  velocities  are  discontinuous.  Denoting 

the  outside  of  the  cylinder  by  u  plus  sign  and  the  inside  by 

a  minus  sign,  it  follows  then  that  the  radial  velocity  on  the 

•  58 

cylinder  itself  (Ic  =  x^,a^  2  ^  a^)  is  given  by 


_ [I  fUM 

[¥h^(Z-a(-Z‘)^  +  2  -Rcos'C(^~(p')P'^\  ' 

(B-8) 


If  the  source  distribution  is  considered  a  constant  with 


respect  to  <f> ,  then,  since  the  flow  is  symmetrical,  it  is 
sufficient  to  consider  that  the  angle  cP  is  zero  in  equation 
(B-4).  In  this  case  the  velocity  component  in  the  j  direction 


is  an  odd  function  and  integrates  out.  If  the  change  in  vari¬ 
able  cos^^  =  sin0  is  made,  equation  (B-4)  can  be  reduced  to 
the  form  of  complete  elliptic  integrals. 


[!  I  J 

-j.  rWfaij, .. ...liMLiA _ \mi\r 


\nj) 


B[eh(z-a.t-z')]B(4,]  ^ 

Hfj  '0^ 


(5-9) 
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where 


_ viii 


m)- 


''yT^sin^& 


■  uv  ) 


complete  elliptic  integral  of  the 
first  kind 


EW  -j  <36, 


complete  elliptic  integral  of  the 
second  kind 


The  axial  component  of  ^ -locity  is  denoted  by  the  compo¬ 
nent  in  the  k  direction  and  radial  component  is  the  component 
in  the  1  direction.  Again  considering  the  properties  of  source 
surfaces,  the  velocities  induced  on  the  cylinder  itself 
a^  2  z  =  a^  )  are  found  to  be  as  follows 


and 


srr 


I 

/  ''iz-at-z') 


(B-IO) 


:±lrjto,z)^= 

Vs  i  eiij  o 


where 


CB-11) 


/ 

^.‘(z-a^  -zf  t  / 


It  is  shown  in  section  (11,3)  that  within  the  linearized 
theory  the  source  strength  is  independent  of  the  angular  posi¬ 
tion  even  in  the  presence  of  a  propeller.  For  this  reason  the 
velocity  induced  by  the  source  cylinder  at  each  propeller  blade 
is  obtained  from  equation  (B-9)  by  letting  z=0.  Since  only  the 
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axial  velocity  is  needed^ it  is  the  only  one  given. 

/ 

_  Jl  fpilit  L2khfh2'JB(-^J  \j  , 


(B-12) 


whore 


_ 

Similarly,  the  radial  velocity  induced  at  the  hub  by  the 
source  cylinder  is  independei  .  of  angle  and  is  obtained  from 
equation  (3-9)  by  letting  x=^^ 


y'('^h.zk  =  -A  _ '^J U,/ 

(B-13) 


where 


/ _ ‘iSkl _ 
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Appendix  C 

Velocities  Induced  by  the  Free  Vortex  ^vstem  of  a  Vortex  Cylinder 
at  an  Angle  of  Attack 

If  the  ring  vortices  making  up  the  vortex  cylinder  arc  a 
function  of  angular  as  well  as  axial  position,  then  a  system 
of  free  vortices  will  be  shed  from  the  cylinder.  This  is  equi¬ 
valent  to  the  vortex  sheet  v.’hich  is  shed  from  a  three  dimensional 
27 

flat  wing.  Consequently  similar  assumptions  will  be  made 
about  the  shape  of  the  vortex  sheet,  namely,  the  free  vortices 
form  a  cylinder  which  has  a  constant  diameter  equal  to  the 
diameter  of  the  bound  vortex  cylinder  with  generators  which  are 
parallel  to  the  f  -axis,  and  which  extends  from  a  bound  ring 
vortex  to  minus  infinity.  The  density  of  the  free  vortices  is 
1  (.^1  ^  ),  where  y'  )  is  the  density  of  the  bound 

vortices,  ^’/hen  the  annular  airfoil  is  at  angle  of  atteck,  the 
free  vortex  lines  are  straight  lines  parallel  to  the  if  -axis 
while  in  presence  of  a  nropeller  the  vortex  lines  are  helical 
in  shape.  In  this  section  only  the  straight  vortex  lines  will 
be  considered. 

The  same  coordinate  system  and  notation  will  be  used  as  in 


Tne  velocity  induced  by  an  clement  of  the  free  vortex 
follows  from  the  law  of  Biot-Savart,  e<r.atior.  (A-1) 


•frr/ij  1^3  'Udf 

From  the  fisuire  the  vector  d1  is  given  by 


(C-1) 


da  =oTto.Ttd^>«5 

and  the  radius  vector  R  fror?  P*  to  P  is 

R  =  (ycosf~njcoscp')(  (c-3) 

The  magnitude  of  R  is 

R=  |R|  V  V(7-fT  fr‘  fH/  -  ^e/pj(f>r^ 

From  vector  multiplication  of  epuation  C0-2>  and  (0-:,),  RXds 
is  obtained 

RXdl  =  ( Ysm  (p.fij  sm{p‘)i dp  -  (rco$  cf  j 

Substituting  into  equation  (c-l)  the  velocity  induced  at 
P  by  an  element  of  the  free  vortex  is 


~  '^'Ra  %  (.<rccsc(>  - /iy  co^cf'}j  \ 


(C-6) 


^■ince  the  free,  vortex  extends  from  -oo  to  f ,  the  velocity 
induced  by  an  infinitesimal  strip  is  obtained  by  integration. 

~  -RjCoSif')! 

r 

•  f 

-o<? 


The  integral  is  evaluated  and  ecjuation  CC-7)  becomes 


(C-7) 
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filAihtH  /* 


.  [ - ! - 

\tV^+  ft/  -  2i'/?jcosCf-‘f')]  '\j((-f)^  ^  ^f?J'  ••  SffijCos(‘r-<P')  j  I  ^ 

(C-H") 

The  velocity  induced  by  all  the  free  vortex  filaments 
from  2  single  vortex  ring  is  given  by  integrating  this  equation 
with  respect  to  (f '  from  0  to  9'tt'  and  ot  the  free  vortex  cylinder 
by  integrating  the  induced  velocity  of  the  single  rings  which 


are  distributed  along  the  chord  of  the  cylinder. 


_ f[&sitttf-K/S}h<f>')c  ~(rc6S(p Jf-n 


U _ _  4.ti>:^cpuc 

(0-9) 


Nnndimensionalizing  this  equation  as  before,  i.e»  the 
axial  component  with  the  bound  vortex  cylinder  length  (a)  the 
radial  component  with  the  ring  radius  (R^j)  and  the  circulation 
with  the  velocity  (V^)  the  above  equation  is  given  more  simply 


..±^1  fSm(p')i-(^eS<P-C0S<f')  yf  J  SHC-Z') _ ^  ^  )X,d(p'jz 

J  i- 1  -  2^c,s (cf.cp'jj  f 

where  z  =  z-a^,  z  refers  to  the  coordinate  system  at  the  pro¬ 
peller, 

Froir.  the  form  of  this  equation  it  can  be  seen  that  the 
free  vortex  sheet  induces  no  axial  component  of  velocity.  The 


eh(Z-Z') 
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j-adial  and  tangential  components  are  given  by  the  usual  relations 
of  velocities  in  cartesian  and  polar  coordinates,  equation  (A-11), 


I  sr 

=  A  f  I  _ 

I" 


I  dlf 

h^J  =  -  h  ff  _ [ _ 

•T  o  i9 


(C-ll) 


On  bound  vortex  cylinder  ^  =1  and  0  =Z  =1,  the  radial 


velocity,  equation  (C-10),  reduces  to  the  following 
(  err 

M=2^f  [ccoH((f -(?')]  -.  ~  -f  / 

ij  J  j  ^  -/'J* 


/  ^  1  •  ** 
vu-i^y 


and  introducing  the  notation 


equation  (C-12)  becomes 


I  £ir 


(C.13) 
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Appendix  D 

velocity  Induced  by  a  Helical  S’laned  Vortev  T■^n» 

The  velocity  component  induced  at  «  point  in  space  by  a 
vortex  line  is  given  by  the  law  of  Biot-Sevart.  In  vector 
notation  and  for  an  element  of  filament  da  of  constant 
strength  this  law  is  given  by  equation  (A-l)  rewritten 
here  as  (D-1). 


where 


Jv  =__/  RXJs 

ifr 


(D-l) 


dVj^  is  the  Induced  velocity  at  a  point  P 

R  is  the  radius  vector  from  the  vortex  element 

to  the  point  P 

8  is  the  vector  tangent  to  the  vortex  ring  at 
the  element 

y  is  the  circulation  along  the  helix  and  is 
constant 


The  following  figure  shows  such  a  helical  vortex  line 
leaving  from  a  point  P',  which  will  be  taken  as  a  point  on 
the  duct  or  the  propeller  blade,  where  both  a  cylindrical 
coordinate  system  (r,^  ,^)  a  cartesian  coordinate 
system  )  are  used.  The  point  P(x',  y',  ^  ')  Is  a 

general  point  in  apace  and  the  point  P'’(x",  y”,^*')  lays  on 
the  vortex  line. 


Figure  9.  Notation  of  the  hellcoidal  vortex  line 
From  this  figure  it  can  be  seen  that  the  vector  R  is 
given  by 


R  =[rco5cf  -'r^cos(Cfp-h<x)]^  4.  [rs'/^Kf  ~ 


■h-  if  -  - 


CD-2) 
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The  tiiegiiifuds  of  R  is 


R  =  /  ff/  ==  [rcosf-  i'  cc>s( (pf,  ta-jJ'  +  (f  -f'^s  '^ifp  t  r.f .  !>'- roar 


R2  =  1-[^-(-r'cO<n.>,/!cY -^t'l'cCOS(f-f^-o<) 

(D-3) 

The  vinit  rectors  i.  j,  ?  d  k  are  in  the  direction  of  the 

Xj  y,  ^  asas.  The  unit  vector  4“  aiso  follows  froai  the 
fisure  a a 


Jf  -s!n((f^t^)T  i-cosCcpp-hot)^  +  tan/3.  ^ 

~  Sec/2o 


(D-4) 


ana  ds  =  sec  /S  ^doe. 


The  contribution  of  the  xchole  vortex  filament  is  found 
by  integrating  eqviation  (D-1)  with.  2“2Sp'3Ct  to  from  0  ‘  to  -*00 


— oc^  Q 

V  --  L.  I  -  -C  /  MI 

-  W  j  ~  VfTj  R^ 


(D-5) 


Substituting  iu  for  R.  eqvati'-n  ;D-2)  and  ds,  equation 
the  induced  velocity  Irom  a  helical  vortex  line  is 
given  as 


HfP  I  •' 
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fCf-t)  +f ^cascf’ ~\rocos(<f^+o()'}tnn/!^  j J 

'  h' 


ro  -  fe.0i(<f-‘pf,  -ex) 


J 


dai  ,  (o<  Uii/i,<  j 


(D-6) 


Tae  axial  velocity  Induced  by  a  helical  vortex  line  is 
given  by  the  component  of  velocity  in  the  k  direction  and 
the  radial  and  tangential  velocities  follow  by  using 
equations  (A-11) 


Ko  -  ^COSC(f~<ff,  -a) 


p^3 


d‘x 


(D-7) 


*'  att'  I  i  ri7  “ 


v/r 


'll 


(D-8) 


-ffe  f  I  r - (f  -  f'J] f/> -o^}  ^  Cr-n cos(f^ ^ 

~y^l  1  f(3  I 

(D-9) 


The  Integrands  of  the  Integrals  for  the  axial  and  tangen¬ 
tial  velocities  are  singular  ^en  the  point  ?(r,  <f  » f  ) 
lays  on  the  helical  vortex  line  itself.  This  occurs  v4ien 
dm  (cp^  cf  )  m  Cf  ~  and  r  »  Tq.  It  can  be  shown 

*  r  Tq 

that  the  integrand  of  the  integral  for  the  radial  velocity 
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equation  (D-S)  is  not  singular  at  this  point  but  has  a 
jxaap  discontinuity,  except  possibly  at  <=■<  "  0.  HaIs  is 
shovm  by  letting  -  I**)  ■  and  r  »  Tq 

in  this  equation,  using  the  series  expansion  of  the  sine 
and  cosine,  and  taking  the  limit » 


\  [B  f  i'o‘^C<e-4  i-oi) 2  -  s'fi- co>(c( -cp^ 


~  [b[i -co5(<p-(pf, -X)]  -h 


The  limit  of  this  equation  as  oc-*  (  cf  -  cp^)  is 
equivalent  to  considering  ( «  0  and  letting 
9  ■*  0.  Considering  the  series  expansion  of  the  sine  and 


cosine,  i.e. 


Sine  =  e-  ~  +  -  ' ' '  ' 

cose  =  ,  _  ^  ^  _ 

'  2/  ‘f! 


The  above  equation  can  be  written  as  follows 


.  .  .j  +  ©2 
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and  the  limit  of  the  integrand  at  9  0  is 


’t'Otpda  ll'hi 

o-*o 


j  .  7/ 

■  -yr)®‘  +  \ 

_  _ 

Since  the  integrand  is  a  odd  function  of  9,  the  limit 
as  -  9  0  must  be  of  opposite  sign.  Considering  9  as  equal 

-«  )  the  limi.t  of  the  integrand  can  then  be  vritten 


1,'Ht  lj_  C-^ine+  ecos6)tAHAo 
e-^to 


.  =  + 


3H^[l  fD-10) 


Since  the  integrand  of  the  integral  for  the  radial 
velocity  is  not  singular >  it  can  easily  be  ^own  that  this 
infinite  integral  for  the  radial  velocity  is  uniformly 
convergent  with  respect  to  r,,,  r,  ^  ^  *  and  -  fp  with 

the  exceptions  discussed  in  the  following  paragraph. 

At  the  point  r  -  f  *»  »  0,  the 

integrand  is  discontinuous.  If  f  limit 

of  the  integrand  is  taken  as  «  0,  obviously  the  limiting 

value  is  given  by  equation  (D-10).  If,  however, 
oi  •  0  and  the  limit  Is  taken  as  f  or  if  , 

<=< «  0  and  the  limit  is  taken  as  <fp,  the  Integrand  is 
singular  as  .  At  this  point  the  Integral  does  not 
exist  and  consequently  is  not  uniformly  convergent. 
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Appendix  E 

Evaluation  of  the  Infinite  Integrals  in  the  Equation  for,  the 
Velocity  Indu ced  by  the  Trailing  Vortex  System 

In  Section  V.3  infinite  definite  integrals  of  the 
following  two  types  arise 


■(T 

r  _ 

=  Isinne 

-It 


and 


(E-1) 


(E-2) 


where 


In  Appendix  D  the  infinite  integral  was  shown  to  be 
uniformly  convergent  if  z  2*  and  sc  the  order  of  integra¬ 
tion  can  be  interchanged.  If  z  =  z',  it  can  be  shown  that 
the  asymptotic  behavior  of  integral  I,  ih  finite  while  I2  is 

;  )  d) 

infinite  as  lim  In  xL  The  integrals  j__  ,  equation 

(5.3-28),  j^g,' equation  (5.3-30),  equation  (5.3-41), 

and  in^\  equation  (5.3-42)}  are  special  forms  of  one  or  the 
other  of  the  above  integrals. 

For  simplification  the  following  notation  will  be 


introduced; 
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2  “  (z  -  z‘)3i 

*  V 


/j  *  X  tan 


A 

z 


tan/?- 


-tanV 


<E-3) 


t(l  +  x)2  -  4  X 


^  2 
4  X  cos  5#j 


n.e  win  b,  dlscw^.d  first.  IntdrcU.„>,e 

the  or<ter  of  Integration  end  wake  the  ch«.ge  of  variable 
1  -(l/2|(ff +o<)  In  equation  (B-1).  then 

ckS>F*-‘< 

I,-2f  [si»att0cosfic(  -  2 ^  \ 

d  4  f.01  f ^  d  i-  - Hcos 

It  makes  no  difference  whether  the  above  Integration 

is  carried  out  over  the  r-rirs&o  ^  ^  ^ 

ir^  ^  4,-  2  ■  ■  =  5®  =  1-+'^  )  or 

<  r-  ?  -J-)  since  the  Integrand  Is  a  nerlodlc  function  of  0. 
Also,  the  part  of  the  Integrand  Uklch  Is  an  odd  function  of  9, 

integrates  out,  consequently  the  above  Integral  can  be 
written  as 

<«  51 

lr--hf  f  +  s)hB^smanSr.o^0,o,\ 

^  j  j  lYg  + 

O-^C  '  I 

Interchange  the  order  of  in^oEv-ation  -ac-w  •  , 

X - g- a cion  again  ana  make  the 

change  in  variable  z’  “^  +  «p  then 
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(E-4) 

Assuming  that  the  infinite  integrals  can  be  evaluated 
in  the  form  of  known  functions,  this  equation  is  much  easier 
to  evaluate  than  equation  [e-i]  since  the  range  of  <;  ,  which 
is  of  interest,  la  restricted.  The  difficulty  in  evaluating 
numerically  the  infinite  integral  of  equation  (E-1)  arises 
in  the  fact  that  the  integrand  is  a  peifodic  function  of  <=<  . 
The  infinite  Integrals  of  equation  (E-4)  can  be  evaluated  in 
the  form  of  Bessel  ftinctions.  This  is  done  by  writing 
sin  n(/)-  ^  )  and  cos  n(^  -  4’)  in  terms  of  each  angle  and 
noting  that  the  resulting  integral  is  in  the  form  of  either 
a  Fourier  cosine  or  sine  transform  whose  value  is  known. 

The  Bessel  function  equivalent  of  the  infinite  integrals 
which  arise  are 

O 


(E-5) 


i  COSl^fi  /  „  ,  X 


(E-6) 


CO 

(E-7) 


(2-8) 


The  ftinetlons  K^(no‘ )  and  Kj^(n<T  )  are  modified  Bessel 
functi.-ms^®  of  the  second  kind.  KQ(nir)  has  a  logarithmic 
singularity  at  cr  «  0  and  Kj  (ncr)  has  a  singularity  like 
^  functions  I^Cncr  )  and  I^(uo  )  are  modified 

Bessel  functions  of  the  first  kind.  r^,(n^)  has  the  value  one 
at  C7--  0  aiiJ.  I^(nO  >  n=£  the  value  zero.  Tne  functions 
L^(ncr)  and  Lj^(n<r  )  are  modified  Struve  functions;  and  both 
are  zero  at  O'  =0.  For  these  identities  to  hold,  tr  must  be 
positive . 

After  integrating  by  parts  and  using  the  above  values 
for  the  infinite  integrals,  the  following  equation  is 
obtained  for  the  integral  I 


'  ^8  +  'b  +  +  I  d 


(E-9) 


where 
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o 


-✓  ! 
o 


/  ^ 
o  i' 


(E-11) 
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(E-12) 
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J^=  cosnCH/f^^^)  -  |^S//J«(?[r,^«<r) df 

(E-13) 

The  only  difficulty  which  wi.il  arise  in  evaluating  these 
integrals  is  when  cr  ■  0  which  can  only  occur  when  y  »  1  and 

sill  ff.  Tue  integrands  of 


yf»u*  rui.  a-jl,  w 


2 


tan^do 

integrals  for  and  I_  have  logarithmic  singularities  at 

0  “  0  wlicn  X  =  1  which  can  easily  be  removed  by  a  <-.hange  in 

variable  as  discussed  in  section  II. 4.  The  integrands  of 

and  I .  are  of  the  indeterminate  form  0/0  and  can  easily  be 
d 

shown  to  have  the  following  value  at  0  -■  0. 

Integrand  of  I^  at  0  =»  0 


^■*o' 


It05<f)5ineti^ 

m 

•1  sin  4> 

-yJC^^  <r^ 

-  2n  cos  n  ^ 

Integrand  of  at  0  *  0 


(E-14) 


(E-l!)) 
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Of  interest  is  the  value  of  equation  (E-9),  for 
?  -  0.  For  this  case  the  Integrals  1^  and  I^,,  equation 
(E-10)  and  equation  (E-11)  respectively,  are  identically 

2ero  and  equation  I,  reduces  to  the  following  very  sijiiplified 
form. 


co^Bipcosctx^  KoLncr)cl(^  -  li)n&ifsin&npK,(na)d^ 


(E-16) 

The  integrand  of  the  first  integral  has  a  logarithmic 
singularity  at  x  «  1,  0^0  which  can  easily  be  removed. 

The  value  of  the  Integrand  of  the  second  integral  at  x  -  1, 

0  -  0  is  obtained  from  equation  (£-15)  as  tan/f^. 

The  integral  I2,  equation  (E-2),  will  now  be  considered. 
The  notation  given  by  equation  (E-3)  will  be  introduced  and 
making  the  change  in  variable  0  »  ll/2)(e  +«<)  equation  (E-2) 
becomes 


0C.OS  P 


(E-17) 

Interchanging  the  order  of  integration  and  making  the  change 
in  variable  ,/o  -  ^  ,  i  becomes 

-V  (  y 
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using  the  Bessel  function  equivalents  to  the  infinite 
integrals  this  equation  can  be  reduced  to  the  follot^ing 


le  +  If  +  Ig  +  Ih 


(E-18) 


where 


r  -  y 


a 

[cm3<^roS2Jifl(l  Jd 


(E-19) 


T  -V-T  l(Sin2ipSi 
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Um  / 


^  r  S 

[stm  (ps'm  p.n  <5)1  facosnf/O-CJ  j  ^  ] 
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(E-20) 


La  =  l(Coszpcos2np)^5inn(;Ko(i^<rJ  -  |"co5^  f  [jpfwcV  -teC««r;]  c/<5 

-t  *•  \ 


(E-21) 


J^=  -s^-y  j S  «  C  [l,  IWt'/  -  i-./  («cr;j  +  5  m  i?  K/  c/<5 

(E-22) 

Difficulty  may  arise  in  evaluating  these  integrals  when 

(T  «  0.  which  can  only  occur  when  x  *  1  and  0  «  0  at  the 

2 

same  time.  As  before  ^ «  Tanh —  ^  vdxen  x  -  1.  No 

can^o 

difficulty  is  encountered  in  evaluating  the  integral  for 
Ig  when  Cm  0  except  if  is  also  zero.  ‘Ihis  possibility  will 
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be  discussed  later.  The  integrand  the  integral  I  has  a 

O 

logarithmic  singularity  at  0  -  0  which  can  easily  be  removed. 
Both  integrals  and  have  integrands  of  indeterminate 
form  and  can  easily  be  shown  to  have  the  follovd.ng  value 
at  0  -  0. 


Integrand  of  at  0  "  0 


Integrand  of  Ij^  at  0  «  0 


(E-23) 


\ 

-  tan/So  sm 

(E-2A) 

The  value  of  I2  at  f  =  0  is  of  interest  especially 
because  of  the  form  of  Ig,  As  before  the  integral  If  is 
identically  zero  and  consequently  for  i!  "0,  I2  becomes 


2irh 

T(H*»/3o 


(ncr)  -  L.jfMcr/'] 


(E-25) 

A<f> 


If  in  addition  t:  «  1>  the  integrand  of  the  last  integral 
is  zero  at  0  =0,  but  the  integrand  of  the  first  integral 
has  B  singularity  at  0  -  0.  In  fact  it  can  be  shown  that  the 

f 

Integral  itself  is  infinite  as  In  00  .  Consider  the  first 

Integral  for  x  «  1 


